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i SYNOPSIS = 
_ method for determining the carrying capacity of restrained columns of “4 
constant cross } section 3 is described. _ The type of column failure considered is 


a 
"plane of bending - is also the plane of the thrust, applied moments, and the re- a 
_ straining moments. The method is applicable when the loads are applied to the a : 

columns with equal or unequal end eccentricities and can be used for cases 


: where the apr moments are not linear functions of the end rotations of = 


a hinged and for columns with equal applied end ‘moments and with equal ro- 
ean restraints at the ends. The availability of nomographs for the latter 
_ two cases reduces the effort involved in designing restrained columns toa <—l 


where such a _ design would be appropriate for - design office se. agg 


, pin- ~ended. They are con- 


translational and rotational restraints. The beams of an \n ordinary tier in, 


for example, serve primarily to loads to the columns. In doing so they 


— 
Note. —Discussion open until March 1, 1961. To extend the closing date one 
ep = a written request must be filed with the Executive Secretary, ASCE, This paper is part ie 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the re 
_ American Society of Civil Engineers, Vol. 86, No. EM 5, October, 1960. Seale pian ' 
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Sa 5 a _ leads to a solution of the true ultimate carrying capacity of such members. It [im 
— 
| 
Columns in framed structures are not, as a rule 
2 


2 also transmit bending moments because of the rigidity of the beam-to- column — 


ad's 
— These bending moments, together with t the axial loads in | the col- t. 


_ straining the end rotations of the column and thereby limit its deflections. ~ wo 
is thus seen that the restrained column must be studied in order toarrive ata 
more accurate estimate of the strength of columns ina continuous structure. 
Such knowledge is necessary inorder that further progress may be ma made in the 

application of plastic design to multi-story frames. 

. la Various aspects of the restrained column problem have been investigated by 

E. Chwalla,2 P. P. Bijlaard, G. 'P. Fisher, and G. Winter,? T. C. en 
F. | M. R. Horne, and and J. w. 6, and J. S. Ellis. 8 


FIG, 1 


% grossen Angriffshebeln," by E, Chwalla, Die Bautechnik, 10. Jahrgang, aot,  PP- 49 and 
“Eccentrically Loaded End-Restrained Columns,” by P, Bijlaard, Fi Fisher, 
and G, Winter, Transactions, ASCE, Vol, 120, 1955. 
Rational Simplification for the Buckling Length of Columns, by Kavanagh, 
Fa Proceedings of the Seventh Tech, Session, Column Research Council of the Engrg. Foun- ase 
a “Buckling of Columns with Equal and Unequal End Eccentricities and Equal and Un- 
equal End Restraints,” by P. P. Bijlaard, F Proceedings of the Second U. U.S. S. Natl. 
of Applied 
6 “The Behavior of Continuous Stanchions,” by J. F. Baker, M. R. Horne, and A. 
Roderick, Proceedings, Royal Soc., Vol.198,1949,. 
aa “The Elastic-Plastic Theory of Compression Members * by M. R, Horne, Journal 


of the Mechanics and Physics of Solids, Vol. 4,1956. 


“Plastic Behaviour of by J. S. Ellis, Journal of the Me- 


(the bending mo- 
h, the sig role o 
n 
™ if the ax mn_end 
t the colu 
ments a 
ie 


"The type of column that will considered is in Fig. The 
-* fully restrained against translation. The loads consist of the forces P “atl 
plied with eccentricities ea and ep: The ends of the column are also subject — 
_to restraining moments that are denoted by fa )and fp (6 to indicate 
a 7 their dependency upon end rotations 6' A and 6',. It will be assumed that the _ 
_ stress-strain diagram and the column cross section are specified so that mo- 
ment- ~curvature curves at several values of the average column stress may be 
constructed by the methods described elsewhere.9, » 10,1 11, 12 It will be further 
= assumed that the shearing force is small compared to the axial force i inacolumn. 
q Notation. — The letter symbols adopted for use in this paper are defined where 
first ‘appear, in the illustrations or in the text, and are arranged alpha- 


betically, convenience « of of reference, in the Appendix. 
A column in n equilibrium end moments and end rotations, as shown in > 
= 1, may be thought ofas a segment of acolumn deflection curve such as the 2 
= 19 in Fig. 2. A column deflection curve is defined as the | shape that | “= 
"compressed member - will take when the axial load and its maximum deflection 
Ym are specified. If unloading of the fibers stressed into the plastic range is 
as will be the case inthe present analysis, then the moment-curvature 


the curvature.12 For this condition the shape and wave length of a column de- 
' flection curve will be unique. It also follows, from the preceeding, that the col-— 


relationship will be such that for every moment there will be but one value of 


, that the entire column deflection curve may be constructed from a known se sola 


tified by the load P ‘and 86 (the ‘roation at y=  O) or rthe load Pand the. curvature 
ae at the point where the deflection is a maximum. In this analysis it was found 
- convenient to use 9% to identify the column deflection curve. . The column de- 


flection curve for a given value of 6, be constructed with the aid of the ap- 


scribed by S. Timsshedes. 11 ‘Computations were carried out foran 8WF31 rol-_ 
led steel column. When the stresses were entirely in the elastic range, checks _ 


;; with theoretically computed cx column deflection yn curves: showed that errors in the — 
wave lengthand in the maximum deflection were of the order of 0. 1%. The length 


of increments used in the numerical integration was 13. 13.88 in in. +» 0 or four times | 


the radius of gyration about the axis of bending. _ i 
9 “Untersuchungen uber Knickfestigkeit, arman, Fors ungsarbeiten, No. “a 


10 “Theorie des Anasermiitig Gedruckten Stabes aus Baustah . Chwalla, 
Bautechnik, 7. Jahrgang, 1934, pp. 161,173, and 180, 


* “Theory of Elastic Timoshenko, McGraw-Hill Co., New York, 


minsky, and L L.S. Beedle ASCE, Vol. 120, 1955. 


9 


Siege The present work systematizes the approach used by Chwalla and extends its ae 
to cases where the rotational restraints at the ends of the column ie. 
al 
> 
— — 


The solution to an actual ee column problem involves the quantities a 
ep, P, L, fa fp (6), the stress-strain diagram of the material, and 
shape of the cross section. The approach consists of fixing all but one of the 
variables and then examining a number of columns that satisfy the conditions 


of equilibrium so as to arrive at the one that yields the greatest value for . 
- open parameter. Since each of these columns may be considered as a segment 
of a column deflection curve, it is necessary, as a preliminary st step, to — = 
several column deflection curves for each value of 


GENERALCASE 
In ‘considering solution of the stability of Fig. the | will 
the variable and L, e,/ep, fp (@'), fq (6'), and P will be specified. Each 


= deflection curve available from the previous computations for the given 


load and the given cc column cross section will be examined to determine the lo- a 
Cation of the segment of length L which makes ea/ep equal to the prescribed 


- value. Fig. 2 illustrates the procedure for the location of the correct segment — 


on a column deflection curve. End A is located at a trial distance Xq fr from the 
end of acolumn deflection curve . and end B is ; found on the curve at a horizontal 
_ distance L from A. The angles 6, 4nd @p areknown at any sectionof the col- 
= umn deflection curve from the numerical integration used to determine the curve. 

The angies 6'A, @ ane and thé the eccentricities e,, € ep are given b by the the following 


vA 


fp 
B 
» 
ea and ep are positive when they are as shown 
a Eqs. 4 and 5 stem from the requirement that at Aand Bthe externalmoment _ 
to P ‘YA and P } YB reapectively, ‘as ca can be seen by a consideration of the column oo 7 
deflection curve (Fig. 2). By plotting ea/ep versus xX, as in Fig. 3, it is pos- a, 
sible to to arrive, at a value of xX, ¥ which locates the ‘Segment the 
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column curve representing. an equilibrium for the 


ratio of eccentricities. 


of column stability problems. *’ “The horizontal tangent gives the maximum 


SPECIAL CASE - 


When eg is equal to ep and fa (6) is is equal to fp (Fig. 5) the solution is 
greatly simplified. For a , est value of P anda given column deflection curve 

_ identified by @,, the slope and the bending moment P y, are known at any dis- ZZ 

_ tance from the point of maximum deflection. The three quantities are plotted - 

fora number of column deflection curves to give curves that are shown —— 

matically in the nomograph of Fig. 6. Each column deflection curve supplies 
- one pair of curves for the nomograph. In the upper portion the internal moment 
¥ is plotted against the slope 6 at the corresponding point of the column deflection 
curve. In the lower portion, the distance L/2 is plotted against @. The distance r. 
L/2 is measured from the low point on the column deflection curve to the point 
wheretheslopeis measured, 
this case, e, f(@), and P may be specified and the maximum column 
length consistent with equilibrium can then be found. _ Equilibrium at the ends © 

of the column in Fig. 5 requires that the internal moment of the column equal . 

restraining momentf (8) ‘minus the applied moment, P e. Int the upper oor 
of the nomograph of Fig. 6, curve A is drawn for the equation M = f (6) - 

a The intersection of this line with the M-6 curve of a column deflection curve 7 
_ gives the value of @ for which the equilibrium condition will be satisfied. A 
vertical line is carried down from each intersection to the corresponding L/2 

Be versus @ curve. This gives the half length of the column in equilibrium for an 
3 end rotation of @. By joining a number of points in the lower diagram, curve B 
is obtained. The point of horisontal tangency on this curve ‘gives | the ‘half - — 


SPECIAL CASE 


. 7 (a)) an equally convenient is Instead of 
= using @ @ as an abscissa it t is | necessary to plot 6" which corresponds to the ro- - 
tation at the restrained end. This is equal to. @ - y/L as can be seen from Fig. _ 
Pp 1 (b). The bending moment is determined as in Special Case-I while the length 
-- measured from the point where the column deflection curve has a zero de- 


flection and therefore a zero bending moment. This point of zero bending mo- a 


| corresponds to the hinged end of the column. The equilibrium of moments 


A nomograph for Special Case- was constructed for an 8WF31 rolled steel. 
column with a yield stress of 33, 000 )psian and a maximum of 9, 


| 
___ against the end rotation of the column deflection curve that gave rise to the value 
— 
‘ 
i 
4 
— 
— 
— 
24 — 
‘ia 
— 
: 


The residual stress pa identical to the one assumed by 
_ T. V. Galambos and R. L. Ketter.!13 This nomograph is shown in Fig. 8. The > 
length and moment ordinates were non-dimensionalized so that the nomograph 
could be used with little error for all wide flange sections.!3 Fig. 7of the Gal- 
-ambos-Ketter paper gives values for the collapse load of beam columns. These _ 


—N ve Moment. 


r, 


COLUMNS, 
nd to those of Special =Pe The construction 
golumns o. For this case M, = P ep. 4 
') is identically equal to zer 
— 


| 
in n Fig. 8 shows | how 1 Mo may be obtained when the slenderness ratio of the col- | 
umn is specified. . Table I shows a comparison of th the results obtained by this 
method with those of Galambos : and Ketter. 13° Asc can be seen, both h methods — 


the more general case when a restraining momentexists at the quan- 
tities ep, fp (@), and P may be taken as | the independent variables. . The maxi-— 
length with remains to be determined. The use of 


TABLE 1,— 


Mo 


from 


end eccentricities and end restraints. 


This is of general investigation entitled, “Welded Continuous 
Frames and Their Components” currently being carried out at Fritz Engineer- : 


ing ‘Lehigh University, under the direction of Lynn S._ Beedle 
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SOMPARON OF RESULTS 
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— , hs for this case follows the description given for the case of equal | —_ 
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of the Department of Civil Engineering. The author wishes to 

gratefully aid of Theodore V. Galambos, M. ASCE, in translating several 
references? ’ 10 and the helpful comments and suggestions on the — 7 


Reese by 4 Akira ‘Nitta, Theodore V. Galambos, and Lynn S. . Beedle. 


= Length of a segment of a column ‘curve. 


eng 0 a restrained co umn, 
length of t d col 
Slenderness 


Bending moment; 
= Applied moment at the _ of a member e-* 
= wes ‘Fully plastic moment value under ‘pure moment, 


Initial y yield moment value under pure ened: and inthe 


= Axial thrust toa 
ofa entire cross 


Eccentricities at which the axial thrust to the 7 


Radius gyration the axis. of bending; 


Distance along a column deflection the point of 


Displacement a a point on on the column deflection curve; re- 
‘Maximum displacement on a column deflection curve; 


between. a line j joining two points ona column deflection 
_ curve and the line of thrust, ayo 
of a a tangent to to the he column deflection 


is zero; 


the end point of the li line; i 


ee ‘Yield stress (assumed to be 33 ksi for AT steel); 
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STRESSES IN DOMES 


_ By P. Go dikas, Vg, _M. ASCE, and M. G. Salvadori,? F. 


SYNOPSIS 


a The bending stresses due to antisymmetrical wind pressures in a hemi- 
Spherical dome elastically built-in into a cylinder are determined by means 
of a Geckeler approach extended to non-symmetrical loads. The simple resco 


the asymptotic solution of Schwerin indicates the satisfactory accuracy ofthis 
approximate solution for all cases of practical interest. The methods used in 


= 


the solution of the present problem are easily extended to sphe rical sector shells 
5 stiffened by rings at the - boundary and to other types ¢ of non- seamen loads. 


» 


of the results allows their use by the designer. A numerical comparison with 


The membrane stresses in and circular cylinders due to 
antisymmetrical wind forces are well known, 3 or can easily be computed from : 
_ the equilibrium equations of membrane theory. 4 The corresponding displace- — 
ments are easily obtained by integration of the differential equations of 1 /mem-_ 


brane theory, but when the dome is built-in into the » cylinder, the membrane 


displacements. do not allow a complete statement of the continuity eaenees dome a 
ais 7 Note.—Discussion open until March 1, 1961. To extend the closing date one month, * 
ay written request must be filed with the Executive Secretary, ASCE. This paper is part | 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of a } 
American Society of Civil Engineers, Vol. 86, No. EM 5, October, 1960. 
 “* Prof, of Civ. Engrg., Columbia Univ., and Assoc., Paul Weidlinger, Cons. Engr., 
“Theory of Plates ”by ‘8. Timoshenko, Hill Book Co., New York, 
4 ‘bid, pp. 373- p. 384 
5 —— und Dynamik der Schalen,” 
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constants. “Thus, bending disturbances must be added at the intersection > 

dome and cylinder, requiring for their rigorous evaluation the solution of the 

_ bending differential equations of the cylinder and dome. This solution has been 

obtained in asymptotic form, and ‘requires lengthy calculations. The only 

‘problem evaluated by this ‘procedure known in literature isthe case of a hemi-_ 

spherical dome rigidly built-in around its edge. 

_ For the case of = serge: eawaee A loaded and supported domes, however, J. W. 

- Geckeler has shown’ that whenever the thickness to radius ratio is small (less 
4 than 1/50), satisfactory approximate solutions can be obtained by substituting — 7 
for the dome the tangent cone (or cylinder) and by evaluating the bending i 

stresses in the c cone. The case is tas reduced to tot the elementary 


yw p 
-PCOSOSING- 


‘FIG, Re FORCES" 


‘The same approach mz may — ‘to the -antisymmetrical case of 
- forces. The membrane displacements : at the intersection of cylinder and dome © 
: > ‘due to wind forces are first evaluated. The - complete differential equations for 


i “Theory of Plates an and Shells, by 8. Timoshenko, McGraw-Hill Book Co, a 


: eee auf dem Gebiete des Eisenbetons, Berlin, 1921, 


“Beams on Elastic baad The Univ. of Michigan 
163. 
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© *Uber Spannungen in Symmetrisch und Unsymmetrisch Belasteten Kugeischaien 

(Kuppeln) insbesondere bei Belastung durch Winddruck,” by E. Schwerin, Armierter 

y J.W.Geckeler, Fors- 

— 


the anti-symmetrical in the cylinder ; are 
then : solved, and their rigid body n membrane and bending components analyzed. a 
_ ‘The bending displacements are expressed in terms of the transverse shear and ~ 
~ bending moments at the cylinder boundary. Continuity between cylinder and 
- dome is finally: restored by means of boundary shears and moments, using cy- 


lindrical ahead for the dome as well as for the cylinder, a and the values of the a 


ave 


shears moments s at the dome and cylinder are thus ex 
uated in finite terms. The approximate solution obtained by this procedure is 
7 accurate for small v values of the rica to radius r ratio. A comparison with 
the fairly thick dome e (n/a - - > 1/40) used! in the example of E. Schwerin8 indicates é 


— 

— 

— 
— 
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October October, 1 960 EM 
‘that the accuracy of the present nt solution is is , sufficient in ficient in all practical ca cases e1 en- 

_ countered in construction, where thickness to radius ratios are of the order a 
1A 00 or less. The method used to solve the present problem can be extended 
to evaluate bending stresses in spherical domes under other types of antisym- 

metrical loads, with opening angles 7/4< taking into — the 


? effect of a stiffening ring at the dome base. ea 


— 
, anormal wind 
p cos 
t 6 = 0, is assumed to 
on a semi-infinite cylinder. The wind pressure p is one- -half the pressure re- 


quired by codes, since Eq. 1 indicates that a suction of - ries p acts al at "* = bal 


x/a 
and the Z of aie’ 1, under the assumption « of a traction : 


(the. zero indicates. the value of a a quantity at give ive the mem- 


stresses: 
MEMBRANE DISPLACEMENTS DUE TOV WIND. IN CYLINDRICAL 
‘The membrane displacements (Fig. 2) are related tothe membrane stresses 7 


by the differential equations: 7 
(1 - 


_ 2a(1-v) 
c 


7 


which hg is th 
of Plates Shells, by S, McGraw- -aw-Hill Bool Book Co., New York, 


«Statik der Schalen,” by w. Fidgge, ‘Springer, Berlin, (1987, P. 
Bas 
16 


| 
a 
4 
integrated for X =¥ 
less boundary (§ = 0) 
— 
wa 
: 


Assuming, in view of the load of 2 
=Ngcos@....... 
‘a 


Neo 


= 


we =weosé . 


"where ‘the amplitudes functions of 5, 6, and 7 


(Ba) 


Ai 
he mem 


the cylinder boundary ( (€ = 0) Eqs. | and 11 give: 


2P sine, Wo = % Wo 


a MEMBRANE STRESSES DUE iE TO WIND WIND IN IN SPHERICAL DOME 


6 


with the symbols : of Fig. 1, a normal wind pressu 


aie * 


— 
| 
— 
7 
— = 4) cos 


N 


= 


. The subscript zero indicates the 


a Se) 
(18¢)_ 


— 


4 

te 


value of a at the boundary» 


i. = 1/2, _ The membrane displacements in the dome are related to the mem- 


= 


gl 


where the amplitudes . 
of Eqs. 17 with the condition that 


w=B sin ¢ log (1 + cos +. 


y 


=. 


, are functions of ¢ only, the integration 


‘the displacements be finite at the top c of the 


(20b) 


— Ng ope Ngo § 
— = 


‘WIND 


At the dome (9 = 


the are: 


(22) 


Once and are known from membrane anal- 
ysis, a general state of stress must be superimposed on the unloaded cylinder _ 
ai. in order to satisfy continuity boundary conditions. The bending stiffness of the 
on . must be considered in order to obtain meaningful results, but the bending» 


stresses vanish rapidly away the boundary, so the 


corresponding homogeneous differential of equilibrium: a 


stituted in 


<a 
and the determinant of the resulting. algebraic equations in U, V, W is set alia i. 
two 4th 


zero, the 8th degree equation for splits into ‘the 


— 
a 
— — 
Qa) 
ne? 22) 
— 
(23a) 
— 
i 
(252) 
— 


he 


mn approximatio 


266-08, Tech. 


cctticients tof Eq. 26a are polynomials in the th 
q. 25a are polynomials in the thickne ratio parameter: 
2. 
— 
Report Ring-Stiffened Cylindrical Shells,” Contract No, NR 
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v_ 


am 
The displacements corresponding to 2 decrease away ; 


>: 


ta can n be ‘inal 4that the bending disturbances atone of the cylinder ends a are > 
not felt at the other, so that only the solutions corresponding to the roots ug, 4 
need be considered to evaluate bending disturbances in the neighborhood of = 
Eqs. 23, 27, and 29, these displacements become: 


cos + Ag sin Be § Je 


in which . Ai, Ag are constants integration. 

At the cylinder top ( f= 0): 


2 


oe © 


— 
that the correspon ratios ares 
4 
— 
large and the cylinder is thi 
— 
— 
— 
— 
9 - Ay) cos Be — 
= . = Be - Ay) Be 
G2 
—........ (2) 


- " 


— 
FF ™~ | 
ag 


3 


and bending moment M¢ per unit of length 
ments (with the notation and 


sign convention of Fig. 3 3): wits 


ne 


¥ 
Be Ais sin « 


7 


and: at the e cylinder boundary "a = 0 have the e values: : 


- 1+ cos é . 


. 
Solving Eqs. 35 for the constants Aj, Ao, and substituting in Eqs. the bound- 
ary displacements in terms of the boundary tractions become: at 


ER, 
Eh, Mo og sin @ 


“Statik und 


=shear, Még = twisting moment) 
Substituting the displacements of Eqs. 31 in Eqs. 33 band c, and neglecting the 
ever k, X1), S¢ 
i 
——— < der Schalen,” by W. Fligge, Springer, Berlin, 1957, p.147, [ii 


= 


moments in become: 


“Po§ 
e sin n Be . 
Non- Sending Terms of General Solution. displacements 


+ 
‘The ratios — ona -V/W are found by substitution of Eas. 38 in the original 


differential equations, 13 and by setting equal to zero the coefficient of each | 
power of Neglecting ke in comparison with unity, the non-bending displace - 


[Wo + Wy + + Wy 
‘Where the Wo, ... 


a ‘rigid body translation in the direction 6=T,. 


rr The terms in W,: 


| 


in Eqs. 34 the shears and 
— 
— 
— 
— oa 
— 
i 
— 
— 
represent 2 ,/a around the axis - 7/2 at 
Substitution of the terms in W 2 in the membrane stress equations, Eqs. 


(4) of the in into 6 gives: 


he 
77 W3 cos 6; 

At the top of the cylinder (§ = 


the traction | N the cylinder base, E=6€ 


‘Thus, from Eq. 39, 44, and 48, the non non- -bending at the 
boundary are in terms of body motions and stresses 


Solving Eq. 46 for ny 


cos 6. 


of 


«54 GENERAL ANALYSIS OF THE UNLOADED HEMISHPERICAL DOME 
Bending T Terms of General Solution. _The ba basic simplifying assumption of 
a Geckeler applied to the present problem implies that the bending ter terms of the _ 
gene eral solution for the ti dome displacements c can be approximated — 


nal stress, varying sinus- 
Nog = Wa cos 83 Nog, 
— 
Né@ is constant in and at the free boundary produces an unbalanced 
resultant and an unbalanced moment equilibrated by the normal traction and 
A 5 a 7 


ag 


by the bending displacements of the cylinder tangent to the dome at its equator, a 


ince bending disturbances are limited to a narrow band at its base. _ — 
_ The ag tog bending displacements near the dome base are therefore 


| 
cos Bg X + Ag sin X 
cos bg X (AL + Ag) sinB,X 


ike 


constants Ag ar are. by the equations corresponding to 35 and 
‘the at the dome base (X= 0) become: 


== aS, 


ed 


9 


Non- -Bending Terms of General Solution. non- of the 
general solution for the hemisphere are only rigid body displacements, repre- : 


=- cos ¢ ‘cos 


Ly 
cos cos 


= 


‘A, + cos B X A, - Ao) sin 8.X| cos 0. .(51a) 
— 
— 
and 
— 
— 
— 
= 
— 
— 
(2) a rigid of amplitude around the axis = 7/22 


Thus the ‘non- “pending displacements at the dome base at X = =0 0 are given by: 


Indicating by the rotations of the cylinder generatrix and the meri- 
-dional tangent to the dome at the cylinder and dome edge, respectively, the 
‘compatibility of at the common edge requires that: 
= (Woe = (eet (Wo = ( \ 


Assuming Ve a of -0 axis | the - axis, ont 

positive for a rotation of the = axis tor towards the = 0 axis 


13, 36, 49, and the ‘edge displacements of the hemisphere given by 22, 52, 
55, four equations are obtained in the eight boundary tractions at the common — ; 
edge ofcylinderanddome, 


7 However, between these the following relationships subsist: 


H 


z 


log Nopg =} + bé 


— — 
— 
— al iii 
bie — 
— 
= Ng 
und- 
— 
ze A _ However, the terms depending on the Np in the expression for Sog and Mog tee 
can be neglected since they are of the order of 1 agninet the remaining 


the take the 


(1. 


as 2. D, . .(60c) 
The constants D, and Do 
P : of the cylinder with respect to the dome. Eliminating one and D2, dropping» 
members of the order of unity in comparison with 62, and and letting: 


and M,, become: 


1+A5+ 


1+v 


kness = 1, 


thick, that for : a. edge 


i 

— 
— 

When dome and cylinder fave the — 


7 The shears and moments at any point in the atin and shell may be computed | 


the cylindrical equations, Eqs. 37, in which and or Bg and X respec- 
tively, ar are used used with the amplitudes Mog of) of Eqs. s. 62. 


_ computes for a dome built-in at the edge with the following aeons 
= 100 kg per sqm,a=10m,h=0.25m,andV=0.25, 
- vist Table 1 gives the values of the boundary shear and moment evaluated by 
‘Schwerin and the alues by means of Eqs. 65. 


v 


— shear Asymptotic solu- 
and nent | 


tion, in Kg per cm tion, in Kg per 


Although the shell considered thick (h/a = 1/40) w pres- 


ent day standards of ‘construction, the accuracy of the present solution is sat- 


_isfactory for engineering purposes. The present ‘solution, moreover, seen 
a Since the accuracy of the present solution improves very rapidly with a de- 
creasing thickness ratio, the present method of solution is proved acceptable 7 


in all cases of significance, , since hf i is than 1/100 in 


| 
— 
— 
— 
TABLE 1,—BOUNDRY SHEAR AND MOMENT 
» 
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. Pister, M. ASCE and L. L. Williams? 


In a recent note te E. Reissner extended the ‘concept of ‘the classical Winkler - 
_ Zimmerman foundation to permit the introduction of a differential shear stiff- 
“ness parameter | - in addition to the foundation modulus. The resulting system of 
differential equations is of sixth order. In this paper the system of equations © 
- is modified to incorporate a different interface condition, and the significance nl 
the shear stiffness parameter is is examined. 
AS an example, the solution for a concentrated line load on an infinite plate ay 
‘supported by a Voigt foundation is along with — curves 


The analysis. of beams or r plates. ona foundation re rests. upon as- 
sumptions concerning the behavior of the plate-foundation These as- 
"sumptions involve: (1) description of the foundation, that is, , whether a solid | 
oma or an 1 idealized set of independent elements exists (2) conditions to | 


4 a written request must be filed with the Executive Secretary, ASCE. This paper is pert. 
the copyrighted Journal of the Division, of the 


, — 1 assoc. Prof., Div. of Struct. Engrg. and Struct. Mechanics, Univ. of California, Ber- 
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rof,, Guggenheim Aeron 


— 

— 
— 
— iia 

— 
— 
— 


at! the plate- -foundation interface (full y of or orfric- 

_ tionless) and (3) system of equations by which the behavior of the plate is de- 

i 4 fined, that is, the Poisson-Kirchhoff theory neglecting normaland shearing de- 

formation, the Reissner theory retaining these effects, or consideration of the = 
plate as a solid continuum. In view of the number of mers combinations of 
these assumptions, many theories of behavior have appeared. No attempt will 

_ be made here to present a bibliography but typical treatments may be found in 4 

Hetenyi, Burmister, Naghdi and Rowley, >. and Pickett.6 A common feature of 

the work citedso far has been the assumption that both plate and foundation be- f 
have as linear elastic bodies. Recent papers by Freudenthal and Lorsch, 7 and 
Hoskin and LeeS have generalized the problem in the sense that a Poisson- 


It need hardly be emphasized that in the treatment of the foundation a com-_ 
- promise between the simplicity of the Winkler- -Zimmerman set of independent — 


springs and the behavior of the solid continuum must be reached. Recently, E. — 


Reissner? generalized the independent spring model by including a differential a 
shear stiffness parameter in addition to the foundation modulus, _ at the same 
time retaining mathematical simplicity in the model of the foundation. The plate ™~ 
deflection was governed by the Poisson-Kirchhoff theory and the foundation in- 
terface was assumed to experience no displacement in the horizontal plane. In _ 
the present paper Reissner’ s work is modified to obtain continuity of displace- -— 
- ment at the interface. This has the feature of emphasizing the role of the ratio. Be 
of foundation to plate thickness. In addition, the physical significance of the 
Shear stiffness parameter is examined. of normal = 
mation in the plate, following Naghdi and Rowley* and Frederick, *Y will be re- 


__ Notation. —The letter symbols adopted for use inthis paper are defined Winte oe 
_ they first appear, in the illustrations or in the text, and are a arranged ional = 
ically, | for convenience of — in the Appendix. 


Foundation. — 
_ is assumed to be capable of transmitting normal stress Oz ‘and shearing stres- — 
Txz al and ad Tyz- The remaining stress com components, Ox, Oy, T xy, are assumed 


Bre. 3 “Beams on 1 Elastic Foundation, by Hetenyi, Univ. of Michigan Press, Ann - 

ig 4 “The Theory of Stresses and Displacements in Layered Systems and Applications — 


Naghdi and J. C. Rowley, First Midwestern Conf, in Solid Univ. 
of Iilinois, Urbana, 
_ 6 “Deflections, Moments and Reactive Pressures for Concrete Pavements, " by G 
Pickett, et al., Kansas State College Bulletin No.65,1951, 0 
rit. “The e Infinite Elastic Beam o ona sage Viscolelastic Foundation,” by A, M. 
an, 8 “ Analysis of Flexible Surfaces over Subgrades with Viscoelastic Material Behavior 
a B. C. Hoskin and E. H. Lee, Proceedings, ASCE, Vol. 85, No . EM 4, October, 1959 . 
9 “Deflections of Plates on Viscoelastic Foundation,” by E, Reissner, Transactions : 
_ 10 “Thick Rectangular Plates on an Elastic Foundation, ” by D. Frederick, Transac- 
tions, ASCE, Vol. 122, 1957, pp. 1069-1087, 
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iq 
Kirchholl beam or plate resting On a SeLOl independent viscoelaSlic Springs nas 
— 
q 
1 
le 
— 
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BENDING OF 


« 
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La 
a4: 
, where partial differentiation of the itn is indicated in the usual m 

stress equilibrium e equations, in view of the above assumptions, th 


_ 


7 
Plate 


— 
— Plate.—The | vlecting the — 
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satisfies the equ equation, 


BOUNDARY AND CONTINUITY CONDITIONS 


‘The interface conditions, Reissner’ s =y assume continuity of 
displacement between plate and foundation. A final condition imposed upon the 
‘foundation is the requirement that the on the cylindrical 
‘Choosing the plate deflection W and pressure p as dependent 
ables, the basic system of equations can be treated as follows: q 
7 


_ 2 show that at T xz» Tyz a are functions of x, y only laid that 


* - where the condition Cz (H) = p has been used and Q=Txz,x + T yz,y - Inte- 


grating Eq. to the conditions of Eq. 4a, for the dis- 
if H 2 


pet 
condition Eq. 4b, )=¥ w y, Eq. Be gives 


deformation. Accordingly, the plate deflection W (x, f 

Where q and p are the 
i terface pressure, res 
j. The usual Kirchhoff plate boundary con a 
With Eq. 3. For the foundation the followi 

a 


‘BENDING OF PLATES 


4 + ) 
Woy 


Eliminating the shearing stresses sesee trom Eq. 8, combining with Eq. 7 and id sim- 


‘Eqs. and 9 constitute the basic system of the dependent vari- 
= ables W, p. In addition to the two boundary conditions involving W, we have the _ 
bw requirement that the shearing stress on the cylindrical boundary of the foun- 
dation must v vanish. From 8a, and 8b this is be 


GE » m 


In Eq. 10 n refers to the outward normal to tl the teundintien boundary. With the — 


exception of the continuity conditions of Eq. 4b and the presence of the term — a 

(1+ +a)! Eq. 9, ‘the | basic system of equations and boundary « conditions are 


= is set equal to zero, s 
equations are recovered. However, while he useda single parameter” topermit 


_ inclusion of viscoelastic behavior of the foundation, the present work will incor= - 
‘porate viscoelastic effects through the well known correspondence principle, ! 7 
connecting elastic and viscoelastic boundary value problems. Briefly, it can 
be shown for a class of problems that a correspondence exists between the 
stress anddisplacement fields andthe Laplace transforms (withrespect 
to time) of the viscoelastic stress and displacement fields, provided that the — 
- elastic coefficients are replaced by the Laplace transforms of the viscoelastic © 


PHYSICAL SIGNIFICANCE OF FOUNDATION PARAMETERS 
As pointed out by Reissner,* the present theory permits description of the 

foundation using two elastic parameters; the usual axial modulus a a Rico 


_ appearing as the term in Eq. 9. Although the physical significance of the _ 


by F. Eirich, Ed., , Vol. L , Academic Press, New York, 1956, Chap- 


Which relates normal pressure | to deflection, and “differential shear 


— 
| 
— 
— 
— 
— 8 
4 


fa 
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well known, it is of interest to exam examine _— term. 
It will be convenient shear stiffness parameter 
Rewriting Eq. 9 using Eqs. 1 


6-8 


he following the significance of ter rm Consider 


the foundation layer by a ap 


g = 
q 
— 
— | 
4 | 

ed directly tothe inter-_ = 


BENDING OF PLATES 37 
‘The solution the condition of finite deflection at and 


’ >0 er 


a factor of e, , hence, it is a measure of the differential shear stiffness of the - 
- foundation. For the Winkler foundation, C = 0, and the discontinuity s| shown in 
“ the figure results. Accordingly, in the present theory k, C Hare the elastic» 
or viscoelastic foundation parameters or a combination of the two. _ —, 


k , viscoelastic 
CH, constant,assume 


to base the value of k on results of a plate ohn test, rather than to base the ~ 
value on its definition | k= ‘E/H. A similar test could be performed to obtaina — 
In this manner the foundation can be 


described by two experimentally determined parameters. 


ILLUSTRATIVE EXAMPLE 
this section the relatively simple solution of the basic system 3 
and 13 will be indicated for an infinite plate supporting a load town in Fg. 


distributed with respect to the y-axis and independent of y as shown in Fig. 3. 
The basic system of equations for the elastic case are Eqs. 3 and 13 where 
We W (x, y, t), p =p (x, y, t), and q = q (x, Ys t). Denote the ‘Laplace transform 


with respect to time of W (x, y, t) by oe ip aya a 
and the Fourier cosine wes of W (x, y, t) with respect to: x by 


ax dx 


y, t) cos 


— 
Eq. 16 is shown in Fig. 2. is the dimension- 
less distance with respect to H in whicn the foundation deflection diminishes by _ 
— 
— 


operators are ‘replaced by their Laplace transforms with 
to time. Accordingly, for-the one-dimensional case under consideration, and 
for “at re rest” transform n equations be- 


D- 


where it. ‘may be no that the plate has been in and 
hence no transform is indicated. a: of the Fourier me transform 


4 
7 


a” W (a, 8) -B (a, 8) 


2 


(s) H2 |= 


ining transform of the plate. deflection 


> 
“= the definition of A is 


- To be specific, assume a step pressure applied along the y-axis together ak 
with a Voigt axial foundation modulus, which is, furthermore, assumed incom-_ 


1+ 


_ The correspondence principle for viscoelasticity states that the same differ- 

— 

— 

SESS: 
— 


‘BENDING 0." PLATES 
and viscoelastic. modulus is 


where the elastic of the m= E/E, and th retardation time 


is 


ak. ot 


- # 


Rig 


and the inversion the along X 


develops that it is more convenient change the slightly in order 


to exhibit the physical characteristics of the sobation more To this end, 


od 


1! 


D(a 


(a 

— 

the condition of incompressibility, v = 1/2, has been in the 

psuedo- ie rigidity Dg. Returning nor now to the definition of A La(z)], ad 


the fact t that =1/9 and 


— 
ia 
ii 
| 
- exp +mA(a) | = 
— 
— 
— 
= D_ at = a’ = 7 
— 
—- 
4 a q 


ny A also © approaches ; a 2 constant, 


( = 


cos (V2 
cos xz)dz 


T) 


= 
_ 


2 


va 

w(0,t/7; 0) 

in which it is noted that the initial velocity is infinite, and t the: fora 
a relatively stiff plate is slower — the opposite situation ge the damping factor oo 


3 


The first of these pertains 
Several limit cases ma — 
namely unity. Similarly, >load may — 
—— at an infinite time under the lo | 
, So that the limit deflection 
— 
rest and may be computed easily. 
The velocity of deformation is also of Interest 
4 


(1/46 )(0¢70) © 


' 


FIG, TIME DEPENDENCE OF F VISCOELASTIC I DEFLECTION, W=W (x, t/7; 


— fo of | | | 
4,—VISCOELASTII 
- iii 


| ymotote 


= 


i 


i wie to present some ‘idea of the variation of the deflection and moment > 


& with time and rigidity ratio, these functions have been | computed and plotted in — 

Figs. 4 through 7. It may be noted that these graphs have been plotted for the Z 

7 A special case of h/H = 0, rather than the more general case given in Figs. 36, ; 
39, and 42. However, by observing the difference between the limit values, as _ 
a functicn. of h/H andthose where h/H = 0, an estimate of the effect of this para- 
meter, which represents the shear bond between the foundation and plate, may 


= 


From a physical standpoint, probably the maindeficiency between the theory 
and the actual behavior of the foundation is a direct result of assuming that the — 
-plane! foundation stresses (o,, Oy, are negligible, whether one con- 
-siders a slip or continuous intertace. It is this assumption that leads directly — 
to the uifear variation of normal stress 0, with respectt to vertical distance, ir-_ 
7 e of foundation depth. This result : appears to be. ‘intuitively contra- 


-dictory R ith a St. Venant type variation that might otherwise be expected. On 


the othef hand it is encouraging to find that even this set of assumptions pre- of 


“dicts an effect depending upon the interface boundary conditions, even though | 
in pract cal civil engineering applications one is tempted to set h/H — 0. eae 
_ It is Rlso interesting to observe that a formal elimination between Eqs. 3 

and 13 an a of Eq. 1 13 to 


ae 
where ldngths have been r non- -dimensionalized on H in the operator. While 

it would|appear | that the series converges rapidly asa consequence of C= 1/3, 

_ further investigation indicates that a truncation of the series after the first 
term, giving what amounts to a beam-column term contributed by the founda-— 

- tion, is insufficient for predicting the actual behavior because the non- -dimen-_ 
9 sionald rivatives in the higher terms become increasingly large. On the other | 
hand it May be possible to use such an approach ComaS wee particularly 

a Ih conclusion it is felt that | useful engineering results can be obtained by a 


—  ° ° ° ° ° ° ° 
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APPENDIX. 


rectangular “components of displacement vector i in 


= rectangular components of stress tensor in foundation; 


= plate and foundation thicknesses, 


= interface pr pressure ‘e and transverse load on on n plate; a 
axial foundation modulus; __ a 


pseudo- -foundation ‘rigidity; an and 
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enency ANALOGS OF VIBRATING BEAMS 


By Frederick L. Ryder} ASCE 


_ For the passive-element electrical analog of a finite increment of a beam © 
_ vibrating in flexure, improvements over the conventional form of the circuit 
are found by making explicit use of the energy as an analogous parameter in oe 
the beam and in the electrical circuit. For a given beam it is found that these - 
_ improvements: lead to a drastic reduction in the theoretical analog error and 
also in the number of electrical components required per beam increment. . It a 
generally be expected that for a given allowable overall error 
_ the number of required beam increments may be reduced, while the cost of 
analog components may be reduced both because fewer components are used — 
_and also because, in consequence of the smaller number of components, greater 
_ imperfections in the components may be tolerated. Alternatively, for a given 
installation, the complexity of problems ‘which may be with a 


INTRODUCTION 


TRODUCTION, 


i Analogs for | the Simulation of the dynamic behavior of beams by means me © 
_ passive electrical elements are well known. 2,3 Ss is generally assumed that — 


the beam deflections are small, that the damping is negligible, and that the — _ 


electrical components (which may cons — of inductances, and 


7 _ _Note,—Discussion open until March 1, 1961, To extend the closing date one month, : 

7 a written request must be filed with the Executive Secretary, ASCE, This paper is part 

_ of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the _ 
> | “American Society of Civil Engineers, Vol, 86, No. EM 5, October, 1960, Sa 
4 ak Member, Scientific Research Staff, Republic Aviation Corp., Farmingdale, N. Y. _ 

7 — of Normal Modes and Natural Frequencies of Ship Hulls by Means of 
; the Electrical Analog,” by E, Kapiloff, David Taylor Model Basin Report 742, July, 1954, 


2 « 

e 

ee. “An Equivalent Circuit for a Vibrating Beam which Includes Shear Motions,” ” by H. L 
-M, Trent, The Journal of the Acoustical Soc, of America, May, ; 


‘ — 
— 
ud 
= 
| 
a — 


2. The equations" ‘of deflection compatibility, expressed | in terms of | the 
_ forces, moments, and elastic constants, are simulated by the equations gov-— 


.. —= the compatibility of the time- -integrals of voltage, obtained with the aid 


the effects of nee | turns ratios on the voltage relationships). 

ts ‘Since the equations s of equilibrium and of deflection compatibility s serve to 
define the mechanical behavior, and Kirchhoff’s two laws serve to define the 

_ electrical behavior, the foregoing requirements assure that the various electri-- 
eal and mechanical parameters are analogous. However, there exists an al- 
.% ternative | lesser-known approach in which the simulation of equilibrium rela- 


Bans is not taken into account explicitly, but results ey from 


the following two requirements: 


©) TRANSFORMER 


IG. 1 1. ELECTRICAL ‘FUNDAMENTALS 


‘ses, mass moments of ‘inertia, and is the Same ai as the capacitive 
energy of the electrical system (written in terms of capacitances and voltages). 
4, The strain energy of the mechanical system (written in terms of de- — — 
flections and elastic properties) is the same as the inductive energy of the 


electrical | system (written in terms of wal time- -integrals of — and the 
4 
_ The familiar Lagrange equations of the mechanical system, based on the 
kinetic and strain energies, are known to. define the deflections, or time- 
- integrals of the velocities; hence the identical equations in the electrical Sys- 
tem must similarly define the time-integrals of the voltages, so that voltages — a 
and velocities are analogous. _ Further, any force (or moment) is the il 
derivative, with respect to its associated velocity component, of the power (or 
_ time- rate of energy) in the mechanical system, while current is the partial - 
power, with respect to its associated voltage difference, in the — a 
electrical system, that currents’ are analogous to forces moments. 
- Hence the summations of current must automatically simulate the equations _ 
force and moment equilibrium, although these are not 


in n steps | 3 and 4 4. This will be illustrated —— 
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ideal. An analog circuit is then found which, in effect, meets 
1. The equations of force and moment Drium, Tertia effects, 4 
_ in the structure are simulated by Kirchhoff’s first law, relating to current j ees 
| 4 
— 

| 
| 

4 


: two methods to the simple case of a flexurally vibrating beam, that — 
€ advantages may result from the use of the energy me method as as against the better- —_ 
equilibrium- 


In preparation for the analysis of the electrical will be well to 
state certain fundamental relationships pertaining to the ideal gprs: ll 
inductance, and transformer, which are subjected to certain currents 


so that 


= 


st for the transformer, 2 ie and Ng are the respective numbers ot turns, all 

considered positive if positive current produces flux in a prescribed direction “a 
in the core. od the law of electro- -magnetic induction: ‘0 on 


Also, the summation of the ampere turns of both windings ‘must be zero, ; 


since the net magnetomotive force required tc sustain the magnetic flux in the 


i i2 = 
“The total power denoted in the is therefore i 


ITE INCREMENT | oF om 


- o n+l... The length of the increment is Ax, and it is assumed that the beam is 7 
uniform throughout this length. The mass per unit length is yp, the elastic | 
: - modulus is E, and the moment of inertia of the cross section about the neutral — 


- flexural axis is I. Purely for convenience, it will further be ascuned that the 
two neighboring increments are identical to the one shown. _ 


— _ For brevity the method of steps 1 and 2 will be termed the equilibrium- es, 
_——" = compatibility method, while the method of steps 3 and 4 will be called the pee 
— 
a 
1 
ee 7 “4 where the dot denotes differentiation with respect to time t. The capacitive a a 
— 
—_— 
eo 
— 
is conventional, the beam is divided into a number of incremental 
4 ay 


accumulated ma mass - to n + namely Ax i is 
concentrated at n, and an identical mass atn + 1. In order to take 


the section, we consider the mass moment of inertia per unit length, namely ~ | 
pb p2 (where p is the radius of gyration of the cross section), and refer the 

_ corresponding value for the whole increment, namely yu | 2 Ax, to the center of 
the increment. The external load forces Fp, and Fy+1 are » assumed to be ap- 
_ plied as indicated. The shear compliance (deflection per unit shear force) for 
the entire increment is Ax/(K A G), where A is the cross-sectional area, G is 
the shear modulus, and K A is the fictitious cross - -sectional (always 
smaller than A) which, if uniformly loaded by shear force, would experience 


ft the same vertical shear stress as the actual maximum, which occurs at the 

neutral axis; this compliance is assumed concentrated at the center of the 
_ increment. The bending compliance (angular deflection for unit uniform 1 mo- 
- ment a at all points of the increment) is eg I), which, together with contri-_ 


r 
4 


2. _—INCREMENTAL LENGTH OF BE. 


and mechanical parameters are for ‘convenience taken as unity. y. The upper 
Bee's of the circuit is devoted, essentially, to the simulation of shear effects, 
- and the lower to moment effects. The capacitances are given the values of 
ening masses and mass moments of inertia, and the inductances c cor- 
to compliances. The voltage at station n is the translational velocity 

and correspondingly for station n+ 1. The voltage at station n +} is the. 


rotational velocity +i, Where 7, is the slope of the neutral plane due 


to bending this be constant the increment because 
_ bending compliance is concentrated at the ends. The supplied currents F, and 


n +1 Simulate the corresponding load forces. 
"The currents V, and and Vp + Simulate the shears immediately left and 


respectively, of n, the shear being positive when ap- 
plied in the upward direction by 3 any section to its neighboring section on the © 
right. The currents Mn and My ; 7 Simulate the moments at stations n and o 
nel, , respectively; moment is caneheened positive when applied ina counter- 
clockwise direction by any section to its neighbor on the right. The voltage > 


— 

— 

‘4 

£4 

™ An electrical analog for this beam increment has already been given* and © a 

: | 

- 
, — 
_ 4 
\4 


- across one side of the transformer is Fo an and across the other side y, 
_ Ax, so that the magnitude of the turns ratio is Ax, and its polarity can be de- 


: duced with the aid of | the plus and minus signs onthe upper transformer winding. | 
The validity of analogy will first be proved by the -equilibrium- 
compatibility method, in a manner which is more suited for the He sewed a 


ering shear in the beam, we have: in 


n+ n 
which by Eq. Lis is the of current summation at n in the 


FIG. 3. ELECTRICAL ANALOG OF BEAM INCREMENT 


: a ing moment equilibrium in n the beam, we e have: oe a 


41 = Mn- Vp Ax- dx 
phil To determine the analogous electrical equation we we first strate the cur- 
rent i with the aid of Eq. Tas 


es that Eq. 10 ) exactly denotes the current summation at stationn +}. | 


— 7 
— * 
— 
— 
— — 
&g 
— 
— 


50 
As to the equations of deflection ‘capability, the he slope lis 1/ax times 


the difference of deflection, due to bending effects omen, between ‘the two ends 
| the increment, making due allowance for the shear deflection, there results: 


Yn * 

sa 2 K A G | 

The corresponding electrical equation, obtained d by « considering the he voltage = 


inst "In EEG n+} 
where use has been made of Eq. . 3; this exactly fits the preceding lait 


(assuming suitable constants of integration). Finally, considering the change 
_ of slope between the increment in Fig. 2 and that at the left, one can write 


nE 


ve 4 
and similarly for the change of slope suiaetiee the increment shown in the fig- 
ureandits neighbor atthe right. 


Since all equilibrium and compatibility equations are simulated t by corre- a 


Sponding equations int the electrical analog, the latter is valid. Consequently, 
if the circuit is connected to other such circuits in the same manner as the | 
beam increment is connected to other such increments, and if the boundary 
conditions of force (and moment) and deflection are simulated by correspond- — 
ing conditions of current and voltage respectively, then the overall 


The kinetic energy of the beam increment, wee for convenience to extend 
om just to the left of station n to a point just to the left of station n + 1, is 


Ax + 


the foregoing ta the deflection | difference due to shear 
alone is exactly the time-integral of the voltage difference across the induct- a 
ance —Ax/(K A G), then by Eq. 5 the strain energy is exactly simulated by i in- 
“ductive energy if the analog. The analogy is ) therefore valid, with « currents x 


- 
— 
= 
= 
7 _ difference between stations nandn+1duetoshearaloneis 
— 
— 
| 


BEAMS 
‘simulating | forces and moments because oft the setiieeis power > consideration ns 
_ stated in connection with steps 3 and 4. Thus the equilibrium requiremen — 
are automatically simulated by current- summation requirements although cur- 


rents are not “necessarily taken into explicit account in devising the 


An ‘ites will now be made to improve the analog by altering the ca capaci- 
= and inductances so that the energy of the electrical analog will more 


provements can by the method, their in- 
— by this method seems less obvious and will not be undertaken here. = 

First t considering the kinetic energy, we shall change Eq. 16 so as to take 
account of the fact that mass is in actuality distributed along 


FIG, 4, —ILLUSTRATIVE EXAMPL 


‘the beam. Ini om sowe shall = the local deforme ati ions caused oe shear 


aa 


ax 2 


2 6 n+ 2 Yn +4 


Next. computing the | strain energy, we assume ‘that the beam is 


_ cept at the end points; this neglects the distributed — force and moment. 


3) — 
— 
— 
— 
— 
— — 
— 
— 
ia 
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= 


the entire beam — 

which is strikingly similar to 21, and is presumably a better approxima- 

tion than Eq. 17. The term Eg may also be — in the 

‘mations, but this is more 


into account that half of tl the two and the 
inductances Ax/(E I) refer to the nt beam increment itself, it is seen with the 
aid of Eqs. 2 and 5 that Eqs. 21 and 24 can be satisfied in the analog of Fig. 3 
_ by changing the capacitance inthe lower portion and the inductance in the upper 
portion to the following respective values: 


only positive can be used, the existence 
rotary mass moment of inertia (the first term in Eq. 25a) and shear compli- 
ance (the first term in Eq. 25b) is actually advantageous in that it it serves to — 


improve the accuracy of practical analogs. 


beam hinged at both ends, and having the constant cross. 


section shown in ar area of each flange is and may be assum ned 


— 
— 

— i 

— 
— 


BEAMS 


concentrated | ata distance r from the neutral axis. The » webs will be consid 


ered weightless for ' simplicity, and are assumed to resist all the shear but | 
Rone of the moment, the value K A 
The problem will be to determine the behavior of the beam after a central 
_ transverse force 2 P, which has” been built up very slowly, | is suddenly re- 
"4 moved. Because of symmetry we need consider only the right half of the beam, — 
_ the left end of whichis subject to zero force (after the initial force is removed) _ 
and is held to zero slope by an unknown, time-varying moment. * ~The right and 
experiences zero deflection and zero moment. ‘Fig. 4(a) shows how this half 
is divided into three equal incremental lengths, and Fig. 5 shows how the con- 
ventional analog of Fig. 3 can be arranged with suitable end- conditions to 
simulate these The following numerical values” have been 


transformer turns are all unity as indicated, with relative wind- 
ns. As can be judged by inspection of © 


FIG, 6. —IMPROVED ANALOG oF HA \LF 


of Eos. 25 to Fig. 3. The upper inductances and lower capacitances have been 
4 mitted because their values are negative, being theoretically equal to - 1/24 
_ The accuracy of the analogs of Figs. 5 and 6 has been determined by an ana- 
lytical means, on the assumption of negligible imperfections of the electrical 


= initial currents in the ae portions of the ; analogs by s separate aux h; 

-iliary means, which is however not difficult. On comparing the initial analog Pos 

_ behavior with the theoretically exact behavior of the continuous noo, is 


a 
q 
| 
— 
4 
-BEAM 
Fig. 5, the effects of direct sheal Compliance and rotary inertia are small but 
components. Considering first the initial conditions, it is theoretically neces-— wilde 
Si ae - sary only to assume that the current generator in Figs. 5 and 6 supplies a 3 oT 
a ‘steady direct current as shown. In practice, because no transformer is even if 
— 
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found ‘that all the ‘initial flexural slopes” half-stations in both analogs 
have the same error, namely 2. 8% of the maximum initial slope o of 4.5 (in es 
magnitude) at the right end of the beam. The initial deflection errors Ay, ex- 

pressed as percentages of the maximum initial deflection Yo, are shown in 
Fig. 7 by curves 3C and 3I for the conventional and improved analogs, respec- 

a tively; the increased accuracy of the improved analog is plainly evident. In- a 
= spection of the relative size of the inertia-force capacitors (in the upper por- 

x tion of the circuit) and inertia-moment capacitors (in the lower portion) in a 
Fig. 5 shows that deflections, in general, are far more important than | Slopes: ; 

_ in determining the kinetic energy, although slopes n may be important in deter- 
mining the elastic emergy. 0 
we In order to get some idea of the effect of an increased number of finite- 

_ difference increments, the initial errors have also been worked out for the 
four-increment analogs ‘corresponding to Figs. 5 and 6. In the improved ana- | 
a it is found that the lower capacitances are still negative, but the upper in- a 


> 


are shown in Fig. 7 by curves s4C and 41 for the conventional and im- 
proved analogs, respectively. The” half-station initial slope are ‘all 

The deflection errors in the dynamic case have been obtained only for the 
_three-increment analogs of Figs. 5 and 6, by comparing the response of these 
_ analogs with that of the actual continuous beam. For a time duration not ex- ¥ 
ceeding one full cycle of the fundamental mode of vibration, the latter response 
can be expressed to within an accuracy of Lae of the maximum initial de- 


a 9. cos 0.2687 t - 0.143 cin 232 cos 2.1188 t 


sin 


| 
‘ 
i 

— 

s where x is measured from the right end of the half-beam. The errors Ayo, 
a _ Ayy, and Ayp at stations 0, 1, and 2 of the beam, expressed as percentages of —— 


BEAMS 


tic. Further, the improved analog has fewer electrical components than the 
conventional analog, for the same number of increments. _ 


CONCLLUSIO 


_ It seems reasonable to expect that these same advantages of the > improved 
_ analog may apply to many other types of beams than the one treated previously, 


4 


RELATIVE DE FLECTION ERROF ERRORS AS AS A FUNCTION: OF TIME S. 


For a given allowable overall error, the number of 
‘The average number of electrical components per beam increment may 
ao Because the total number of electrical « components is reduced asa re-— 
“sult. of the above two considerations, the sensitivity of the analog to unavoid- — 


able imperfections of the electrical components may be reduced, thus rther 


| the maximum deflection in the initial condition, are shown Fig in 
siting fram tha nea af tho imnrarad analagw ic canan 
— 
— 
— 


ing their cost since this is usually a sensitive the allowable imper-— 


ery way a expressing these advantages is that, for a given analog in in- 
Stallation, the complexity of which may be handled with a given < 


is given to Melvin | Zaid, of Teknik, Inc., Garden | City, N.Y, : for 
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edie PLASTIC ANALYSIS OF TRANSVERSELY LOADED PLATES 


_ The flexural behavior of thin plates in the post elastic range is discussed. - i 
elastic, purely-plastic relationship between moment and curvature is as- __ 
‘ sumed, and general solutions developed from small deflection theory, for axi- 


q symmetrical yielding over an area, and axi-symmetrical and unsymmetrical _ 


Applications include the inelastic analysis | of metal plates, ‘pavement slabs, 


The study of inelastic an in plates might be considered to hewn i 
_ oped along two paths, the classical approach used for metals, and the rather 
more phenomenalogical basis of the yield line theory for reinforced concrete. - 
_ The classical study of plasticity appears to have originated in 1864, when 
Tresca (1)2 published a preliminary ‘account of punching and extrusion tests, a 
= and deduced that yielding occurred when the maximum shearing stress attained 7 
_a critical value. Among various other criteria suggested, was that of von — 
_ (2), derived ‘mathematically, and shown by Hencky (3) to imply that the by 
elastic shear strain energy is the critical factor. Considerable progress has 
' ~ since been made, particularly during the 1950's, in the case of plates, but the © 


"criteria of Tresca and von Mises remainasabasisfor the theory. © 


4 


_ Note.—Discussion open until March 1, 1961. To extend the closing date one month, = 
- a written request must be filed with the Executive Secretary, ASCE, This paper is part — 
- of the copyrighted Journal of the Engineering Mechanics Division, ’ Proceedings of the 
_ American Society of Civil Engineers, Vol. 86, No. EM 5, October, 1960, = =~ 


1 Senior Structural Engr., Commonwealth Dept. of Works, Melbourne, Australia, 
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materials is based on the work of Bach (4), who determined the yield moment > 

for a square slab or plate, and of Ingerslev (5), who extended the theory to 

: rectangular plates. The theory has been subsequently developed by Johansen 
(6), (7), (8), to the point where it may be simply applied to the determination 4 


a of ultimate | flexural capacity for a wide range of loading and ‘Support conditions. 


i A description (in English) of the yield line theory for reinforced eabieele ; 


plates is presented in a paper by Hognestad (9), and ercgnnee in the classical 


theory has been summarized by Hodge (10). 
nl The emphasis, in general, has been centered on the niet of methods © - 


for limit design, where yielding ata boundary or ‘throughout | the plate normally 
the problem to simple statics. _ A knowledge of behavior in the 


> 


a 


g 


‘FIG, ELASTIC (a), AND SIM-_ FIG, 2.—MOMENT—CURVATURE 


Y _ Notation. —The letter symbols adopted for use inthis paper a are defined where = 
they first appear, in the illustrations or in the text, and are arranged alpha- = 
- betically, for convenience of reference, in Appendix I IL. 


mal deft 
pe to deflections in 


and the case of linear strain hardening has been treated separately elsewhere 
a The distribution of pure bending stress in the elastic range is assumed to be 
_as shown in Fig. 1(a), and in the post elastic range, in Fig. 1(b), and conse-— 
quently, the relationship between moment and curvature may be awe 


t2 


EMS 

— 
— 
- 
— 
— — 
"particularly in thin plates, may be the most severe criterion of failure. 
l—> 

_General.—The present paper is limited in scc = 
ideally elastic-plastic plate with no strain harde1 — 
7 | sidered initially to be uniformly thin and to be composed of a homogeneous, ~ 7 _ 


~ in which oy is the yield stress in simple tension or compression, and t is the ¢ 
_ thickness of the plate. The magnitude of the yieldmoment, when other stresses 
= present, may be obtained from the yield criteria of Tresca or von Mises, : 
which have been found to des the simplest and most satisfactory for © iectrop 
_— Tresea’ s yield stress s criterion may be written in the form 


in which and are stresses. The constants ky and ko 
, and are ky and 
for yield moments from these stresses, the e orientation 
of the axes of reference with respect to the directions of principal stress must 
& be considered. It follows also, that the distribution of bending and shear stresses 
at the critical section must be determined or assumed. Integration over the 
depth of the plate then gives the relationships required (see Appendix I). a 
the particular case of axi-symmetrical bending, using polar coordinates, 
_ if radialshear is assumed, for simplicity, to be carriedonly by a central core, 
= 
the Tresca criterion for plastic radial moment 


where bending 1 is in the second 2a moments are of opposite sign, 
For convention used see Fig. 3. Yield moments in the third and fourth quad- -_ 
‘rants are obtained by substituting - “hed for My in Eas. 4 through 10. The sal ; 


determined from | the appropriate yield criterion. 
o#, for comparison, the shear stress were small and were arbitrarily assumed 


to be distributed uniformly ery the depth of the plate, the Tresca criterion 


a 
— 
=f 
— 
—— 
— 
Tun iaurant, that is lor moments or the same sign, OF 
— 


Introducing the further that the stress for are rec 
tangular, the von Mises criterion becomes 


- M, M 
In particular cases of axi- symmetrical bending, a a distribution of shear and 
tS, bending stresses may be determined from the equations for equilibrium of an 
element and the Tresca yield criterion, by extending the method of aun 
described by Drucker (11) for a rectangular beam. For the limiting | cases 

_ considered in Appendix I, the resultant radial yield moment is given by — 


nore bonding is in the first quadrant, 


— | might, are to bad maghocted. 


vided these conditions are reasonably satisfied, as discussed in " Appendix L ; 
_ Eqs. 9 apply to yielding at an internal, or external, circular boundary, and in 
the limit where the radius of the critical section approaches to one- 
The stress ‘@stributions. determined in each case lay the two 
_ tremes previously assumed and produce a greater moment of resistance than 
Generally, effects will be insignificant. The slight moment redection 
due to shear tends to be offset by a higher yield stress in bending than is ob- 
tained in simple tension, , apparently due to non-uniform strain anda smaller 
, i volume of yield involved. Furthermore, at points of maximum shear the mom- 
ent curve is generally steepest and the above factors are The 
moment criteria of Tresca and von Mises reduce ‘respectively to 


way - Mr Mg + 


4 
i 
— 
4 
(2) 
= & for bending inthe second quadrant. The assumption is made in each case, that 
loading in the regio 
and that the powers | 
— 
— 
ax 
— 
4 hee 
oor 
— These relationships ar 4 
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_ rectangular coordinates (u, v) by substituting u and v, in which u and v are 
directions of penaciees moment, in the case of f Eqs. 4 through 1 11, for r and ¢ in 
_ ss ielding.— —Where the requirements of statics or the geometry of the plate 
prescribe that moments should vary in the directions of yield, yielding will tend 
to be confined to lines, and the plate will remain elastic except at these lines. q 


- Where no such requirement exists, yielding may extend over an area of the plate. 
___ _ In the case of first quadrant bending, yielding overan area may occur either 
= n one direction of principal moment only, in which case the plate behaves elas- 
_ tically in the other, or yielding may occur in each direction, causing the plate 
~— act essentially as a membrane under additional loads. When bending is in 

the second quadrant, yielding occurs simultaneously in each direction of prin- 
moment, but ina predictable manner. PES 


_A simple example in which yielding is confined to lines is a uniformly londed = 


te way plate, clamped at each end as shown in Fig. 4. In the elastic range of 
7 stresses, maximum moment occurs at the supports, with maximum stress il 
_ the extreme fibers; it is here that yielding commences. Under increasing load, 


(Fig. 2) is reached. Further load produces r rotation a at the ends without in- 
crease in restraining moment, so that, in effect, a “plastic hinge” is formed © 
at each support. Under additional loads, the plate behaves as if it were pinned : 7 


at each end, the moment at the center increasing until it too reaches Mp, anda 
third hinge is formed. An arbitrarily small load increase produces unlimited — a 
‘rotation about these hinges causing the plate to fold and collapse as a mech- 
De The zones of yieldare shown in Fig. 4(c). ‘Since the width of yield is greater 


- at the midspan than at the supports, because the ‘different slopes of the eyed : 


“the (actual) maximum elastic moment Me is exceeded where 


“plate. For the purpose e of analysis, each band is considered to be 1 reduced, in 
width, to a line, and is termed ayield line. This is consistent with the approxi- 
introduced in ‘Fig. 2, in which the plate is considered to remain elastic | 


> 


The lines at the the plate (Fig. 4(e)) are produced | by negative mom-— 
ents and are termed negative yield 1 lines. The line at the center is correspond- 
ingly termed a positive yield” line. Further distinction might be made on aa he 


basis of whether or not shear forces are present. 
‘The three lines together form a yield line pattern. Yield line patterns might 
"themselves be classified by the degree to which they contribute to failure by | 


collapse of the plate, and on oe basis, a a pattern might be divided into one of 
_ three groups 


— 
= 
q 
— 
— 
— 
} 
ere 


‘The first group consists that form a = 
ism, allowing collapse by folding of the plate as in the previous example. ketal = 
‘The ec: nd group 0 of patterns requires some contributing cause such asten- 
sion, compression, or shear failure to produce collapse. A uniformly loaded 
square plate, simply supported on all sides and producing the yield line pattern : 
- shown in Fig. 5, is typical of this group. Here a purely flexural failure is im- _ 
possible and tensile yielding across the diagonals or buckling along the edges 7 
is necessary to produce collapse. _ 
_ The third group consists of yield patterns that are active only during a cer- 
: tain stage of deformation in the plate and do not contribute directly to collapse. 
The line of distinction between groups two and three! is not always particularly © 


5, -YIELD LINES, UNIFORMLY LOADED, SIMPLY 
SUPPORTED, SQUARE PLATE (NEGLECTING 


Patterns of the first and second groups are Fo capable of analysis by 
whereas group generally a1 are ‘not; 
_ Thus, for bending in the inelastic range , yielding ee be confined | to a 
or may extend over areas of the plate. If the moment curvature relationship 12 
is assumed to be approximated by ‘Fig. 2, the lines may be considered of — 


width. For the same reason, the yield areas may be considered to have def- — 


— 
— 
— 
— 
— 
— 
— 
q 
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. Method of Analysis. —For the purpose of analysis, the plate may be divided — 
into elastic and or areas, and each zone may be considered 
analysis of a on type elastic foundation will beconsidered 
7. initially, as a general case. The solution for an edge supported plate may then 
a be obtained from this as the limiting case where the density or reaction mod- a 


i ulus of the a liquid approaches zero. In certain cases, such as ol 


but finite value in order to reduce the interaction of loads and reactions, and — 
Considering, firstly, the elastic portions s of the p plate, the differential e 


for deflection, w, is 
4 


= 
which q is the intensity of lateral y is the density or reaction 


of the foundation, D is the flexural of the plate, and is the 


‘a point load on the plate 


EB A, J, (x vi) +B, 
es where J, ( ) and Hy ( ) are Bessel functions of the anit and third kinds re- — 
"spectively, of n, and i= =Vv-1, The constants Ap, and Bn are complex 
are evaluated from the boundary conditions. functions J, ( ) (1), Jj ( 
)( ), and H,(1) ( ), or their mathematical equivalents, are tabulated 
pe - elsewhere (16 6), (17), (18), and functions of higher order may be obtained from 
_ Where there is more than one internal boundary in the plate, it is conven- 
ient to consider the two functions of x in the solution separately and super- ; 
impose the results. Thus, the plate may be considered to be initially infinite 
by and the conditions at the positions of the boundaries satisfied by superposition. Ra, a 
_ The external problem of deflection produced in the surrounding plate by a we 
deformation at or * within an 1 internal _ boundary, whether due to loading, temper- 


owe} 


) 


cosine series only, by choosing a suitable origin for ¢. 


2 “a _ Likewise, the internal problem of deflection within an external boundary 
produced by forces or 7 yielding at or outside the << + is given by the series 


jcos no my 


— — 
ae 
= 
a 
am 
z= 
il 
and ¢ are in 1884. +) Intro st the solution 
is the radius o 
a : 
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> 
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a 
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‘Thus, | Eq. 16 is used to describe the conditions at « each sesieiiiaitiiaaiaiaille in 
_ infinite plate and Eq. 17 to satisfy the conditions at the external boundary. = 
_ The constants in each equation ‘ are evaluated with = loads and boundary e ef- : 


Where the boundaries remain elastic, the fora clamped boundary 
4 are that w and = pr 2re zero, ino which n is the normal to the edge = Similarly, 


for. a a simply edge, w w= - 0 and M, ata free ¢ edge, M M, = 02 and 
0. (Vy = Qy - Myy/ay, ‘the effective shear at the edge.) In the case of 
an integral elastic support the conditions of continuity must be fulfilled, and 7 
in the case of a non-integral elastic support the displacements must be com- 
_ Where a yield line occurs at a clamped boundary, the boundary conditions 
are w = 0 as before, and My = Mp, where Mp is the appropriate plastic moment. ; 
The subscript n, in this case, is the direction of yield, which does not neces- — 
sarily remain normal tothe (original) line as loadincreases. Where a straight a 


_ yield line develops within the plate, , and not at a boundary (of a support, load, 


My = Mp ‘and Q, = 0 satisfied, except atan intersection of more than 
q ye, two lines, or at points where the line meets a boundary in the plate. Deflec- 
tions, in each case, are assumed to be continuous at the line. ie all 4 
- An incomplete line may be considered as an internal boundary and be de- 
a scribed by the series of Eq. 16. Alternatively, the moment condition M, = Mp 
- might be considered tobe satisfied by a rowof coupled point loads, or concen 
trated moments, acting along each edge of the line, in the normal direction n. 


The deflection produced by wa cad pair of these concentrated moments is 5 


where C is a real constant, wi @ is measured from the direction — 4 
tion along the line of yield gives the solution required. =  — cm 
yield line at a boundary is incomplete, the appropriate elastic and yield 
conditions are appliedat the respective sections involved. The extent of yield, 
in the general case, is necessarily determined bytrial, = | ae 
_ Hence, where yielding is confined to a line it may be treated as a boundary — e 
condition. . Expressions for moments and shears in terms of deflections, in the 7 
elastic zones of the plate, are available elsewhere (15), (20). a ee : a 
Where bending is axi-symmetrical, and yielding extends over an area, ex-_ “ 
 . pressions for deflection may be obtained by modifying the procedures used for — 


a purely elastic zone. The maximum shearing stress criterion will be adopted — 


Jt in this case, and for yielding in the first three different conditions are 


Mg = My > Oz). 192) 
and plastic strain is to the zis the direction normal 


to the plate. The direction of the symbol > in Eq. 19a refers to moments in " 
the first quadrant of Fig. 3, _ stresses in the lower half of the plate. For 


— 
5 
ay 
3 
— 
| 
— 
q 
‘ 


= 


_ behaves essentially as a membrane under additional loads. Hence, the bending 
_ the first quadrant, only the case of tangential yielding (Eq. 19a) need be con- 


sidered further. For bending in the me petenat, the yield criterion is * 
Mg = My +My (a ( ) 


plasticity | is restricted to plane. 


first quadrant and tangential yielding only, the corresponding 
equation, from Appendix] is 


_ yielding i in the second quadrant, the deflection, w, is s conveniently divided 


(q - yw 


Primes in the e preceding equations denote differentiation with respect to 


in which represents Poisson’s ratio. jo. Introducing the variable x = r/s where. 


1s. 21, 22, and tively 


- _ Here plasticity is confined to the rz plane and the only case of practical im-— > 
[am 6@€6—FStiérrttzncce’ is thatt in which the yield area is reduced in width toa line. For the 7 _— of 
\ of h each case are derived in Ap- = 
pendix I, and it is of interest tocompare them with the corresponding equation 
fo an elastic zone. For symmetrical bending in the elastic range, the differ- _ 
— 
in which We is produced by the strain and Wy Dy 
— 
— 
— 
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w we 


where x, in this case, is introduced only for dimensional consistency with Eqs. 


the w wis 


ma are independent s solutions and are real constants. 


2. a(n- (n- 2) ).(n- 4) (n- -5)* 37.2. 


5,9,13.. -2) 
an 


Values of Rj and their diravatives| R" R™)) in the range 0< x <5 
are given in Table 1 and elsewhere | eee ha values of x inthe range 0(0. 1)10). 


5 
— in this case, and from 
q 
— q 
— 
— — 
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‘We we= Cy » 


4.3. 


4) (n-5) -6): 
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For the particular r case of an edge ‘supported the solution of Eq. 


log x loge x 


for and first ‘quadrant bending, the e deflection from 


w= +e + x? + C4 Loge x 


| 


Eqs. 37 to 39a might alternatively have been obtained from the condition that, © : 
in the limit we 0, x(= r/s) steo approaches zero, and that only the first term 


= 


For yielding in the second quadrant, Eq. 28b gives 
| 2 S2 + C3 Sg + C4 Sq + C5 S5 - Cg = Cj 
which 5; through 54 ar ionsto thehomogeneous equa- 
_tion (in we), and (C5 S integral. i 
4 
= a ) ‘ q 
“Ne 
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in each series on for the e plate. on a liquid support is ciate 


Rg and Ry loge x 
so that the ‘solution for w in Eq. 33 approaches the solution for \ w in Eq. 38. — 
Cases of nearly symmetrical yielding may be approximated by the symmet- “a 
rical case, but less symmetrical yielding is not so easily analyzed. _ Fortunately, — 
i where bending | is unsymmetrical, yielding tends to a greater extent to be con- 
fined to (single) lines, so that a wide range of problems may still be readily — 
a For ‘symmetrical bending in the elastic range, moments and shears may be 
expressed as derivatives of deflection. As example, 
dr2 


= 


dr 
In the case of first and tangential ‘yield, t the corresponding 


relationships a are, from Appendix L 


and the value of Mg from Eq, 192 ho holds. 
} 


(45) 


and the value of M, given in Eq. 20 holds. 
? Hence, for the case of a eee on an elastic foundation and symmetrical ben i- 


— 
De 
‘eo 
which 


Cj + C5 loge x| j= 1, 


“4 - 


-(1-y) 4, +C4 


| 


C; + C4 x loge 


— 
For the case (58) 
clastic range, setting 
ii 
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and in in Ge second 


4(2 x + 


Plate on Elastic Foundation ‘with Concentrically Loaded Circular Insert 
_ Consider the case of a . large « or infinite plate on a liquid type elastic founda- 
tion, with a rigid circular insert, concentrically loaded (Fig. 6). er 
_ In the elastic range of stresses, the deflection is given by Eq. 48; the con-— 
stants may be evaluated from the conditions of zero deflection and slope at in- 
-finity, zero slope at the edge of the insert, and from the shear at this point. 
Hence, from Eq. 49, the maximum moment occurs at the edge of the insert, in 
the radial direction, and it is here that yielding commences. Yielding about 
“2 the edge of the insert allows rotation of the plate, and a yield line is formed. 
Further rotation, with load increase, allows tangential moment to build up until “7 


a reaches the yieldmoment. 


Under increased load, yielding occurs in the tangential direction at or near 
4 the edge of the insert and spreads concentrically outwards. At any particular q 
load, the edge of yield may be determined and the constants evaluated, as fol- a 
lows. Since bending in the region of the insert is in the first quadrant, deflec- a4 
— tion within the yielded zone is given by Eq. 52, and radial moments and shears 


by Eqs. 53 The conditions, at the edge of the insert, 


My =M from Eq. 10 and Qr = a oe » in which rp is is the radius o of the in- 


_sert, and, at the outer edge of the yielded zone, Mg = My, and continuity in the ‘A 


plate. In the ‘surrounding elastic portion of the plate, the constants are evalu- 
_ated from these conditions at the edge of the yielded zone, and the condition that 
_ slope and deflection approach zero for large values of r. For the particular 
case of rp = 0.4 is and p= ( 0, _ the ‘moments at at _ various» | stages of of yield are as_ 


Where the plate is thick and yielding is considered to be confined to a line, 


effect of shear deformation in the may be included using 


Re| 


arcos 


and the complex constant B may be evaluated in this case by introducing the 
plastic moment = in which 7 by E 9, and from the 


— 
— 
4 
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B Mp=!.5P, 

P=I10.78Mp=215R / {| ! \ Mp=2.8P 8 
FIG, 6,—MOMENTS AROUND A LOADED INSERT IN A PLATE ON 
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— of the shear at the insert edge. General expressions for elastic 


‘Simply Supported Circular Plate with Concentric Ring Load.—The case of a 1 


circular plate, simply supported at the edge and with a ring load concentric 
about the center (Fig. 7), will now be considered. This | case is treated ed by 

pes the elastic range the moments are greatest within the loaded area where 
a state of pure bending exists; that is M, = Mg = Mo (constant). Under in-— 
creasing load, Mo i increases to Mp so that the moment everywhere within the 
ee circle is at the yield moment, and pure, two dimensional yielding m may 
occur. Rotation of the plate at the edge of the loaded area, under increasing 
load, prodhoes an increase in tangential moments in the surrounding plate, and 
=. tangential yielding extends gradually outwards to the edge of the plate. ie a 7 
: _ The bending is ee in the first quadrant, so that, for yielding within the 


loaded area, M, = M iy = = My, and for tangential yielding in the surrounding =, 


“FIG, _7,—SIMPLY SUPPORTED CIRCULAR PLATE were RING OAD P. 


M . = My. Deflection and radial moment are given by Eqs. 62 and 63. In the i 
case where the tangential moment just reaches Mya at the exterior edge, the 


constants cy through C4 are » determined from the conditions M; = My and Qr 4 


at the edge of the load, and w = Oand Mg = DD, - — at wae 


=0 this 


— | 4 
— 

— 

— 


‘in which. a i is the radius of the plate and b s), the radius of the 
loaded circle. The magnitude of the load to produce this condition —< 


7 = Uniformly Loaded Circular Plate. nil iain loaded plate 1 may be treated 

_ similarly. _ This case is already treated by Tekinalp (14), but is included to 


tion due to the uniform load is given by integration (of Eq. 105 in Appendix -_ 


w = At for first quadrant wendy and this when added to the homoge- 


_ of a plate on an elastic foundation, . The load, q, may then be considered to ex- - 
tend to ed no only 


~x? 


+ Cg x? 


‘The first two terms of series must be here since 


_ is another constant term in the equation and = 72 not be in com- 


4 
1 


Bron with \ their sum. Hence, for a uniformly loaded plate, phil supported, 


teal th an internal panel of the plate, ,under a steadily i increasing, uniform : 


plastic moment, M, is constant in each quadrant. 
‘The plate pote, Sy readily analyzed in the elastic range by introducing a - 
small but finite liquid reaction and using superposition techniques, as described 
elsewhere (19), (20). In this way, the plate may be considered as initially 
; - infinite in extent andsupported only by the liquid. Each load and reaction | may 
_ then be separately applied. The uniform load alone produces no moments or — 
shears in the floating plate, only a uniform deflection. Superimposing the ef- ; 
fects of each support reaction effectively eliminates the liquid reaction, leaving 
only an arbitrarily ‘small reaction due to the actual, net deflections 


n€OUS SOLUTIO q. al deliection ii the plate. Ihe same &g 
may be obtained by considering the edge supported plate as the limiting case 
w= =, al 7 
civing 
4 
[3(3 + 2v) - 4(1 + - 6(1 - v) rg2r? + (73) = 
_ Continuous Plate, Intermittent Supports .—Slightly more complex, but still 
readily handled during certain stages of yield, is the caseofa continuous plate 
, 
— 


is, in n effect, the actual moment in the plate. Shear and deflection 


‘Ai: _ Let the supports be arranged in a square grid of spacing, L, and let each 
“have an effective radius, rp, of 0. 05, ‘If the supports are assumed to | com- 


Bending in the region of a is ‘essentially about ts 
center. Maximum elastic moment occurs at the edge of o support, in ee ra- oa 


-Loadonpanel=P 


IG, 8.—CONTINUOUS PLATE, UNIFORMLY- LOADED SQUARE, INTERNAL PANEL, 
"MOMENTS AT CENTER-LINE OF SUPPORTS, (V = 0) 


_ formed around the edge of the support (Fig. 9(a)). Rotation about this line, 
ae a under increasing load, allows the tangential moment at the edge to gradually — 
_ inerease from zero (for p = 0, Fig. 8(b)) to Mp (Fig. 8(c)). Further load 
produces tangential yielding which commences at or near the edge of the sup- 
_ port andspreads concentrically, or nearly concentrically, o outwards inthe | sur- 
ee Bending is nearly symmetrical only in the: region of ies supports. As yiel 
7 ing spreads outside this area it tends toextend farthest along the lines of max 
mum curvature, which, on the assumption that the shear is uniformly distrib- — 
uted around the support, are along the panel edges (Fig. 9(c)). _Hence, the 
zone of yielding tends to narrow to a line, eventually reaching midspan. a>? 
Meanwhile, positive yielding commences s midway t between the supports, whe 


A} positive moment is a maximum, and A epreads with — load, —— form a 


— 
—- 
— 
— 


iy line along the panel center. Thus, a a yield suennnate of th the , third kind is completed ip? 
oe The plate will be assumed at this stage to develop a preference (to an ar- 
bitrarily small to yield in one direction of span. Since, in 


‘straight and thus capable of unlimited rotation, the negative line curves — oa 

face of the supports where radial shear resistance prevents collapse. 

= j. Limitedyielding about these lines in the preferred direction, causes cai 
bending in the region of the support to be less symmetrical. This produces a 
significant variation in the distribution of shear force at the edge, and a con-_ 

_ sequent change in the directions of yield. ‘Thus, a straight yield line tends to 
form along the edge of the supports (Fig. 9(d) ), and completion of this line re- 
sults in a pattern of the first group. Further load causes unlimited ha 
along this line and along the positive line line at the re hints the plate x 


Moments in the elastic range are “determined from Eqs. 5. 49 and 50, tor the 


conditions, at the edge of of the Support, that radial slope is zero and 
in which P is the panel oad. The load, at which radial viel commences, 


are found f the Qp=3 


My = My - = at edge. 


During the range of symmetrical, or nearly symmetrical, tangential yield- 


_ ing around the supports, moments are determined from the conditions My = Mp 

= “Tam at the edge of the s supports, and the conditions Mg = Mp, , and 
4 of at the outer perimeter r of yield. Bending within the yielded zone 
is described by Eqs. 52 through 54, or by Eas. 62 through 64, provided me ef- 


_ of the distributed load is is separately included in 1 the latter c case. 


at 1 which a pattern ‘of the third group is ‘completed, and ‘the load P4, at which a 

_ pattern of the first group is formed, and consequent collapse occurs, are each — 
gee statics. Hence, moments in the elastic range, up to load Py, are as 
shown in Fig. 8(b). Moments at the load Po, at which tangential yield i a 
‘mences, are given in Fig. 8(c), and moments at a load greater than <a 
-sibly near the limit of reasonably symmetrical are 8(d). 


q 


7.46 M 0. 377 


— — 
— 
at this stage of yielding, its eifect will tend to become significant. A purely 
— — 
at which tangential vield commences, moments in the loading range P1<P<PQ 
— 
— — 
— 
— 
g 
— 


== Radial yield 


hinges 


7 
NEGATIVE. 
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— 
\ 


1960 
0 ter a square panel, and, from statics, 


ix 


which is the ultimate flexural capacity of the panel. 
_ Whereas load Pg is based on the assumption of a uniform shear distribution _ cs 
around the edge of the support, load P4 requires the shear to be (largely) con- 7 
centrated at the point on the edge to which the straight yield line is tangential 
(in Fig. 9(d)). If, in an actual plate, a preferred direction of span were to 
develop significantly at a load less than Pg, shear variations would occur at 


Loads P, through P4 are computed for the condition that Mp is ¢ constant ine 
- each quadrant. For the Tresca yield criterion, Py and Pg would remain i a 
- same, but Py and P4 would each be reduced, since part of the yielding, i 
case, is produced by bending in the second quadrant. ~ 
baal The same procedure of analysis may be used for thick plates in the elastic 
_ range of stress (20). It might be extended also to the inelastic range, for the - 
condition that yielding is confined to an axi- symmetrical line around the ll 7 
ofthe support, asinthefirstexample, 
Plate on Elastic Foundation, Effect of Local Thickening.— —Consider the — 


"surrounding and concentric with the loaded rigid insert. 
a _— in the elastic range are determined from the equations | presented. a 


ion: ro slope, at the edge 


- 


- Due to the different plastic moments of the plate and the flexible insert, ra- 
_ dial yielding may first commence at the internal boundary of either, depending © 
on the relative proportions involved. Under increasing load, yield lines will 
eventually form at both these boundaries. Rotation about these lines allows — 

i tangential moments to increase, and at greater loads, zones of tangential yield 7 

are also produced, commencing: at the radial yield lines and spreading outwards — : 
i. position at which yielding first commences is determined from the elas- 

= Mpi is substituted 


= 
in which k = = 
— 
— 
| 
— 
q 
— 
— 
yielding commences at the other internal boundary may then be similarly de- 
termined. Under additional loads, the plastic‘moment conditionis used ateach 
edge, and moments are again determined from the previous equations. Sim-_ a 
ilarly, the loads at which tangential yielding commences may be calculated and 
plate analyzed during tangential yield, a 


ELASTIC- PLASTIC ANALYSIS 


preciated. 


_ The above procedure has application also to the corresponding problem in 
a continuous plate with a ‘supports, as will be readily appr 
i _ Thus, yielding may be _ confined to lines or may occur over an area of the 
plate. Yielding over an area may be in one direction or intwo. Yield lines, © 
it is shown, may be straight or curved, and a curve might be considered the 
general form. Lines producing collapse, however, must be seamieenineal 
_ - Yielding patterns have been nominally classified into three groups. The 
_ formation of any of these groups may be nominated asa design criterion using 


factored loads, although since they contribute in different degrees to os 
collapse, it would be reasonable to use different load factors for each. Sih | 
_ - Yield line patterns of the first group produce the ultimate flexural capacity 
and collapse of the plate as a mechanism. _ Yield patterns of the second group 
generally also represent the ultimate flexural capacity of the plate, but collapse 
in this caseis resisted bythe additional capacity produced by in-plane or shear 


Yield patterns of the third group do not normally represent ultimate flexural 
capacity and do not properly fit into the yield line theory as such. The examples | 7 
> considered herein were generally of this group, and the case of symmetrical 
__ bending, in particular, may be readily handled by the solutions presented. ae 
In the case of symmetrical bending about an insert, the condition that yield- 
4% ing is confined to aline might be used as a design criterion, with ‘suitably fac- 
tored loads. In this case the plate may be analyzed by conventional elastic 
"methods except that at the insert edge the plastic boundary condition M;r=M 
=M is substituted for the elastic condition n of : zero slope. Moments and 
4 flections in the ‘surrounding plate are the same as in a continuous, purely elas- 
- tic plate with a rigid insert merely bearing on it, so that utilization of this de- _ 
gree of plasticity is by no means unjustified. If symmetrical tangential yield a 
= is permitted, the solutions in this region are as simple as in the elastic case. z 
a _ Although in the preceding theory a homogeneous, isotropic material is as- om 
_ sumed, use is made of the assumption only in the calculation of the moment 
7 ordinate Mp; many of the solutions presented may be applied to other cases = 


gy te plates, or slabs. In 


including 


struction in particular, the relationship between moment and curvature is ap- 
proximated fairly well by the idealized curve (4), (5), (6), (7), (8), (9), of 
‘Fig. 2. The moment “ordinate, represents, in effect, the point at which 
yielding commences in the reinforcement, which may be estimated with rea- - 
_ sonable aang and may be considered to be independent of the other stresses r 
int the plate. . Hence, , the equations» presented for first quadrant cael en 


‘metal plate ‘may be used in each quadrant of the »slab. 
_ The fact that deflections may be determined in the inelastic : range isa a 


“of considerable in which deflection the most 


small areas. 


— 
— 
q 
q 
Fig 
tes 
= 
rt 
the ranges OF reintorcoment rating nracticanie tor clahe and tor tliat 
— 
— 
— — — 
advantage of utilizing the plastic properties of materials in plate design a 
tial, particularly when loads d — 
an particularly when loads or reactionsare concentratedover 


i “at the case of metal plates, yielding through the depth of the plate i increases : 
the moment capacity per unit width, while yielding along a line or over an area 
-: "eliminates moment concentrations producing a more favorable distribution of 


stress in the other two dimensions. Thus, considerable material savings are 


- ents in plan again provides significant savings . Instability does not enter the 
problem as it does in the case of steel and the effect | of local — 


: lished, permissable stress methods of design is a always so great, because 
_ in many cases renee for plasticity has already been made either from ex- 
_ Inelastic analysis” in cases at least allows design | practice to be re- re- 
examined, providing a justification for many of the reductions and simplifica- 

_ tions introduced, and allowing their extension outside the limits of previous — 
experience, and to new fields and materials of | design. Connected closely to 

: ‘economy is the question of safety factors. Permissable stress methods of de- - 
te generally result in inconsistent factors of safety for different elements of — 


a structure. , Using inelastic methods, the in the post elastic — 


ability of failure may be retained, as in conventional design, savings 
in materials and costs. The same material quantities might also be used, but | 
_ proportioned to provide both a larger marginfor safety and better 


Solutions are presented for small deflections due to axi- symmetrical, and 
"unsymmetrical, yielding along lines ina plate and axi- yielding 

- obtained are applicable also to isotropically reinforced concrete plates, ‘since 

similar moment curvature relationship may be assumed for each. 

— _ The solutions for axi-symmetrical yielding are relatively simple, and 7 


suitable for use in design. They may be extended also to cases of nearly sym- G 
metrical including yielding around cylindrical st pports ina continuous 


_ gtresses at a line of yield may be determined from the yield stress criterion 
and tl the e equations of of an element. For axi- -symmetrical bending, 


— 

| 

&g — 
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} 
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— 
% — 


convention 


“The: wide of coordinates is at the middle surface of the plate. If o, in Eq. 7 


is considered identically equi ual to zero, the term r Trz is constant with respect 
_ tor, and hence the term r in Eq. 75 is similarly constant. ial 
‘Consider initially the specific case Oy = %% for values of r in the region of 


4 ana for 07 = 0, integration we respect t tor - gives. 


Og, er 


proximately so, asa ‘consequence of the condition =0. 
“= From the | Tresca criterion for yield in the first quadrant, a 


_in which Th is the radius of the critical section. . Eqs. 80 are applicable in the 


Tp loge 


the equ 
—- 
— 
a 
— — 
— 
which rg is the 
— 
— 
ing for f the plate), 
itive z (in the lower half o 
— 
— 
a 
= 


Octot 


‘the 


the ex ion for gives the 


i 7; 
over a core, and) varies sine wave on each each side of this, 


A A in’, 


although this is » eutatds the 1 range of - — that might be e expected ina aplate. ; 


The case of of axi- -symmetrical, first quadrant bending, for any distribution of 
0, with respect to r, and with a loading of uniform - intensity, q, between the - 


[1- (for z 
into and integration with tor gives 


q 


Tr 


< 1, the Tresca yield” may be expressed 


iminating o, from, Eqs. 86 and the > simplifying approxima- 


ries form and substituting Q, into 
— 

— 

— 

— 
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ri 

for Limiting caseo, = 6 
to z, and oz a cubic ehetiok of z, the resultant differential equation for Trz is 
a special case of Riccati’s equation for ' which a solution is aL but — 


in a slightly less 3 convenient form. Integrating Eq. 1, gives "ss oe 


eee ‘cosh | Va a 


truncating the series obtained, gives Eq. 9a in the range < 


1, and; again a 4 rectangular stress block for oy 


tag 


EMS 
— tion that e in the limiting case 
% = | 
— ++ (88>) 
j 
— ne yield 
)| (9%) 


| 


‘Eqs. 9 apply to at and external circular boundaries and, in 


limit r, ~ ©, to one way bending in a plate. Eqs. 82 and 83 may be derived 


_ The Eqs. 82, 83 and 9, being basedon states of stress that satisfy both equi- 
librium and a yield criterion,are lower bounds or approximations tothe lower 
pounds for the cases considered. They may also be shown to be close tothe __ 
f =: upper bounds, however, and hence, _ be consideredas approximations to the 

ae The independence of Eqs. 9 of the magnitudes of the constants a and b, that 

is Te, (radially), , and q, tends: to suggest that they may have 
wider application. The stress. Aeictetionn derived, in fact, are in the neigh- 
___ borhood of a stationary one « on the interaction curve (for st significant values 
7g eM of — in the range considered), and the effect of any dit distribution of shear stress eo 
over the depth | of the plate, from a linear (triangular) vai variation, to one ap- | 
-proaching a half sine wave, _ is approximated by Eqs. 9. Hence, where approx- 
-imations for ar and Oz ar are the effect on moment is to 


- directly, for the latter case, from the results of Drucker (11), for a beam. 7 


shear distributions previous range, the of 


uihbeuane to ‘lit di in the limit where e shear i is carried only by a central a 


or to one-half if shear is uniformly distributed over the depth of the section. 
All ‘practical c cases would ‘appear to lie between the linear d distribution and the 
ten one so that the of would be 3 and =. Over 
a substantial portion of this range, the value zis directly applicable, and pro- 
vided 


Deflection.— 

ferential equations for are e determined from “the requirements of 
- equilibrium and from the relationship in the — direction between stress — z 


Fora plate on an elastic foundation, moments, ‘shear, and the intensity of a. 


Re Eqs. 9 might be considered to have reasonably general ap- a «i 


ae and in Eq. 9b is obtained in the range oe 

‘i 

— 

— 

— 
— 

= 

| 

— 


may be by with to of 16 ‘ 
spectively. Eliminating Q,. from Eqs. 94 and 95 results in the useful relation-— 
m,)-2 


For bending in . the first quadrant and tangential yielding over an area 1 of Mo 
to . 19a an 


a and ene M, for M 


The moment and deflection is given by Eq. 43 in 


g in Eq. 96 gives 


which D, = E t3/12, and primes denote differentiation with respect to hh 
relating moments to displacements, no previous history of plastic strain: is 
assumed. The term yp represents ‘the elastic value of Poisson’s ratio, and 
_ yielding is restricted to the oz plane. 


=. 


Substituting for M,. in Eg. 97, results | 
in Eq. 22. ‘Similarly, for yielding in the radial direction ection only, Eq. 1$ 9b results sults — 


and substitutin 


(98) a 
g for My and nd Mg in Eq. 96 gives 7 ae ae 
of first yielding in Eq. 19c_ 
~ results substituting for M, and Mg in Eq. 96 leaves © a 
= 


For yielding in the second quadrant, we utilize Eq. 20 and euiiiinais for 
1 


_ Yielding occurs in the ro mene and the bending strain € may be divided into two. 


components, an elastic ‘component ,anda plastic component €, 
7 ponding deflections, produced by these strain components are denoted by We 
and Wp respectively. Ri 


. The corres- 


Radial moment is is related d only to th component, 
which results in Eq. 45 in which D2 - E 
Eq. Eq. 24a, 


& 


ull 


4 
— 
q — 
— — 
— 
— 
we 
— 
 &§ 
— 
— 


that Eq. 24b hold true. 


45 assumes €,,and €py proporticnal to z, which the Linear re- 
b lationship of momentand curvature shown in Fig. 2. Alternatively, if only total 


strain is assumed proportional to z, if ‘the Tresca criterion for yield etress 


curvature to the radial strain in the 


a For an edge supported bi y= , ae — 22 and 23a reduce respectively 


Expressions f for shear are obtained in each case for and 


we" = 


= moment per unit width of plate: 


ary = force per unit width of plate a 


= = effective she: shear per unit width 


coordinates 


= rectangular coordinates _ 


subscr ripts elastic yield in tension, we 


q 
ial moment and curvature are an ves — 
a — 
— 
= @ 
4 
— 
tress 
e, y 
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= plate thickness 


= flexural rigidity _Et 
12 (1 12(1 - 


reaction — of elastic foundation 
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| 
_ An analysis of the response of beams traversed by two-axle loads is made 
based on the assumption that the dynamic deflection is proportional to the static = f 


nt due to the weights of the beam and loads. Numerical results show 
quasi 1 resonance due to successive of loads and to wavy 


wel Problem Considered.—This paper deals with the elastic. deflections and bend- 
ing moments in simply supported beams: under the passage o of two-axle loads. 4 
The system considered is shown in Fig. 1. It consists of a simply supported | 


R beam spanned between two rigid supports and approaches, and a two-axle load = | 


moving from left to right at a constant 


flexural rigidity. In other words, the surface unevenness of the beam is as- — 


sumed to be small as | compated to its depth, so that its effect on the distribu- 


that the elastic deflection is ; due entirely to bending, the deflection of any point y 
ina cross-section of the beam is equal to the deflection of the neutral axis. a 


i< Note. —Discussion open until March 1, 1961. Separate discussions should be = 


mitted for the individual papers in this symposium, To extend the closing date one 
- month, a written request must be filed with the Executive Secretary, ASCE. This pa- 
= is part of the copyrighted Journal of the ore, Vol. 86, X Mechanics Division, al : 
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“unsprung masses” through two linearly elastic springs. ah sprung mass is 
e: ‘simply a mass that : is suspended on a spring; an unsprung mass is not. Tj. 
P Bs The usual method of dynamic analysis of elastic structures, » by expressing 

the deflection configuration of the structure as an infinite series involving its 

normal modes of vibration, is ineffective for this problem if the mass of the 
; Seo load is to be taken into account. This is due to the fact that this meth- 


od requires the normal “modes to be functions of the | space coordinate only, 


as on the Space or coordinate. 
> 


LOAD 


UNSPRUNG 


AXIS OF BEAM’ a 


hs aie to arrive at petine® method of analysis, it has been found nece 


sary by investigators of this ys of moving loads on beams, to make some 


system, a linear differential with variable coefficients. An exact 
_ solution of the equation was reported by Stokes. 4 The system of a simply sup- Se 


2 «Report of the Commissioners appointed to inquire.into the Application of Iron to Pa. 
‘eee Structures, ” Appendix by R, Willis, H. M. Stationery Office, London, 1849, 


3 “Vibration Problems in Engineering,” by S, Timoshenko and D, H, Young, 3d Ed, 


Nostrand Co., New York, 1955,p. 359, 
4 “Discussions of a Differential Equation to the Railway Bridges,” 


Fig. 1, the depth of the me line designated as “neutral 
of beam” completely defines the elastic deflection curve. = 
: _ are obviously dependent on the position of the mass and, hence, on time as wel oe 2 

iim 
— 
har 

‘Under the first category, the system itself is modifiedso thatan exact anal- 
is obtainable. The system of a massless beam anda single unsprung mass 
— 

> 
= 


3 ported b ene with uniformly distributed mass subjected to a moving force oe 
ee #3 load 1 devoid of mass), a system with (an infinite number of degrees) of ade q 
dom, was analyzed by Timoshenko. 5 In this case, the method of 
; 7 ‘ normal modes mentioned previously does work. Although an infinite number of _ 
q _ differential equations are involved, each equation can | be solved exactly, because 
it is linear and contains only one dependent variable. The complete solution is, "| 
‘The second category of simplifications includes those which prescribe the 
beam a given shape at any given time during the passage of the load; this, in. = 
effect, reduces the number of degrees of freedom of the structure to one. = i 


Inglis6 solved _ problem of moving mass J making use of the assumption 


its instantaneous static deflection curve due to the moving load (the deflection 
curve if the load on the beam at that instant were e applied statically). _ The case zy 
of a moving load consisting of an unsprung mass supporting a sprung mass has 
‘been considered by Tung et al. 8, and by Biggs et al. 9 who had respectively ie 
used Hillerborg’s and Inglis’ assumptions in to the of the 
“an All of the work cited dealt with the case of a single-axle line In the case 
of multiple-axle loads, if the moving loads are massless, the principle of super-— 
ie "position applies. That is, the total effect of several massless moving loads on 
% the beam is equal to the sum of the separate effect of each load. This super- __ 
-_-: position, however, does not apply if the loads possess mass. Hence, the — i 
ton of aproblem involving multiple moving masses cannot be obtained by com-— 


did Looney, !9 using a numerical proontene, has obtained solutions of the prob- 
lem involving multiple unsprung masses. - Hillerborg? has also indicated that the p 
of multiple loads could be handled use of his theory. 


£2 ng the ; same ‘assumption regarding -the mode of vibration of the beam as Inglis 


~ stant and the weight of the beam itself pata statically applied to the beam. _ or 
a It is an improvement of Hillerborg’s assumption. Consider the case of : a a 
_two-axle load with an axle spacing greater than the span length. During the 
interval after the departure of the front axle, but before the entry of the rear 
axle, there is no load o on the beam. According to Hillerborg’ Ss assumption, | the a 


_ 5 Vibration Problems in Engineering,” by S. Timoshenko and D, 3d. Bay 
“A Mathematical Treatise on in CE Binge, 


“Highway Bridge Impact Problems,” by T.P, Tung, L. E. Chen, 
N. M. Newmark, Bulletin 124, Highway Research Bd., 1956.00 
_ 9 «Vibration of Simple- Span Highway Bridges,” by J. M. Biggs, ‘iW s. Suer, and J. _ 
Louw, Transactions, ASCE, Vol. 124, 1959, 291. 
—— “Impact on Railway Bridges,” 


4 
pe 
Ricrere _ that the (dynamic) deflection curve of the beam at any time is proportional to iim : 
a 
— 
the present analysis is that the dynamic deflection curve of the neutral axis 
| jj. ofthe beam is of the same shape (not necessarily in actual values) as that — 
Hillerborg, Inst. of Structural Engrg. and Bridge Building of the Roval Inst. of Tech 
- — 


of course, does not exist ‘under the present a assumption owing to the ‘inclusion e . 4 


of the weight of the beam itself in the determination of the basic shape of the 


_ Another feature of the analysis reported herein lies in the considerations — 
- a more general type e of moving loads, that is, loads involving rotary inertia, | 
and possible unevenness of the surface of the approaches and the — 
In the following sections of the paper, the differential equations that govern > 
the behavior of the beam-load system are first derived. This is followed by a | 
a 7 listing of the dimensionless parameters of the system. After an explanation — 
of the computations of the dynamic responses, some numerical results are 
presented. Among them are those pertaining to the effects of axle spacing on 
unevenness. It is shown that a synchronization of the applications of 
x axle loads with the fundamental frequency of vibration of the beam could aug- ; 
_ment the dynamic response inthestructure. Resonance phenomena due ‘to sir sinus - 


Notation. The letter ‘adopted for use inthis pa paper are defined where 


DERIVATION OF DIFFERENTIAL EQUATIONS 
‘The general precedure of deriving the differential equations consists of first - 
a writing out the energy expressions for the and d then 


_ a with the beam at all times. The usual beam theory is used, , and the effect 
rotary inertia on the motion of the beam is neglected. 
“Profile of Travel Surface.—For a given point on the surface of the beam or — 
Zz the approach, the term “profile of travel surface” is defined to be the vertical __ 

distance measured from ahorizontal line throughthe end supports of the beam T 


is in its unloaded position (not even loaded by its own and is 
“thedeflection curve of the neutral axis measured with reference to the unloaded 
o of the beam. All three quantities in Eq. 1 are considered positive in = 

the downward direction. It should be noted that for any given n problem, whereas 
ei a is an unknown function to be determined from the analysis, the function p 

Dynamic Deflection Curve. of the previously mentioned assump 


tion regarding the deflection = — axis, Eq. 2 may 


= 


time, but not 
af “4 


— 
“a 
— 
— 
— 
— 
q 
— 
— 
qq 
— 
in which the symbols w(x,t) denote the fact that w depends on both the length 
+ 
a a coefficient of proportionality which may v _ 
e; Y(x,t) denotes the deflection curve of theeutrar a 


axis if the of the beam a the ween were stz tically applied to 
the beam. The function Y depends on time t (as well as on distance x), because vale 
the positions of the loads depend on t. Letting the time t = 0 correspondtothe — 
instant of the entry of the rear axle to the beam, and xj and x2 denote the dis- ; 
tance from the left support to the rear and front axle, respectively, one sees — 

_ that xj =vt, andxg=vt+s, , where v denotes the horizontal velocity of the 
loads and s the axle spacing. Since both v and s are constants, for a given i 
the function Y may be computed by satiate In other words, it is a oe func- a 

ie function Y can be written as a polynomial (this is done i Annentin D, 
it can be expressed as a series}! as 


The subscript, jr refers to the rear, or front axle according as it is equal .. 
 lor2 _ Wj denotes the weight of the jth axle; the symbol €;, defined as: €j= e 
7 for 0<xj <L, g = 0 otherwise, is introduced to ensure that any axle load out- 


‘The s summation in n the second bracket extends cover odd values es of n only, oe 
cause the static deflection of the beam, due to its own weight, is cymmetrical a 
In order to avoid writing lengthy | in later 

may be written in the following form 


Tat 
n 


nn sin 22 
and, for odd values of Wits p58 


in which | Ko is some constant depending on the choice of Kj. saad even values — 
of fn, the second term inside the brackets in | Eq. he should bi be ® dropped. Pee 


_ 11 “Advanced Mathematics for Engineers,” by H, W. Reddick and F, H. Miller, John F 
Wiley and Sons Co., New York, 3d p. aA. 


> 
aA 
| 
— 
— 
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- variable that specifies the dynamic deflection curve since W is aknown function. 4 . 
pe = Potential Energy. —The initial energy level of this system is defined by the — 
- - following : conditions: (a) The beam is not strained (not even by its own weight © 
‘ ; as yet); (b) the weight of the beam andthat of the unsprung loads are at a hori-— 
zontal line having the same elevation as that of the unstrained beam; (c) the 

ae springs of the loads» are undeformed (not even by the sprung weight as yet) :: 

a -_ and (d) - ‘sprung weight is situated at the upper ends of the undeformed 


we the 


each x a differentiation with respect to x. 
— For alinearly elastic spring its strain energy is equal to 1/2 (k d2), where 
denotes the spring constant and d the deformation. For in t 


The the 
vertical a mk of the centroid of the aprung 2 mass ‘(positive if downward); 
ae u denotes the rotation of the sprung mass about its } centroidal axis (positive if if 
* clockwise) ; a; is the horizontal distance from the centroid of the sprung mass” 


& Wj a2) 


which M denotes the total sprung mass, the deflection of the centroid 
M when the sprung weight | is statically : saehied to the § springs, and mj Ss 
_ the unsprung mass of the jth — 
16,11 11 and 12, and noting the following identities: 


4 
— 
— 
— 
— 
— 
Wa 
&g 
al 
tates ror te denotes w evaluate 
tU eight static energy of both sp — 
ees pring constant ofthe jth axle, and the sym 


RESPONSE 


Kinetic —The Kinetic energy of the the 


"energy of the beam; the persia and the third terms represent the kinetic en- 
ergy of the sprung mass due to the vertical and angular motion, respectively; 
ann the fourth term represents the kinetic energy of the unsprung masses. _ 


tem are obtained by substituting Eqs. (14) and (15) into Lagrange’s samiitent bes) 


tt 


¥ (-1 


| 
6 
— a 
which the On with respect to time t, and = 
| 4 ri represents the polar moment of inertia of the sprung mass about its centroid. ni Nal 
— 
— 
— 
—— 
— 
— 


y n neds! for which the terms due tothe weights 


 the“fastest” converging series is 


the beam the loads decrease in as fast and n 
est” conver erging ser 


Since this series and those s series s involving Mn in . Eq. 17 converge rapidly, to Ge : 
facilitate numerical computations, it seems reasonable to keep only the first | o 
_ term of each of these series. , Consequently, Eq. 17 = reduced t vo the the following — 


we 
j(z+(-Dba, 


Ww; 


in in which M;- iy is called the “sprung mass of the j th axle” ” and n = 1. 


am It may | be noted that Eqs. 18, 19 and 20 constitute a system of three linear — 
differential equations with coefficients that are functions of time. 
Dimensionless Parameters and Variables.—To obtain the dimensionless 
‘parameters and variables of the system the preceding differential equations _ cS 
could be transformed into dimensionless form. This may be accomplished by 


first Eqs. 18 20 by the quantity my g L, and Eq. 19 by the 


2 Tp/L) M aj’, in which h Tp, denoting t the fundamental natural period 


of the beam, is equal to (2 L2/n) )/mp/EL Next, let the following dimension-— 
1 fundamental natural period of 


1 axle spacing 
mass of j th axle 


1 mass of bea 
tota [mass of beam 


_unsp rung mass of j th axle 
total mass of 


ro lar moment of inertia about 
4 =Mj sum of moments about centroid 


requenc of j th 
“(Frequency of beam 


4 

afore-_ 
__mentionedterms decrease in magnitude atthe rate of and n-4, respectively, 
i 

— 

.. 


DYNAMIC RESPONSE 


deflection due a Sores equal to to ‘weight 


Dependent Variables of System. . 


‘rotation of spr mass x a spacin 
distance of rear — left suppor 


23a) 
l 


that the differential quantities, be written fol- 


in which each subscript T represents a differentiation with respect to7, and ‘ 


the letters A, B, and C (with subscripts) represent known functions of the di- = 
-mensionless parameters and the independent variable T . These functions 
COMPUTATIONS OF DYNAMIC RESPONSES 


of the neutral axis, is measured from its original position when the beam is. 
. not even loaded by its own weight. However, in the case of such practical en- 7 
gineering problems as the dynamic behavior of highway bridges under moving © _ 
vehicles, (for which the present beam-load model was devised), the weight of c 
the beam (representing the bridge) must be considered as already applied to 
ra the structure at the start of the problem. - is this case that is considered in 4 
3 Before describing the computations of response, it may be pointed out that - 
_ for the prior- ‘-mentioned case the initial value of y, for the integration of the 
"differential equations, should be set equal to the static deflection due to the e 
weight of the beam. Furthermore, if the beam is cambered in order to have a 
q ¥ 4 a surface under the action of its weight, the initial profile function 
p(x) must be set — to the negative | of the static one due to the dead 


weight. 


' 


— «| 
4 
(222) 
— 
— 
(20) 
— 
— 
— 
. 
— 
In the derivation of the differential agquationse the quantity v 


ber, 1960 EM 
may be noted also that the expressions of given subsequently 
or do not include the static effects due to the weight of the beam itself. ee 
Dynamic Deflections .—Denoting by yg the dynamic due to the 


Plications of the axle loads, it is seen — 


Dynamic Bending Mome vts.—The | dynamic ‘bending moments are e computed 
ee . treating the dynamic reactions between the axles and the beam surface and an 
the inertia forces of the beam as forces statically applied to the structure. ae 


Denoting by the reaction between the jth axle and the beam | sur- 


first term in brackets represents the weight of the axle; the ‘second ie 
represents the dynamic force in the ae ageing; and the third term the inertia in 


‘The dynamic ‘shear at x due to the inertia of the beam 


ix. ax... (30) 


represents the reaction at the left support; 
‘second | integral represents the total load to the left of section x. It has been — 


noted previously that y= = sin , and the series in Eq. 
rapidly. Since the inertia forces | don not constitute the primary part of the 
= loading when the total maximum dynamic bending moment is considered, it 


L 


q 


The total eaiuatcas bending moment at section x, (Max, is then equal to the 


i} 
| 
— 
— ay 
| 
5 
{ 
The bending m beam is then given by the 
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mentioned previously the governing differential equations, s. 19, 
d 20 or Eqs. 24, 25 and 26, may be integrated numerically. | The method of — 
integration used for t the solutions presented hi herein is duet to 


_ The numerical problems considered here pertain to the case mentioned 4 
previously: that is, as an initial condition of the system, the weight of the beam hi 

has already been applied to it. Furthermore, the response of the beam is con- 
sidered for the time interval between the entry of the front axle to the span and “a 
the departure of the rear axlefrom the span. Unless otherwise stated, the re- 
sults are obtained for the following conditions at the instant the front axle en- 

beam: (a) The aprung part of the load has no vertical or angular 

ti and (c) the beam surface profile is perfectly 
In the following, the dynamic response are presented as values of the ratio_ 

_ of the dynamic response to certain reference static response. ih each case, 

_ the reference static response used may be computed from theories of deter- : 
The solutions were obtained with a computer program! 3 prepared for for use on 

the ILLIAC, the digital computer of the University of Illinois. noe aS 

_-«- History Curves.—Ahistory curve is agraph of the response of the beam a 
a function of time. Fig. 2 shows a history curve of the bending moment at the | 
mid-span. This is due to atwo- axle load in terms of the maximum static bend- 

ing moment at the same section of the beam, due to the same load. ‘The hori- i 

_ zontal axis in this figure denotes the dimensionless quantity v t/L, which rep- : 
7 resents the distance over which the rear axle has traveled on the beam asa | 
- fraction | of the span length. The numerical values of the parameters of the 
problem ‘solved are given as follows: a= 0. 15500; 8/( 30110; 
12600; 0.11820; = 0.01498; 0.02282; = 0.30400; op = 0. 27000; 
These numerical values correspond to a 45-foot highway 
7 4 bridge and a 20 ton two- axle vehicle with an axle spacing of 30 feet traveling — : 
Also plotted in Fig. 2 is a curve representing the static response, that is, 
“the response when the speed of the moving load approaches zero. It may be : 


tha about 15% larger than the maximum static moment. — res 
the same as s that considered in Fig. 2, the dynamic and static 


12 “A Method for Structural Dynamics,” by N. M. Newmark, Journal 
we of Engineering Mechanics Division, ASCE, Vol. 85, No, EM 3, July, 1959. ans eae 
a peed 13 “Dynamic Behavior of Simple Span Highway Bridges Traversed by Two-axle Ve- *) 
hicles, R.K. Wen, thesis presented to the Univ. of Dlinois, at Urbana, in 1957, in partial 
; fulfilment of the requirements for the degree of of 


curves (in terms of maximum static mid- 
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2, HISTORY CURVES FOR ‘MID-SPAN BENDING MOMENTS. 
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4.—MAXIMUM RESPONSE AS FUNCTION OF AXLE SPACING PARAMETER, 
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MAXIMUM DYNAMIC BENDING MOMENT 


ABSOLUTE MAXIMUM STATIC BENDING 
= 


= 


— 
— 
— 
— 
— —_ 
| 
te 
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span deflection) are shown in Fig. 3. In this case, the maximum dynamic call 
__ flection is seen to be about 25% larger than the maximum static deflection. _ 
«Spectrum Curves.—A “spectrum curve” is a plot of the maximum dynamic | 
response against some parameter of the beam- “lens system. Hence, each point 
ofa spectrum curve represents the maximum ordinate taken from some his- | 
_ tory curve. Unless otherwise stated, the numerical values of the various para- 
meters for the name presented below are the same as those previously, 


_ for mid- span bending moment as ‘functions of ‘the « axle spacing parameter” 

s/(2 L), for three values of the “speed parameter” 
~~ seems logical to expect that when the applications of the moving axle 

loads on the beam are synchronized with the vibration of the structure, the — 

dynamic response in the structure may be augmented. This appears to be born o 
out by the results shown in Fig. 4. The condition of the afore-mentioned = 

_chronization may be represented equation: 

| 


which n is an integer, or by the equivalent, - Se 


_ It can be seen from Fig. 4, that for a given curve the peak rer responses occur 
at values of s/(2 that are very close to multiples of the value of to 
which the. curve corresponds. For example, take the curve for a = 0.105. The = 
‘first, second, and third peaks ‘Occur respectively at values of s/(2 L) equal _ 
Lene ‘Fig. 5 shows two iain curves for bending moments at sections | in- 
stantaneously underneath the axles as functions of a. In this case, the maxi-— 
7 mum dynamic bending moments are presented as values of the ratio of these = 


general shape of these curves are e similar to those obtained for the case 
, of a single axle load.’ They indicate that the values of the maximum reapense 
an “oscillatory” appearance, but the amplitude of the “ oscillations” in- 
creases with an increase in the values of the speed parameter a. 
‘It may be observed that not only the peak values of the response se for the 
_— underneath the rear axle are appreciably larger thanthose for the front — 
axle, but the former response is also more sensitive to the variation of the 
parameter a. ‘This behavior may be explained by the fact that the beam has 
already been excited by the front axle before the rear axle enters it, whereas ‘ano 


it had been at rest when the front axle first came onthe span. 


in ‘which ba and m m denote, ly, the amplitude and the ‘number of halt 


sine waves representing the surface of the For the problem 
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DYNAMIC RESPONSE 105 
in Fig. 6, the value of m, represented by ‘the horizontal a9 axis, is varied, while — 
b is held at a constant wanes of 0.0032 g hw" . The values of the other para- Ps 
meters for this problem are: a= 0.10000; s/2 L = 0.10000; Ty = =0.07184; 
= yg = 0.01112; 04 = = 0. 82180; and = = 0.25000. dig 
ae it may be noted that the graph in Fig. 6 has the characteristic appearance 


of a “resonance diagram” for some oscillating system. It shows thatthere are 


two “critical” values of m: m= 8, 5, and m = 12.5, at which the responses at- 


: tain a relative maximum. It is of interest to note that these values may be pre- | 
dicted approximately by a simple consideration of possible “resonance” be- 

Picturing the crossing of the load over the wavy surface of the beam as a io 
|p load resting on a horizontally moving surface, and neglecting the a 
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OF HALF SINUSOIDAL WAVES REPRESENTING | 


‘FIG, .6. RESPONSE AS FUNCTION OF NUMBER OF HALF SINUSOIDAL 


‘ “flexibility of the beam, the load may be thought of as being subjected a a peri- = 
odie disturbance with a period equal to (2 L)/(mv). _ Setting» this quantity 
- equal to the natural period of vibration of the moving load a“critical” value of 3 
It can be shown that if o1= which is the case for the system: being 


— | 
— 
: 
— 
f angular oscillation, he fundamental natural 
| and that o (m v) equal to t : ained as follows: 


surface. 


beam subjectedto a moving two- ‘The analysis is based on 
_ _ assumption that the shape of the dynamic deflection curve of the beam is sim- 
‘moving loads. It may be noted that this } assumption would lead to an exact solu- 
a _ tion of the problem if the s speed of the moving loads approaches zero, andthe 
load system is initially atrest. 
‘The question naturally arises regarding the validity of this method of anal- 
«ysis for | finite load speeds. . On this question, it may be mentioned that. an 
analysis, developed for the case of a single axle load based on the previously 2 
mentioned s assumption, yields results that agree very well | 
_ experiments. 7 Since, as discussed under the heading, “General Description 
of Present Analysis,” the basic assumption used herein is a more rational one 
~ than Hillerborg’s, it would seem reasonable to expect the results obtained a 
the present analysis to be as accurate as those obtained for asingle axle load > 
_ Among the few numerical solutions presented, the results showing the effect fe 
_ of the variation of axle spacing seem to vindicate the belief that a synchroniza- a 
_ tion of the successive applications of axle loads with the fundamental period al 
- vibration of the beam would augment the dynamic response of the oye 
_ In this paper, only a single unit of two-axle load has been considered. - 
pi damping in the system has been assumed negligible. However, it would not - 
difficult to consider in the analysis multiple units of two-axle loads and damp- 


ing effects. In fact, the differential ae. the motion of such a 
system have been derived elsewhere.12 


_ This paper is basedon part of a thesis, by tl the author, submitted to the Grad. a. 
. College of the Univ. of Illinois, in partial fulfillment of the requirements for the : 
degree of Doctor of Philosophy in Engrg. . The thesis was _wrtiten under the a ee 
_ supervision of Professor A. S. Veletsos. 
a The investigation was sponsored by the Illinois Division of Highways at a a 
. the U.S. Bureau of Public Roads. The statements made in this paper, ee 
y do not necessarily represent the views of these organizations. 
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7 or large values of m (greater than about 16), the response a, er a * 
- ractical purposes, the same as that for m = 0, that is, for a smooth beam ae =a ee 
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APPENDIX COEFFICIENTS oF DIMENSIONLESS DIFFERENTIAL 
| 24, 25 and nd 26) “ 
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‘The functions Wj (in polynomial form), n, and their derivatives that soar 
in the preceding are as follows: 


“in which Yp a and Yj denote the static deflections due to the weight of the beam 7 ‘ 


d the j th axle, respectively, Eq. 39 may then be written 


woe 


Choosing 


24 EI 


y 


— 
gg 
a 
— 
be derived in a similar manner. 
function yj can d its derivatives, evaluat — 


in which the subscript k, representing an axle, is a “dummy” dex used in 
“summation; and represent Xx; L) and 1x, /( 2L),  Tespectively. 


«Eg. 41 holds true only w when < For the case in which 
symbols and Should | be interchanged. be noted that 42. 
and 43 are symmetrical n €; and 
 Then- function is by Eq 6. In to with the pre-- 
+ vious choice of the value of Kj, the value of f Kp in Eg. | 6 ‘should be | be set equal t to oo 
— 
The first. and second derivatives of are, respectively, 
j . 


The following symbols were adopted for use in the paper ont for 


function of time defined in 7: 


= gravitational acceleration; 


= polar moment of inertia of M about asian axis; ea 


= axle ir index; j=1 2 tothe rear rear “and: front axle; 


spring ‘constant of jth axle; 


— 


=total sprung mass; 


7 

= 

| 
| 

i 

— 

ye 

— 

— 

— 


ass of jth 


 =aninteger; 


= ef initial profile of travel surface; 


= kinetic energy; — 


= fundamental | period of f beam; 


ential energy; 

=p+ 


coordinate, see Fig. 


a distance from left support t to - th axle; 


.. static deflection of beam due to its own m weights 
Mection axis of beam; 


mins = . static deformation of spring in in jth axle; 
= defined in Eq. 21; 
i) defined in Eq. 21; 
j = defined in Eq. 21; a 


defined in Eq. 
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‘BEARING C CAPACITY OF FLOATING I ICE 


By G.G 1 


collapse of floating ice is estimated on the basis of 


a rigid-plastic plate on an elastic foundation for various loading and boundary = 
conditions. The results are compared with field observations in the epee PY q 


; _A large floating ice sheet under load can be represented by a plate on an ‘ . 

- elastic foundation of which | the reaction pressure at every point is directly pro ‘ 
u - portional to the deflection of the plate at that point. This approach using athe 
Be 2 - ory of elastic plates in the problem of the bearing capacity of ice sheets was 
‘ first made by Hertz some 75 yr ago (1884) and has been used ever since. In- 
vestigations have shown that the estimated strains and deflections under short-_ 


_ time loads, much less than the ultimate, agree well with observations. How- ‘ 
ever, as the failure load is approached, or under long- time loading conditions, : 
the of ‘becomes less applicable. It is believed better 


Note.—Discussion open until 196i. Pe the closing date one 
i a written request must be filed with the Executive Secretary, ASCE, This paper is = 
3 f the copyrighted Journal of the Engineering Mechanics Division, Proceedings o of the 4 
of Civil Engineers, Vol. 86, No, EM 5, October, 1960, 
of Civ, Engrg., _ Nova Scotia — College, ‘Halifax, Nova Scotia, Can- 
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estimate of the ultimate load can be obtained from a limit analysis on 
a CENTRAL LOAD ON LARGE ICE S SHEETS (INFINITE PLATES) © 


Circular Load.— -When a central concentrated load, much less than the ult. ae 
; ~ mate, is applied over a small circular area on a large floating ice sheet of P18 
_ uniform thickness, the stresses and deflection of the sheet can be computed as 
for a thin, elastic, and infinite plate on a discontinuous- -elastic foundation 
(Fig. 1(a)). The maximum positive radial and circumferential bending mo- — 
be ments per unit length, My and Mt respectively, are equal under the center of _ 
_ the load, and the corresponding principal tensile stresses at the bottom of vad of the — 


st = 1.1 (1 + u) (log L/a + 0, 26) P/h2 


in which L is the characteristic length or radius 0 of ‘relative stiffness, and 


is the of ray of the plate, is applied load, a a isthe con- 


dation modulus and is Poisson” Ss ‘ratio of the plate. The maximum 
— bending moment occurs at a distance of about r = 1.9 Lfrom the load 
is, Mr = ~0.02 P. The maximum deflection 

ine the of deflection dish, below the evel sur is 


a 


terms of the maximum or ‘yield moment per unit length, 

= fp h2/6, where fh is the flexural strength of material. In the limit 

a point load (a= 0) the theoretical yield load is zero. Since the maximum neg- wee 
s . ative radial bending moment is small compared with the moment under the a 

= x ences. elastic ice sheets may be classified as infinite when they have a oe, 


radius R23 L. Further, a concentrated load would have to be more than 3 ee 


With increasing load the » radial cracks in the ice sheet increase in length 
and number until the bending stresses along a circumferential section of the 
x sheet become equal to the flexural strength of the ice anda circumferential " 
tension crack is formed on top (Fig. stm). ghee ultimate bearing capacity or 


ie 
7 
4 
[| 
“ 
a4 
4 
estimate rom 4 and is Shown in Flg. 4a ior various ratios OF 
: 
= 
4 
_ cal collapse load it is assumed that the ice sheet is a thin, rigid, ideally- _. tm 
| and infinite plate, which can, without cracking, resist a full plastic or 4 
2 Numerals in parenthesis, thus (1), refer to correspondingitems inthe Bibliography a 


ICE SHEETS 
limit — moment per wma? length, = fo 2/4, in n which fo is the yield 
2 strength of the material in tension, At collapse radial yield hinge lines extend © 


_ from a plastic section below the load to a circumferential yield- hinge- -circle — 
aia The load carrying capacities of simply - supported ; and fixed i ideally- plastic 
_circular plates under central concentrated and uniformly d distributed ed loads 
have been discussed by several authors (5),(6),(7). method of analysis 
; = can be extended to plates on a discontinuous-elastic foundation by neglecting, 
as previously, membrane action and shearing deformation. Full details are- 
a given in Appendix I, for a concentrated load on an infinite plate of which the 


at material obeys the Tresca yield condition (shearing strength equal to one- 
half of yield strength in and the flow rule. theoretical 
A 


ps 


hown in Fig. 3(a) for various ratios of a/b . While local — 

‘bearing capacity under the load governs the minimum contact radius a, | as 
shown or the collapse | value for a point load is found» to be independent well 


te 
— 
— 


_ The radial deflection of the plate is in the form of an inverted cone near he ’ _ 
load and passes into a logarithmic shape at greater distance. 
‘The radius b in the present analysis can only be expressed in terms of th e > 
characteristic plate parameter L by making a complete elastic-plastic analy-— 
sis. _ When a concentrated load exceeds the yield load, the plate material _ 


the load passes from the elastic to the plastic state and the central deflection 


load and the total reaction pressure below the plate is ‘maintained, the radius 

of the deflection dish at collapse is less than the value of 3.9 L of an elastic a 

‘Plate. Similarly, the position of the maximum negative radial bending moment 45 
= to move towards the load so that at collapse the critical circumferential _ ; 
section with the yield hinge circle is likely to have a radius of less than the os 
me 9 L obtained in the elastic state. Within the above limits, an effective rates 5 


_ increases more rapidly than the load. i. Since overall equilibrium of the apled a 


sidered subsequently) as an an approximate and safe value in estimating —— 
lapse load from Fig. 3(a). The corresponding results are shown in Fig. 4(a) _ 


for) various ratios of a/L. While Eq. 4 is unaffected, the collapse load from — 


3.39 (1 + 3a/2 L) Mo (for 0.05<a/L<1) . 
in which L is given by Eq. 2. "Since the yield moment, My, of ice 
equal to the limit moment Mo, comparison of Figs. and 4(a) ine 
_ dicates the theoretical collapse load of large, floating ice sheets ‘under a cen- 
ore tral concentrated load is, at least, approximately twice the load at which — 
* ing begins or the first cracks form. 
‘Strip Load.—If a number of equal central concentrated loads are applied in 
ee on a large floating ice sheet, the analysis is simplified by ignoring longi 
_ tudinal bending and using an equivalent strip load, P, per unit length (Fig. 
1(a)). In the elastic state the maximum positive transverse moment, 
per unit under the center of the load is 


ne is ine half- -width of the load and L is given by Eq. 2. The maximum > 
negative transverse moment occurs at a distance of about 4 = 2.2 L from the 
load and is, approximately, My -0. 07 PL, The maximum is 


ko and the half-wi width alae the deflection di dish is about r= = 3. 3 L. The yield load at 


¥ ong which the first longitudinal tension crack develops at the bottom of the ice 
- gheet below the load can be estimated from Eq. 7 and is shown in Fig. 2(b). 
_A-strip load has to be more than 3 L from the edge of ead sheets to satisfy 
the above equations and develop the estimated yield load. 
he As the load is increased, the longitudinal crack in the ice sheet yong 
until a longitudinal crack on each side of the load is formedon top of the sheet, 
and the collapse load is reached. At collapse, longitudinal yield hinge lines 
7 run along the center-line of the load and on each side of the load at a distance _ 
“ from the axis (Fig. 1(c)). A solution of the collapse load has been obtained 


— elicadl I and is. ‘shown in Fig. 3(b) for various ratios of a/b. While the 
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_— bis likely t to be between the limits of. 2. 2 Lan and 3. 3 . a pre- 


aw 1960 


4 
viously, latter the collapse load per unit length is, ap 


The results n in Fig. 4(b) for various of a/L and it is found 
_ that the theoretical collapse load is of the same order of magnitude as the 
yield load (for My = , as before . . Comparison of the present collapse loads: a 
with those for concentrated loads indicates that point loads spaced more than a 
- about b/2 or 2 L apart should betreated as individual loads, while foragreater = } = 
contact width the governing distance is about 3 Lto 4 L. The transverse de- 


4 


- flection of the plate at collapse under a strip ~_ is in the form of an inverted 2, 


ae aie LOAD ON ICE SHEETS OF LIMITED EXTENT e+ 


ae _— Load.—The size of floating | ice sheets is sometimes less i that 


required to develop the limit moment, Mo, along the ‘circumferential section 


at r=b (Fig. 1). . An analysis of simply supported ideally-plastic circular 


of R = 2.7 b to develop the radial limit moment. Because it was indicated pre- _ 


= plates with overhang shows (6) that plates must have a minimum total —_ 


&g ——— Viously | that the radius of the critical circumferential section is likely to be © 


less than 1.9 L, large plastic ice sheets may be classified as infinite when _ 
they have a radius R >5 L, which is nearly twice the radius suggested for - 
large elastic sheets. ee concentrated load would thus have to be more than 5 


4 from the edge of large sheets to develop the collapse load of Eqs. 4 and 6. 
_ As the plate size decreases, the maximum stress under a central concen- 
trated load on an elastic plate decreases below that given by Eq. 1 so that the 
_ theoretical yield or first cracking load of an ice sheet ‘is then greater than 


that of an infinite plate (2). A solution of this problem for the collapse load — 


eS of a plate of medium size (2 < R/L < 5) is given in Appendix I andthe results — 


are shown in Figs. , 3(a) and 4(a). The theoretical collapse load decreases — 
with the plate size (expressed as ratio of R/b or R/L, using | b= =1.9 L - 


ee, with afurther reduction of the plate size the maximum stress | under a cen- 
al concentrated load on an elastic plate decreases, while” ‘small plates 
el with a radius of about R < Lthe reaction pressure is practically uniform (2). 
On the assumption of a uniform reaction pressure below a centrally looted 
plate of small size (R < L) the collapse load has been determined in 
land the results are shown in Figs. 3(a) and 4(a). . The theoretical — 


7 aaa way to that of a | medium plate, but it is generally smaller, except for bs 


a a point load when Eq. 4 applies. The collapse load for a larger contact radius _ ut 


and ee 9 Las larger plates is, approximately, 


(or 0.41 <a/t< (10) 
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Strip ‘Load. ~In onder to develop the limit moment, Mo, along the longitudi- — 
nal section at r = b under a strip load, the cantilever overhang of the plate be- 

_ yond that section must have a minimum length of to, so that a large plate must 
_ have a minimum total half-width R=b+to (Fig. 1(c)). Since the critical 

a section is likely to be at a distance of b less than 2. 2L, as pre- _ 
viously indicated, a safe estimate of the minimum half-width is, after substi- 


2.2L 
in which q q is the unit weight of the material. 


fi strip load would have to be 
-more than the distance given by Eq. 11, from the ono of large sheets to de- 


_ For a smaller plate width, the yield load increases while the collapse load 
decreases. _ The collapse load of a plate of medium size (0 <t <to hasbeen 
_ derived in Appendix I and the results are shown in Figs. 3(b) and 4(b) using 


b= 3.3 Las for the larger plates. | The collapse load per unit unit length is, ap- = 
proximately, 


LOAD ON LARGE SHEETS (SEMI-INFINITE 
Circular: Load. 


— When a large floating ice sheet is the bearing 
capacity is reduced and it is assumed that a crack cannot transmit bending — 7 


moments or shearing forces. For a concentrated load applied over r a small 
circular area tangent to the edge of an elastic semi-infinite plate, the maxi- _ f 
mum tensile stress at the edge under the 
sr = |8.8 lo + 0.09 - 


of about 2 L from ‘the load and is, approximately, My = -0. 075 P. ‘The cor- 
as responding maximum moment at right angles to the edge occurs at a distance 


«of of about 1. 4 L from the load and is, approximately, Mr = -0.04 P. ‘Both mo- 

s ments are found on about the same deflection contour of the plate, as shown by ; 
model tests on elastic plates (4). The maximum is 

‘(1 0.89 a/L) P/2.1 k L? 


and the radius of the deflection dish is about r=3.9L ‘along t the edg edge e 


hag 4. 9 Lat right angles to the edge. ‘The yield load from Eq. 14 is shown in | 


a distance of the load trom. ‘the free edge of the sheet et increases, the yield wal 


— 
— 
— 
— 
— 
Po (for 0 < a/L< landO<t < to) (12) q 
For small plates with R< L, as before, the deflection is approximatelyuni- _ 
| a ' _ form, The corresponding collapse load has been derived in Appendix I and is ae 
er 7 shown in Figs. 3(b) and 4(b) for b = 3.3 L. The collapse load per unit length — oa 
a3) 
— 

| 
4 


j Fe. of the sheet in the central zone and spread outwards until a pt dont 
crack is formed on top of the sheet extending to the free edge thereby 
leading to collapse (Fig. 5(a)). state of plastic equilibrium and the de- 
_ flection shape in the central sector are similar to those of an infinite plate * 2 
with radial yield hinge lines anda circumferential yield hinge arc. Within the 
er, two outer zones, however, the radial andcircumferential bending moments are 
less than the limit moment, which obtains only along the circumferential sec- ‘= 4 
tion along yield lines continuous with the yield hinge arc of the central sector. | 
An approximate analysis of this problem is given in Appendix I and indicates 
that at collapse the outer zones are roughly triangular with an angle of 45° 
along the edge, and that the central sector has a center angle of 90° (Fig. 7 
st 5(a)). The results shown in Figs. 3(a) and 4(a) for a semi-circular contact 
~ area of the load can be applied to a circular area by using an equivalent con- = 
tact radius of The theoretical collapse load of an edge- loaded plate is 
thus about one-half of that of a center-loaded En | 
With increasing distance of the load from the free edge of the sheet the col- 
_ lapse load increases and approaches that of a central load when this distance 
is about 5 L. = An approximate solution of this problem is given in Appendix I 
by extending the analysis of an edge-loaded plate, which indicates that at col- 


lapse the central sector lies — two outer zones extending 
the free (Fig. 5(b)). 


Ja + 1+ e/4L) )Mo 


n= 17 show that at the upper limit of e = 5 Lthe hecseuabed collapse load is 
- about one-third greater than the value for a central load from Eqs. 4 and 6. 
‘ This difference is due to the approximate (uamer bound) nature of the analysis. 
|. a Strip Load.—If a strip load of width 2 a is adjacent to the edge of an elastic © 
= semi-infinite plate, the maximum negative transverse bending moment — 


at a distance of about 1.1 L the and is Ber 


7 and the width of the nettinthins dish is about r r= 22 2 L. 


The yield load derived a 
from Eq. 18 is shown in Fig. 2(b) and is roughly the same as for a center- 


a collapse of an edge-loaded plate when the deflected shape is a plane, a a 


longitudinal yield hinge line develops at a distance b from the load. The col- 
_ lapse load has been derived in fesse I and is shown in Figs. 3(b) and hae 


| 
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ay 
hese estimates may be seamen with th the solution using elastic theory (Eq. 


elastic- plastic analysis ‘and exceeds that estimated from. Eq. 20. The dif- 
ference decreases rapidly with increasing contact width and is reversed for 
gente widths (Fig. 4(b)). Hence, the theoretical collapse load of an edge, 
loaded sheet is generally above one-half that of a center-loaded large sheet. a 
Comparing the present collapse loads with those for concentrated edge- loads, : 
_ it is found that the latter should be treated as individual loads if they are ry 
mere 
Ss With increasing distance e > of the s strip load from the free edge of the sheet = 
_ the collapse load increases 3 and approaches that ofa central load at e = R (Eq. 
11). A solution of this problem is given in Appendix I and, using a safe value 


a b= 3 L as for a central load, load unit length is, approxi- 


approaches the value for a central load (Eq. 9) at tates 


CORNER AND TIP LOADS ON LARGE ICE SHEETS 
(WEDGES AND SECTORIAL PLATES) 
: te a load is. ‘applied at the > corner or | tip of a wedge formed by a large - 
floating ice sheet that is cracked along several radial sections, the stresses 
_and deflection in the elastic state and the collapse load in the plastic state can _ 
be determined as for an edge-loaded sheet. For a concentrated load applied 
over a small circular area at the corner of a large elastic wine! (90° = of 
infinite extent) maximum tensile stress is (3), 


stress occurs sil it the top of he plate at a distance from the corner of 


5 - the radius of the deflection dish is about r = 3.6 6 L on the ¢ edge and about 
L along t the bisection (4). The yield load | from Eq. 23 is shown in 


— 
— 
q 


wedge or sectorial plate with an apex angle of less than 90° does not appear to 
_ available (9). However, an approximate solution can be obtained by ex- 
_ tending the analysis leading to Eq. 23, which was based on the assumption that 
_ the corner-loaded plate acts like a beam of variable width without transverse . 
bending. _ This assumption is closely supported by the results of model tests | 
on tip-loaded elastic plates with a 90° angle. It is in complete agreement with 
: tests on sectorial plates with an apex angle of 60° or less, when the deflection 
F 7 contours remained straight and the radius of the deflection dish along the bi- | 
. ector was about r = 2.7 L (4). If the original circular contact area of radius a 
a of the concentrated loadadjacent to the corner of a 90° sectorial plate is re- 
placed by a ‘a quarter- circle of equal a area with the center at the corner of the > . 
_~ (or by a similar equivalent triangle), the center of gravity of either lat- _ 
_ ter area is ata distance of about 1.2 a from the corner, compared with 2a for 
the ag contact area. that a concentrated at the tip of 


u close to that of the original circular contact area of a a plate with @ = = 90°, Ac- 
- cordingly, the magnitude and position of the maximum stress and the deflection — 
of a tip-loaded elastic wedge or sectorial plate can be obtained with little, if 

any, loss of accuracy by modifying Eqs. 23 to 25 for the apex angle (@ < 90°) 
and the distance g from the tipto the center of gravity. of the concentrated load, 

applied over a contact sector of radius a or similar triangle 
a Thus, the maximum tensile stress at the top of an infinite elastic wedge 


late with 90° is (using Eq. 23) 
— 


L, measured 
bisector Eq. 24). The distance g varies from 0.60 a for 6 = 90°to. 7 


the upper limit of 2 2 a/3 for a very small a (< 10°, say) or narrow wedge for 
which the maximum stress, from Eq. (26), ‘edhe 


1 - 0.88 


For a concentrated load applied over a small sector of aat the 
of a plate the yield load, from Eq. 26 for 8 = 90° is somewhat less than that — 
given by Eq. 23 because of the > smaller contact radius. The theoretical yield bf 
load of a large floating ice sheet cracaed radially into numerous long narrow 

and supporting a common concentrated load at their tips can be esti- 

_ mated from Eq. 27 for 6 = 27 and is shown in Fig. 2(a). This yield load is 
= twice that of a large uncracked sheet. The great difference of the theo- 


A complete ana S for 4 Confentrared (040 41 the tin oF ininite els 4 | 
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the 
43 = sector of the same radius 2 a as before, then the center of gravity of the con- $$ 7 
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distance r - , the maximum deflectionis 
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, 1960 
ay bending moment near the concentrated load on an uncracked a q 
, The e collapse load of a corner or tip-loaded long wedge or sectorial plate 
ean be determined as for a large edge-loaded plate. As the apex angle. of the * 
_ plate decreases, the outer triangular zones (Fig. 5(a)) are assumed to re- 4 
; ~~ unchanged but the central sector decreases with the apex angle, and the : 
- collapse load decreases accordingly. For the lower limit of a large | 90° wedge : a 
: é central zone vanishes and the theoretical collapse load is the same as for os 
a beam of variable width (Fig. 6(a)). The solution is given in Appendix I and 
shown in Figs. 3(a) and 4(a) for a corner-loaded plate (@ = 90°), using b = 2.7 
as mentioned previously. load for 6 < 90° is, approximately 


Mo tan 6/2 (lor 0.2 <1). 


While the theoretical load of a a corner- loaded 90° wedge is about one- 
quarter of that of a center-loaded large plate, the collapse load for 6 < 90° 
decreases in direct proportion to the width of the as as would be expected 

“= The theoretical collapse load of tip- loaded ‘itil plates" (Fig. 6(b)) is 
directly proportional to the angle 8,as shown in n in Appendix I I and, using b= 2.7 


‘These results are shown in Figs. 3(a) and 4(a) t together with an ay approximate 
_ relationship for a/L < 0.2, as indicated in Appendix A ‘It is of interest to note 
“= this collapse load is only slightly less than that of a large uncracked ~ 
_ sheet. The deflected shape of ¢ a wedge or sectorial plate with 6 < 90° at col- 

= is generally in the form of a surface. 


These results shown in ‘Fig. 4(a). While the latter ‘collapse load fora 
.- load is the same as using plastic theory, for small contact areas the col-_ 7 
_ lapse load using elastic theory is up to approximately 30% smaller than the g 3 
ns in value from the limit analysis. Similar results are obtained for 
edge-loadedand corner-loaded radially-cracked sheets, since the correspond- 
_ ing collapse loads are respectively, one-half and one-quarter, of the loads - 
by 32 and 33. tt may, | bec concluded that the actual col- 


> 
* 
: 
td f 
— 
[a j the total collapse load at the tip is, from Eq. 31 for 6 = 27, ere: oe 
— 
of the collapse load and may be compared with the solution (Eq. 27) of 
the came nro m ucing plastic theory. hen 2 er mit of the collansed _ 
4 
> 


ets could be approximately 10% or 


s&s In the preceding analysis concentrated loads were assumed to be distributed 
cover a small area to limit local shearing stresses and ensure the full limit - 
moment. The local ultimate bearing capacity under the load is governed wpeall 


i sistance to local circumferential tension failure, whichever is smaller “ail 
mm _ Surface bearing or penetration resistance can be estimated from a 


area 


tent and great h is ao) 


in which f is of the in lower limit of 
Eq. 34 applies to a perfectly smooth contact surface and the upper limit to a 


fora ‘concentrated load Pp is, 
= 0.33 (for a/h < 0. > 


If the contact radius is greater than about one- -fifth of the thickness of the — 
’ - bearing material or sheet, the surface penetration resistance decreases with 9 
smaller ratios of h/a until the lower limit of punching resistance is obtained 
a at about a = h. Punching resistance can be estimated from the vertical —— ff 
ponent of the total tensile strength of the sheet on a conical ———ee 


section inclined at 45° to the vertical and having an suite radius of a + h/2- 


7 the minimum co contact in order: for a load P to > avoid punching failure 


However, because this minimum contact radius the tensile and ‘shearing 


: strength of the e ice s sheet has been reached under the load, the corresponding | 
local limit moment is zero. With distance from the load the shearing : stresses” 
- decreases rapidly and the limit moment increases accordingly. Analysis has 
shown that for a concentrated load with a radius a =h the collapse load is 
Sow 105 less than that determined on the assumption of negligible shearing ~ 
stresses, while for a = 2 h the difference is very small. The contact radius of 
the load should, therefore, at least be equal to the thickness of the ice sheet. — 


Since: both the theoretical collapse” load and the local bearing 


lapse load of large radially-cracked she 
20% smaller than indicated by plastic theo 
= 
— 
— 
— 
— 
i 
A 
— 
— 
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should also to the thickness of the ice sheet ote 


APPLICATION OF. OF ANALYSIS AND COMPARISON WITH OBERVATIONS: 
ii order to apply the present : analysis to an estimate of the bearing capacity 
of floating ice sheets and compare this with the results of observations, it is » 
“necessary to know the yield strength of ice in tension and the moment- rotation — a : 
"capacity of ice in in bending. The most important physical characteristics 
_ fecting the strength of ice are its structure, density, and salinity (11). Ice is 
~ “an anisotropic material having needle- like crystals with a. vertical axis in 
_ floating sheets when the horizontal strength is about 25% less than in the ver- 
P tical direction. Stratification and a greater air porosity reduce the strength 
ol roughly in direct proportion to the resulting change in density. | Similarly, a -s 
_ greater salt content and higher temperature reduce the strength. ‘Thus, sea- 
ay water ice with about 3% of salt has approximately one- -half of the strength of . 
river or lake (fresh-water) ice at the same temperature, and, at approximate. 
ly 10°F, the strength of both types is about twice that at their melting points. 
‘The strength also decreases considerably | - with slower rates of loading and a 
- larger size of specimens. While the static strength with a few minutes to fail- 
: ure is less than the dynamic strength, the creep ‘strength deduced from labo 


ratory tests and the movement of (fresh- water) | glaciers was found to be less 
than one-quarter of the static strength and amounted to a yield strength inten- 

_ sion of about 15 psi to 45 psi (12). Moreover, the strength of about 6-in.- 

_ thickspecimens is only about two-thirds of that of 2-in. .-thick specimens (11), 

and even smaller strengths are likely to be given by very, 'y large specimens. > 
This is shown by cantilever bending tests on lake ice (13) and sea ice (14) in 
the field, when the flexural strength of 4 ft to 6 ft deep beams was only one- 


F rs] third to one-half of that obtained on small beams of 2 in. to 3 in, depth, — 


= At very slow rates and long duration of stressing the strain capacity of ice 
is large, as shown by creep tests and observations on the flow of glaciers 
(a2). It may, therefore, | be assumed that in this case the moment- -rotation 
a capacity of ice in bending is adequate to enable the limit moment My (using , 
“the plastic section modulus and creep strength) to be mobilized at the critical © 
ary Hence, the collapse load of floating ice sheets under long- time load- 
ing | can probably | be estimated as proposed for ideally-plastic plates. , How- 
ever, at normal rates of loading (a few minutes to failure), ice is fairly brittle 
_ and its moment-rotation capacity after yielding is limited. In the absence of 
experimental knowledge of this li limit it may be assumed that the most highly 7 
_ stressed or first critical section 1 of an ice sheet begins to crack in tension 
when the yield moment, My (using the elastic section modulus and on 
- strength) has been reached at this section. On further loading of the sheet the a 
moment capacity at this section is reduced, and the bending moment is trans- 
ferred to the second critical section until the latter is subjected to the yield 
moment and begins to crack. Again, the bending moment is transferred to the 
third critical section until the yield moment is reached | here, resulting 
~ eracking an and complete collapse under the ultimate load. ‘When there is only | - 
one critical section, as for a concentrated load at the tip of a wedge- shaped 
_ sheet with an apex angle of 90° or less and for a strip load near the edge of a 
- large sheet, it follows that the collapse load can be estimated as proposed for — 
ideally- plastic plates using the yield moment. For sheets with — angles 
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- capacity leads to premetere cracking under the yield moment at the first cri- 

- tical section (under the load). On further loading these cracks rapidly spread : 
to the second critical (radial) sections and subsequently atensile crack at the 
third critical (circumferential) section causes complete collapse under a 


concentrated load, while tensile cracks at the second critical (longitudinal) ‘ 
1 


sections produce collapse under astriploadg 
_ Because the first critical section is of small lateral extent (of an order of 
_ ‘magnitude as the ice thickness), the corresponding yield moment will _ not gov 
2. the collapse load under short-time concentrated loading. A safe lower es- 
timate of the collapse load under this condition can, therefore, be obtained by fi 
- taking that load at which the second critical (radial or longitudinal) sections | 
are subjected to the yield moment. Any third critical (circumferential) sec- 
tion is then still under an elastic radial bending moment, which is small com- 
ame with the yield moment and could well be neglected. By assuming that 
the ice sheet is cracked at the latter section (with radius or half-width R=b _ 
=2 mw approximately), the solutions for plates of medium size in conjunction — 
: Aime o the yield moment can be used to estimate the collapse load under short- _ 
time concentrated or strip loads c on n large ice sheets. Alternatively, a safe 
loads can be made a as previously indicated for wedges. it is 
that the first and second critical sections are cracked and unable to resist any ( 
_ bending moment, but the yield moment obtains at the third critical section. It _ 
is of interest to note that both methods give, sensibly, the same — load 
_ within the practical range of the contact radius of concentrated loads (0<a/ L 


4 0.5) (Fig. 4(a)). Moreover, this theoretical collapse 1 loadis about directly 
proportional to the ‘apex or wedge angle of the | sheet. 


_ It may sometimes be desirable to . have simple conservative formulas for 
_ estimating the collapse load of floating ice sheets under various short- a 


ora zh) .. 


is the apex aad of the sheet, “ear the stiffness radius o of the sheet (Eq. f 
: : and M,, is the moment of resistance of the section per unit length. Because the 
7 = ratio a/L is usually small and can frequently be neglected and M, = — : 

ee in which fy ts the flexural strength of ice on ‘the e basis of an elastic ie section 


- Pos! fp 
a. From the same considerations, the collapse ne per unit length of a i strip load ¥ 


Central load = = fp! h 2/21 


in which x is given by Eq. 2. In these « design equations Pom may be. taken as as the a 
total of an individual unit as a vehicle or because 


in particular fo apse 10a 
— 
— 
| 
— 
(40) 
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“the usual wheel spacing within a unit has no significant effect on the collapse 7 
ir practice. a row or line of such units (load train), however, the 
j total collapse load is given either by the sum of the ‘collapse loads of the in- 
} dividual units or by the collapse load of an equivalent strip load, whichever ‘is 
=, ee... it is known that a number of full- -scale loading tests to failure have » 
been carried out on floating ice sheets (15), results of such tests together 
_ with roe corresponding properties of the ice have not yet (1960) been published. — 


i Surveys have been made of the collapse loads and the — — maximum 
working loads on floating fresh-water ice sheets used as roads, bridges, and _ 
i landings of pulp and paper companies in eastern Canada (16), (17), and the re- 4 
7 sults are shown in Fig. 8. It is found that the observed | collapse loads vary 
mainly from about P, = 100 h? lb to 500 h2 lb, while the customary safe loads _ 
of traffic with speeds limited to about 10 or 15 mph vary from about P = 70 

h2 Ib to 250 h2 lb, where h is the thickness of ice in inches. The ert 
collapse load for central concentrated loads on large ice sheets can be esti- 
_mated from Eq. 39, so that the flexural strength of the ice deduced from the 
observed collapse loads would vary from about fp = 100 psi to 500 psi. a 


thicknesses, of fresh-water ice (11). “Moreover, the average 
i of safety deduced from the working loads and these flexural strengths 7 q 
is thus about 2, which appears reasonable for slowly moving loads when the 


strength of ice would govern. On: the other hand, an elastic analysis: 
a these flexural strengths (Eq. 1 with a = h) gives yield or first cracking 

— loads of about Py = 50 h2 lb to 250 h2 lb. These estimated yield loads - 

equal to the working loads used in practice, therefore, the corresponding fac- 

tor of safety is unity. Hence, some local cracking on the bottom of large ice 


‘Sheets under slowly — loads does not appear to be detrimental, as such ; 


mentioned previously, the theoretical creep load from Eq. 4 varies 
about P, = 25 h2 lb to 70 lb. This estimate indicates that long-time 
_ storage loads or parked vehicles should be restricted to about one-third of 
7 a safe slowly moving loads, since a factor of safety of unity could undoubtedly | 
be used in such extreme cases using the creep strength. On the other hand, 
_ vehicles moving at about 20 mph or more on ice over deep water and more © 
on ice over shallow water may set up resonance deflections and col 
lapse of ice sheets, which have been considered elsewhere (18). (ree q 
A further comparison of theoretical collapse loads with pen (emer- 
gency) and normal working loads on the basis of average flexural strengths 
is possible from recommended minimum ice thicknesses for landing aircraft : 
on floating airfields (19). This information is shown in Fig. 8 for single 
_ slowly moving or briefly parked aircraft on freshwater or far from the edge 
& of the sheets | and in the absence of resonance and other deteriorating condi- 
tions at two extremes of ice temperature (14°F ‘and 31°F). Using again Eq. ; 
ea 39 for the theoretical collapse load, it is found that the emergency loads cor- _ 


respond | to a flexural stress of about 100 ‘psi to 200 psi « or r assuming } a factor 
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‘aa = 150 | seb to to 300 Opel. ‘This strength is within the range of Sent test re- 

_ sults on fresh- -water ice and it is of the same order of magnitude as deduced 

_ previously from pulp and paper operations. Further, the theoretical factor of 

_ safety deduced from the normal loads and these flexural strengths is about 3 

_ against collapse and about 1 5 against - cracking, both of which ‘seem reason- 
able. a similar of minimum thicknesses recommended for 


by Black (20 hile 


Theory f, =11 500 
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FIG. 9. RESULTS OF F LOADING SECTORIAL GLASS PLATES 


tors of safety as deduced previously. 
as In the absence of full-scale loading test results of edge and corner loads on — 
i floating ice | sheets , laboratory t tests with concentrated loads on glass plates of | 


various shapes on a rubber foundation (20) support the ‘present theoretical re- —_ 


= 


—— 
— 
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— if 
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‘and of edge- loaded plates” is about one- »-half of that of 
-center- loaded plates (Pig. 9). Moreover, the lines of constant deflection of 


with the estimated age: at collapse and indicated a central sector between two 
roughly triangular outer zones with an edge angle of about 45°. Collapse oc- Nee 4 

_ curred by rupture along similarly shaped yield-hinge sections. A further anal- 
ysis of these loading tests can be made onthe basis of subsequent flexure tests » 
on the glass used (conducted by the author on small beams cut from an noriginal 
os kindly supplied by the National Research Council of Canada), which gave 

an average flexural strength of fp = 11500 psi with a coefficient of variation of 
about 40%. ‘In the plate loading tests the ratio of the effective contact radius” « 

7 of the load (measured at mid-surface of the plate) to the characteristic length 

or stiffness radius of the plate was about a/L = 0.03. Thus, the theoretical — 

_ collapse load of a (medium) center-loaded plate (Fig. a(a)) is about P,/h2 
= 15000 psi, which compares well with the observed average value of 15240 psi. 

7 Similar agreement is found between the estimated and observed average col- 
lapse loads of edge-loaded and corner-loaded plates as well as for tip-loaded 
sectorial plates with a smaller apex angle, as shown in Fig. 9. The scatter of 
the individual loading test results from the average values can be explained by 


_ the relatively large variation in the strength of the glass, as indicated by the 


If floating ice sheets are represented by plates on an elastic foundation, ythe 
yield or first cracking loads and maximum deflections ca can readily be deter-— 
. “mined from elastic theory. A method has been developed to estimate the cor- - 

__- responding ultimate or collapse loads from plastic theory and solutions have 
been obtained for various loading and boundary conditions. Simplified 


formulas for use in practice are supported by the results of field observations 
and some small- scale tests | 


uthor is indebted to D.C. Drucker of Brown University. and L Ww. Gold 


LOAD ON INFINITE PLATE 
Circular Load.—Aconcentrated load P is distributed uniformly o overa 
circular area of radius a on the surface of a th thin, rigid, ideally-plastic, and 

infinite plate of ‘uniform thickness h; the plate can resist a full plastic or limit 

moment, per unit length and rests ona discontinuous - elastic foundation 
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> 
ae 
3 
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E am previously. If membrane action and shearing deformation of the plate 
2 are neglected, the collapse load, Po» of the plate can be determined as follows. 
Let r and 6 be cylindrical coordinates, My and M; be the radial andcircum- — 
ferential bending moments, respectively, per unit length and Q be the vertical 
shearing force per length of arc with the directions as shown in Fig. 10. Then 


for the present case of rotational symmetry, the of are 
(5): 


in which p is the pressure on element of plate with downward forces and mo 
ments producing tension onthe bottom of the plate being taken as positive (Fig. 


Using the Tresca yield condition for the plate material (Fig. 11), it can be 7 
shown (5) that the radial bending moment My varies from Mp at the center of — 
the load P(r = 0) to - Mo at the circumferential yield hinge circle (r - b) and — a 
that with M, = 0 at an intermediate radius r=c, ‘say, the the following 1 relations a 


— 
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10,-ELEMENT OF PLATE UNDER LOA 
t 
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— 
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ICE 


M, 


and M; = My (Tresca regime B), 
3 O> -M, and My =M, + Mo (Tresca regime BC), 
at r- = - Mjand M = (Tresca regime C), 
Moreover, the rate of the deflection Ww tat plate has to satisfy flow 


: rule associated with the yield ‘condition, and in the wenn? case the the principal 
rates of curvature are given by the = 


which C;, Co, Cg and Cy, are constant of in the 
boundary a and equilibrium conditions. 


wig +10) b/e)] (for 0< < <a. = 


(in + In b/e) (for ¢ 


in which wa is the deflection “ie (at r= 0). 
ome radial deflection is thus in the form of an inverted obtuse ‘right. cone 
for 0< r< _c with the apex at the load. The conical shape merges at r= ey 
into a logarithmic shape which obtains 


The overall equilibrium condition requires that the load P is balanced by the 
- total reaction pressure onthe bottom of the plate, which is given by k times the 
volume of the deflection dish, where & = foundation modulus. Because this re- 
. 


the is 


w 


Putting Eq. 45 into 47 the gives 


nkw (b2 - 62/3) 


Eq. 49 into 41 and integrating with the boundary ¢ conditions 


— 
(44> 
ay 
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ed 
October, 


2 


4 


c2/3) 


2 : 


2 
d. 
oz a Substituting Eq. 50 into Eq. 42 and integrating with the boundary conditions 
(Eq. 43) continuity requirement of Mr atr=a ayields 
Kat 


rape 


In b/r) - c%(1+In - b? In e< r< 51c) 


je) b 
which is is the. ‘equation | for Wt 
_ tuted into Eq. 51b, and using the boundary aun of Mr = Oat r=c,a solu- 7 


d for various ratios of a/b. It is found that 3 . 


- the assumption of c > a holds for a/b < 0.2, but the above equations can readily — 

be modified for c <a, when a/b > 0.2. These results, which satisfy t both | static © i, 

and kinematic requirements, are shown in Fig. 3. 

} In the limit haa a i load when a= 9, oa 52 shows that c = 0 and the col- oe” 


The corresponding deflection of the plate outside ‘the loaded zone itself 


3 = by an upward reaction pressure on the bottom of the plate has not increased 

the collapse load. It may, therefore, safely be assumed that Eq. 53 is independ- _ 

ent of of the -sizeand — of the plate and that for a plate mater material 


a 
2° & 
— 
— — 
= 
“4 
= a 
= 
| deflection under the load is large at collapse while the overall equilibrium § [im 
= condition (Eq. 48) holds. This conclusion is borne out by the observation that _ q ree 


load would be 


in accordance with the result for a Same sas on a fixed circular plate (7). sit 4 
As the radius of the contact area of the load increases, the collapse load 
rapidly increases (Fig. 3). Within the practical ‘Tange of the contact a 


= 3.8.4 ‘ (55) 
‘contact areas the collapse increases more rapidly and ap- 


proaches infinity for a/b = = 0.8, approximately, when the rate at which the load 

_ does work on the deforming plate approaches zero. 2 Oe 
- ® Strip Load.—A load P per unit lengthis uniformly distributed over a narrow 
strip of width 2 a andgreat length resting on the surface of an infinite ideally- _ 

7 plastic plate, as before. An elemental strip of unit width along a transverse — 

section of the plate can then be considered as a beam of unit width, uniform — 
depth, and great length ona discontinuous - elastic f foundation (Fig. 1). Poe 
The deflection wof the beam satisfying the flow rule is given by Eq. 44a or, 


to the deflection under load (at r= 0). ‘reaction 
_ the bottom of the beam is given by Eq. re ee 


tis 


(1 /b + + P/2a ( for 0< ba 


p dr 


, after substituting Eq. , 58 and integrating with the boundary condition a M, — 


-1/3b) /b - 1/2a)] (for0< r<a). 


N 2(1 - fora< r< 

- Eq. 60b, gives the collapse load per unit 


— 8 | obeying the Mises yield cond 2 circumscribed 
— (56) 
a 
— 
— 
— 
— 
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Circular bicaailee —When the size of the e plate is less than that ‘aia to de- 
rm velop the radial limit moment My = - Mo at r= b, the collapse load will be 
_ less than that of an infinite plate, except in the case of a point load as men- : 
tioned before. The analysis previously given can 1 readily be extended to ob-| — 
tain a solution for a central concentrated ona circular plate of medium 


My = - Mp > - Moat r = —K 
this way a 


creases roughly in direct proportion to ‘the magnitude of Mj, above the lower 
limit for > Mb = O and that for a point load the same value (Eq. 53) is obtained © 


— bd 


_ For small circular plates when the deflection may be taken as uniform, the 


grating with the boundary condition of rQ=Oatr= 


(1/b2 - 1/a2) ) (tor r<a). 


2 ra< <r<b). 


Putt ing M 


(for 0.05 < 


: 
» 
ee which holds for a/b 0.2 if Mp/M,) = 1, increasing to a/b S 0.5 if M)/ M, = 0). ; “a 
The results, which satisfy static and kinematic requirements, are shown 
— 
PYESSure Delow the plate under a Central concentrated load 1s given. 
p Pp 
... (64a) 
the continuity requirement of My atr=a furnishes 
= Mg + (1/2 (for OX ....... (65a) £7 
2a/r - 3) (fora< r<b) ...... (65b) 
5) Eq. 66 satisfies 
static but no long kinematic requirements. The deformation of the plate will 
a g then not be uniform under the load but be similar to that of a medium plate, so : oa a 


533 for a point load. The corresponding re results are shown in. 3 and indicate 

a for greater contact areas the collapse load increases rapidly andapproaches 7 

infinity for a/b - -1, as would be fora loaded 


in which q is the via china of the plate. The half-width of the plate is then 


For a Scat load on a medium-sized rectangular plate (deflecting 
in a ‘wayto an infinite with an overhang of length t < t o» the mo- 


h 


by substituting ‘Eq. 67. Et to last Eq. of at r= b, 


or or 0 afo <2/3 and 0 <t<te) 


(1/o - (for 0< r<a) (71a) 
which b is the half- -width of plate. Substituting ‘Eq. into (59 and 
grating with the boundary and continuity conditions, yields 


a) (fora<r< b). 
Putting M, = =0 = =b into Eq. 12b, load per unit length 


Th The various solutions sien satiny ‘static kinematic r 
and are shown in Fig. 
2 EDGE AND CORNER LOADS ON SEMI-INFINITE PLATES, ponent sa 


Circular Load. —Let a concentrated load P be distributed over a 


— 
3; 
= 
ra strip load resting on a large plate the limit moment 
of My = - Mo can be developed at the section of r= b only, if the cantilever 
_overhang of the plate beyond r=bhasalengthoft,suchthat 
— 
— 
Foracentral strip load on small plates when the deflection may be taken 
4 
| 


edge of the plate with the moments My and Mj and the shearing 
force Qas before (Fig. 5). 
"are are 
modified to the extent that the bending moments and shearing forces are — 


along the free edge (@ = 0). Since the problem is symmetrical about @ = 90° cas 
omy one-half of the plate about the axis of symmetry wil). be considered sub-_ = 


‘Thus, setting M; = 0 along the free edge gives a conical deflection shape at 4 
° the edge (@ = 0) according to the Tresca flow rule (see above). “y As the angle 
_ @ increases, Mt increases until Mt = Mo (Tresca regime AB for 0 <r <c) 7 
and Mt = M, + Mo (Tresca regime BC for c < r< b) at 6 = = 6, say, as well 

as for larger angles. Accordingly, for 6 > 6, the combined conical and log- 

rithmic | deflection shape similar to that of an “infinite plate obtains, which ap- = a 

"proaches the purely deflection shape the axis of symmetry — 


4 but an approximate analysis can be made by using the flow field shown in Fig. y 
‘This assumes that the stressand deformation conditions in the central sec- 
_ tor (89 < 6 “ 90°) are identical to those of a similarly loaded infinite plate, 
Mle in the outer triangular zones (0 < a< 6 ) the deflection and the moment > 
M; vary monotonically with r and 6 as ane been indicated. For the whole cen- 
. _ tral sector an approximate solution can, therefore, be obtained from the pre- 
: = equations for an infinite plate by using an angle of 7 -2 6 © in place of 
2 7 for the overall equilibrium of the concentrated load in Eqs. 47 and 48. The | 
collapse load of the central sector from Eqs. 53 and 54 is then found to be di- 


in particular the critical circumferential section atr equals b, 
‘radial bending moment Mro equals - Mo. Further along any radial section of a 
the central sector (@9< 6 < 90°) it was found from integrating Eq. 51 and 
using Eq. 43, that -~ circumferential bending moment My equals Mo, very — 
: nearly. Hence, in particular Mto = Mo, along the junction between the central — 
and outer zones (@ = 6,4), where the radial and circumferential yield hinge 
= lines meet at a right angle | ( Fig. 5). Accordingly, from Eq. 74, 69 = 45°, which | 
‘ 4 indicates that the outer zones deform as _wedge- shaped plates with : a width of 
bat r= b, where b is the radius of circumferential yield hinge arc of central 
ss Phe wedges are subjected to a twisting moment since they remain st1 
no along the free edge while their deformation along the junction with the central "ia, 
a sector follows the deflected shape of the radius of this sector. Since the twist- 
ing moment reduces the collapse load of the outer zones below the value for hs 
bending.alone, anu upper limit of the collapse’ load will, in the first instance, ybe 
B _ obtained by ignoring twisting. A solution of this problem is then identical with 
of a an infinite ite wedge or without 


&g 
(| a 
— 
a 
= 
q 
4 
— 
— 
.—l 
= 
— 
fi 
— 


bending a and with a load at the tip, the load distributed o over 
a small og (or triangle) with its apex at the tip and having a radius (or | 

a (Fig. 6) . The collapse load of such a plate can be determined as : 


- Consider an 1 infinite wedge of variable half-width x ata distance r, , measured 


xer r tan (0/2) .... (75 


a in which @ is the apex angle (Fig. 6 6). 


i 


n which wo » refers to the deflection under the tip (r = 0) and b is the distance y 


of yield hing line from tip. y 
= Since the reaction pressure on the bottom of the e wedge is given by E Ea. “6, 


Hence, Eq. 46 6 becomes, after substituting Eqs. 76 and 


tan (0 /2) tan (6/2) 


05 < a/b < 3/4). 
a 


— 
— — 
(76) 
— a 
= 
bending moment per unit width | 
_or, after substituting Eq. 79 and 
Putting M, = - Mo atr 
— ae: 
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for a sectorial < 90°). These results both static and kin- 


wt. It is of interest to note that similar computations using Eq. 45 for the 4 val 
7 5 flected shape of the wedge have shown that the collapse load of the wedge _ 
me Closely given by Eq. 82, which, therefore, appears to be sensibly unaffected 
4 ¥ _ by the deflected shape. On the other hand, trial computations of the plastic 7 
: _ State within the contact sector (or triangle), where the bending moment is ae 
positive 1 under the load and conditions are similar to those of an infinite plate, 
have shown that the collapse load, estimated from Eqs. 82 and 83 is excessive , 
: for small contact areas with about a/b < 0.05. In that case the collapse load 
decreases rapidly with smaller ratios of a/b to one-half of the above loads at i : 
* 0. A similar decrease applies also to a plane wedge as a point load is ap- 
; proached, because the reaction pressure on the bottom of the plate below the | 
- load is then negligible, as shown above infinite plates. The results are shown 


a, a a | semi- infinite plate is obtained by combining the ‘collapse load of the cen- 
tral sector with a center angle of 90°(P,)=17 M,/2 from Eq. 53) with the col- 
: lapse lo load of both outer zones with a c¢ combined wedge angle of 90° (Po = 2 Mo, 


ie that is one-half of Eq. 82). Hence, the total collapse load of the semi- infinite ; 


= 


plates (90° < @ < 360°) is, , from 55and82, 


y= [1.65 (8 a/b) + 0.05 <a/b<0.2). 
from which the collapse load of an edge- -loaded plate is obtained for 6 = 
| with the results shown in Fig. 
: Se Comparing the results of Eqs. 84 and 85 for different ratios of a/b with the _ 
corresponding values for an infinite plate (Eqs. 53 and 55), the collapse load 
of a | semi- Somes plate is found to be about one- >-half of that of an na plate. — 
On the 
other hand, Eq. 83 shows that, on the assumption of zero transverse bending, 
collapse load of a sectorial plate increases indirect proportion to the 
angle of the plate, so that the collapse load of an edge-! -loaded semi-infinite 
plate would then be one-half of that of an infinite plate. age 
_ Strip Load.—If a strip load of width 2 a is applied adjacent to the free edge 
of a semi-infinite plate, the deflection is given by Eq. 56 and the reaction pres- 
ai sure under the plate is similar to that under one-half of an infinite plate. The — 


- maximum transverse re moment occurs at r= b and is, for the he cantilever — = 


of agent b, ‘given by 


@ < 90°) the collapseload oe 
for a wedge (0 < 0° 
ak 
> 
— 
ae 
— 
— 
| 
— 
a7 and setting M) = M>, the collapse load 
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= - infinite plate, but assuming that at collapse the stress and dihertention condi- 4 
7 a tions in the central sector are identical to that of a similarly loaded infinite 
aos : plate, while the moments and 1 deformation inthe outer zones vary monotonically © 


from the junction with the central sector (6 = = 9 as to the free edge. The cor-— 
responding flow fieldis shown in Fig. 5 for one- the plate axis 
of symmetry, as considered subsequently. 
Measuring from an axis parallel to the free edge and 
- load, the collapse loadof the central sector (6, < 6 < 90°) can then be deter- — 
mined from Eqs. +53 and 54 in direct proportion to the center angle o) 
as before. In the outer quadrilateral zones (0 <6< 6 o) the requirement of 
zero resultant bending moment along the free edge is, approximately, satisfied 
by the following expression containing the terminal bending moments along the = 
yield sections of the outer zones (with “right angles at junction with ‘central — 
and at free edge on axis which is a yield- as 
= 
_ which b is the radius of circumferential yield hinge arc of central | py ie 
_ Ignoring twisting moments in the outer zones and applying to them the prin- . 
ciple of virtual work, an upper limit of the collapse load is obtained inthe first 
_ instance. In this way it can be shown that the the collapse load of each — 


-4 oer b)(b - 


= _ yield hinge line (r= m ta DAx = first moment of outer zone about outermo 
yield hinge line, and the ratio of these terms represents of the 


7 ‘The total collapse load of the semi- infinite plate. is obtained by combining — 
" the collapse load of the central sector (using Eqs. 53 or 55) with the collapse _ 
load of each 1 outer zone (Eq. 90). Hence, for a point load (a= = 0 and negligible — 
reaction pressure below plate, « as for an infinite plate) the total collapse load _ 


icfie. c requirements, is shown 
AT SHORT DISTANCE FROM EDGE 
Circular Load.—Let a concentrated load P be unifor 
Circular Loa -—Let a <pbf the free 
circular area of radius a ata centroidal distance of e 
edge of an ideally-plastic semi-infinite plate, as before ( Fig 
— 
— 
a 
st. 
the 
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[ose + 68/0) + 


1-4a/ 3b 


0,08 < b< 0.2and0<e 

- very “nearly. An allowance for twisting moments would reduce these values 


| 7 somewhat, as woulda statically admissable lower bound solution in the outer 


Strip Load.— if a strip load of Width 2 a is at | a centroidal distance 


load (at r = e) becomes equal to the limit moment, Mo. The upper limit of the 
collapse load is given by an equation analogous to Eq. when the reaction 
- pressure ‘below the plate between the sections at tog 0 and r= =e becomes | tri- 


angular. In this way & can can ee shown that the collapse load per unit length for 


Po (tor 0.4 <e/o<o. 6) . (98b) 


(fo rb< b< eS b+ to) 
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E STRESSES 


The computation of membrane stresses and of the corresponding deforma- 


Et _ tions for a hyperboloid shell of revolution under the most general loading Se a 
ra reduced to the evaluation of simple integrals. For the determination co 
effects, diagrams capable of simplifying computations are presented and an 


- 7 


Reinforced concrete shells having as middle surface a rotation hyperboloid 4 
_ currently used in the construction of cooling towers. Wind action consti- * ; 
tutes the severest loadingfor such structures. If the distribution of wind pres- 
sure along any parallel circle is assumed to be proportional to the cosine of r 
; = angle that the radial direction makes with the meridian p plane containing the | 
wind direction, the of the stresses can be performed according to 
known methods. 2 However, it is recognized’ that the wind pressure is 
distributed | quite from the ‘mentioned law and that much 


SR 
7 -Note.—Discussion open until March 1, 1961. To extend the closing date one month, a 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of me 
ss American Society of Civil Engineers, Vol. 86, No. EM 5, October, 1960. 
ae Professor of Applied Elasticity, Polytechnic of Torino, Italy. «ee ai 
_ 2 Theory of Plates and Shells, ” by S. Timoshenko, McGraw-Hill Book Co., Inc. N 
“Sulla teoria elastica della parete sottile di rivoluzione e l’applicazione agli in- 
volucri ” by P. Cicala, Corso di Costruzioni in Cemento 
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yr in the stress analysis of the shells having this middle surface. “This 
_ was pointed out by Fliigge* who indicated a homogeneous solution; that is, a 
stress system produced by a special edge loading. Here the problem is dealt 
with in general terms. The solution for membrane stresses and 
a erry under any arbitrary loading is obtained. These results can be = 
rived at using the theory of characteristics for the resultant equation. This is 

a second order partial differential equation of the hyperbolic type. However, 7 

an elementary deduction has as leading, more immediately, to 


the final results. | 


Analysis of Stress. ~The of the ‘membrane stresses 
or a surface of revolution are written? as follows: 


Ng sin | 


Ng 9 9 are the ne normal and forces per unit arch of the 
are the normaland tangential forces per unit arch 
Tr of the meridian; , Y, and Z are the components of the external force per unit | j 
surface on the ee to the parallel, the tangent to the meridian and the nor- _ i. 
mal to the surface, respectively; r is the radius of the parallel; Ty the = 
of curvature of the meridian; 6 is the angle defining the position of the merid- 
ian plane; and 9 is the angle made by the normal to the surface and the axis 
of rotation, . Along this axis the coordinate z is measured, the elemental 


arch the meridian the equivalent hold 


cos 


sing gL to denote a constant length, we notation 
al 


stresses arise than computed from the cosine distribution giving the same re- 
fst ise tha ted from th ine distribution giving th ——— 
j sultant thrust. A more exact computation is therefore inorder. The peculiar 
= 
“Statik und Dynamik der Schalen,” by W. Fligge, Springer, Berlin, 1957. 


stituting Eqs. 5 and 6 into Eqs. 1 and 2 and intinie account of Eq. 4 leads to fi 
the equations defining the functions: 


(Ta 


or rdinates 1 r, z may be expressed in the 


Here Ris throat radius of is the Je by th as- 
ymptote s of the hyperbola with axis of w° ranges 


_Lower redge 


- -90° to +90° while the point running the hyperbola goes from z= - @ to. 
+@ af varies: from 90° + ae to ‘90° - by means of 


cos’ 


— — 
— 
— 
— inction, according to Eq. 10we get 
a generic function, according ‘neil 


_ The surface of the hyperboloid may be mapped into the rectangle ABCD in 
the plane of coordinates 6,W (Fig. 2). To every point of the hyperboloid there 
7 a corresponding point within this rectangle. Parallel circles and meridian 


‘ lines are transformed into horizontal and vertical lines, respectively. Gener-— 


a 
one > such line as PR, the ‘slope dw/dé = 1 have 


the derivative b taken along sueh a Integrating from P P to R yields” be 


is 


P and T are points on a free edge, there is Hy = 


"stress resultants at | resultants at Rare written | 


@ 
\dw-} (Fy - Fp) Jaw. 


« 


using Eqs. 10 to 12, if L = R2/ 
— 
ia — 
a) 
.—l 
+ a) = Fy + 
— 
0. Thus the 
a 
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If the is symmetrical with respect to the plane | 0, that 


then, for any couple of integration as oil 


‘Parther simplifications arise if the load distribution | is also | 


with respect to the w= 0. That is, u 


Then the relationship (See Fig. 2) 


ported, 

is obtained. In conclusion, in the of double ‘symmetry, the ‘integrals: in 
Eqs. 20 and 21 need not be computed on both systems of lines dw/dé and 
suffices to integrate along one system, on one side line 


)= 


in in which the arguments of the integrand are We “Hence, f 


=- w'is a free edge, the stress resultants are expressed by on 
6)= Tw, 6 )+ -8)- S(w" w' +w-6)- w' (21) 


path, such as as HI (Fig. 2), intersects the contour @° =90°, it 


to be prosecuted as indicated by line 


: 4 the hyperboloid is subjected to normal pressures p=Z, with X= Y=O,and_ 
if p does not depend onw, then, by Eqs. 14, replacing ®p/2@ by dp/dw ba ac- 
cordance with Eq. 15, ‘Eq. 26 be written 


Jo sin® )/cos” 


— 
— 
— 
— 
— 
— 
i 4 
— 


e corresponding to the argument 6+ 


ure d 

: ch m isa positive integer, the stress resuitants 2 are ‘given by 

= Hy cos m = Hy” sin m 

in which H,* * are fi e functions of m, a,w "and w. For m = 


tanw" 


and Ww, may be deduced from Eqs. 21, 28, 30, 31,and 32, or, more directly, 
iat Analysis of Deformation. —The deformation of the membrane shell of revo-_ 
lution is expressed by the relationships 


in which €,, €,, and are the strain components of the middle surface (cir- 

‘2 _ cumferent 1, meridional and shear strain, respectively); E is the Young’s 
modulus, the Poisson's ratio; v, u, are the components of displacement 
pe - tangent to the ‘meridian, tangent to the parallel circle and normal to Ge aa 


5 
d f 1 
respectively; and t is thickness © fe) shel = 


stands 
flos 
— 


. fed ~ the particular case of the hyperboloid shell, the sa same Cintra 
+ — 


€. 
2 w sin” o 
Hence the unknown functions. , Hy “can “be computed byt the. same pro- 
cedure as suggested for Hy 
Boundary Conditions.—The situation for arbitrary enc end-conditions on the two 
‘parallel edge circles has been dealt with in the mentioned note where | 
of strain represented by sinus functions of 6 have beenconsidered. In the gen- 


case the solution is obtained by combining the following: 


ot 1 da particular solution of Eq. 13 for stresses eat of Bas. 42 for displace-_ 


ments taking account of the applied loads; 


2) +two homogeneous ae solutions, that is two independent solutions 


3) two homogeneous inextensional solutions, that is, two solutions obtained 
‘The solutions of of type eta * are well known | in the case of wien os symmetry. 2 They 
a can be determined? in the general case for a shell of revolution by using a — 


power expansion of a thickness parameter. 


‘conditions along each edge. With upper edge free, the three components of the : 
stress resultants and the bending moment are to be equated to zero. Along the 
lower edge, the shell forming the cooling tower is supported by trussmembers. _ 

- The constraint condition thus afforded may be assumed to nullify the compo- 

7 “nents u, v of displacement as well as the transverse shear and the bending mo- 
ment. In this situation, if solution! is by prescribing that H, = Hg = 
ml along t the upper edge =-y' and that H,' = Hy’ = 0 along the other edge w = w" SS 
then no correction is needed, as previously shown for any component in = 
Fourier expansion of strain configuration. In other words, the 


“The industrial al Cooling 7 by K, ai and d M. | Brooke, Elsevier, Amster- 


— 
tm 
= 
= 

=. in which 
a" 
— 
= 

Se 
tm 
4 

— 
— 
— 
— existing information concerning the distribution o 
lete. Consequently, McKelvey and Brooke? 
bodies of revolution is not yet complete. 


1960 


“used the the _— as given by French regulations for long cylinders, to compute 


Bs y wind effects on conical cooling towers. If these same data are assumed to be 
valid for the hyperboloid, the results plotted on Figs. 3, 4, and 5 are obtained. 7 
i The abscissa @ is the angle made by the meridian plane with the direction of — 
7 wind velocity, 6 J = 0 being the upstream meridian. The reference pressure Po 
is the dynamic pressure (half air density times square velocity) . The data 
Wy referred to indicates that in the downstream sector (that is, 9°<- 90° or > 90°) 
the pressure is constant and equals -0.52 pp. Therefore, the computation is 
carried conveniently for pressure values increased by 0.52 Po, the effects of 
the additive constant being separately accounted for. _ The line a of Fig. 3 rep- 
_ resents the adopted pressure distribution. The values marked near the curve, 
iis a spacing of 5° or 10°, are taken from the aforementioned regulations or 
interpolated. They were used to compute J; and J9 according to Eqs. 31 and ~ 
= The integrations are ee by the Simpson formula, over intervals of - 
For a constant pressure, me 31 yields J, = 0. Therefore, the values of Jy 
aren not warp by the constant added to p. They are plotted against 6 in Fig. 
, 60°. Each curve is to be read according to the scale 
7 P marked at its pe The curves should enclose zero area, except for slight im- 
precision due to the integration process. 


the values of Jo = - 0. 26 W, , to | be read according to the scales 
at the end of each curve. 2 The mean ordinate should equal - -0.125 tan2 wW since | 
0.25 po is the mean ordinate of curve a in Fig. 
aa of the interval - -90° < 6° < 90° + w® the values for Jy and Jo are zero. 
Both series of curves, in Figs. 3 and 4, show more or less abrupt changes in 
slope for = 90°. + and for @° = 90°, in addition to their end-points. 
By using the charts of Figs. 3 and 4 for double interpolation, the stresses 7 
- at any point o of a hyperboloid shell c can be € computed | for any set of data, tine 
the value w® = 60° is not exceeded. This is usually the case in design practice 
_ of cooling towers. ‘The geometry of the surface is assigned by fixing two out 
the quantities a, anda, related by Eq. 9. The between 
# and w is established by means of Eq. | 8 and the value w" corresponding to the — 
_ free edge is computed. Fora given point of the shell | w= = Wy, 6 = 6, the values = 
for J; and Jo = Jo + 0.26 w are determined from keg charts of 3 and 
4 for the arguments w = = 61 and @ = =w W, + - 6, and 
6 = +W" + Oy. Then the corresponding J values are computed from Eq. 
id and suntieiis in Eqs. 27 and 28 to obtain H, and H2 at the point considered. 
Finally, the unit forces are computed from Eqs. 5and 6, by the shell 
7 an it is to be remarked that the angle a enters explicitly in Eq. 30. This mal makes — 
indicated charts valid for any value of a 
‘These diagrams have been deducedunder theassumption that p dependsonly 
o on 6. A generalization is readily feasible. If the pressure distribution is rep- 
byalaw suchas C(z) ple )in which C is a coefficient having a stepped 
& diagram when plotted versus z, then the value of J may be expressed by the 
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summation © in which AjJ is is the Aifference of the the values taken 
“the ends of the generic interval where C = Cj. 


particular example representative of usual cooling tower design has been 


considered. The values fixed are (a/R) = 2.4, w’ a 22°, w’' = 56°. The cor- 
related shape of the shell is represented in Fig. 1. The line a of Fig. 5 shows 
a the values of Hy computed along the lower edge of the shell according to the 

graphs of Figs. 3 and 4. The dotted line b indicates the values obtained from — 

= same data when wrnsiind into account only the Seat harmonic in the Fourier 


FIG. 5 —NORMAL STRESS DIAGRAMS AT LOW ER EDGE F 
__unconservative if it is carried by assuming the wind pressure to be propor- 
tional to cos: 6. this is the from which the resultant wind 


7 ae 7 “Etude aérodynamique d’une tour réfrigérante en forme d’hyperboloide de révolu- 
by G, Golubovié, » Int. ‘Ass. for Bridge an and Sti Eng. XVII Vol., 1957, P- 


— 
— 
q _ The same structure, in the pr by line a in Fig. 6. This 
en from the data obtained by — — 
7 a is taken from the F 
tm 


October, 1960 

‘ments on a hyperboloid model having a smooth ‘surface. In this case, 

a strong depression exists in the neighborhood of the meridian plane perpen- Be 

 diculartothe wind direction. On the contrary, measurements on models having © 
rough surface effected by Golubovié and others indicate diagrams akin to 
the one to which line a in Fig. 3 refers. The computations were carriedon the — 

_ basis of the dotted line of Fig. 6, that follows the original graph closely but a 

corresponds to a simple diagram for dp/dé , as shown by line b of Fig. 6. On 


this basis, the integrations were analytically. The values ‘for Hy 


6.—PRESSURE DISTRIBUTION ON AHY PERBOLOID 


7 _ the ‘the lower edge of the structure were commie. The most significative 


branches of the diagram obtainedare represented by line c of Fig. 5. Of cauie, 
a the pointed peaks, such as the sharp corner st C, are actually rounded off by 


neighborhood of the (positive or negative) maxima of the pressure diagram are 
not be be reduced greatly by the said Therefore, comparison of 


We 
+ 
— 
second case considered, — 


The humerical investigation showed that the maximal values of stresses are 


a mainly « "due to the component proportional to cos 2 6 inthe Fourier expansion 
of the pressure diagram. This fact suggests the use of Eq. 35 for an approxi- — 
mate evaluation of these maxima. According to actual design trends for cooling — 

towe towers, containing sin? a@can be Thus the modified 


is suggested. This formula furnishes the value of the maximal posi- 
tive and negative direct stresses. The unit Ap is the interval of values within 


which the pressure ranges. The remaining symbols are explained by Fig. . 


_ Application of Eq. 45 in the cases considered shows excellent agreement | 

with the results of the rather lengthy exact computations. In the case of pres- 

_ ~sure diagram a of Fig. 3, for which Ap = 1.8 po, Eq. 44 indicates the range of 

Bi for Hy denoted by Ag on Fig. 5, in | agreement with line a. For the pres- _ 

sure diagram of ‘Fig. 6, ging Ap = = 3 Po » the ne range of values for 


computation of the membrane solution to be reduced to simple integrations for | 
any distribution of loads. Further simplifications arise in the case of eae, 7 


| peeeonres | distributed according to laws that depend only on 6. The charts — 

- “presented greatly reduce the computational labor. They are based on a stand- 
ard diagram of pressure and allow one to take into account the influence of a 
the: ‘geometrical parameters of the shell. Moreover, an approximate formula 
of immediate use is suggested for the evaluation of maximal direct stresses. 
- Concerning the quantitative results of the present investigation onthe effects 
of wind action, it is to be remarked that the consideration of the actual distri- 
» bution of pressure leads to much higher stresses than the ‘simplified analysis 
based 0 ona distribution proportional to cos 6. Furthermore, considerable dif- 


For instance, | precise quantitative indication onthe influence: of surface ‘rough- 
-ness should be offered to the stress analyst. As long as the information eS. 
ual lacking the calculation cannot be carried out with sufficient reliability. The 

3 diagrams presented here are intended for demonstrative purposes, however, 
they inated ‘show the practical importance of the problems dealt with. 


diagrams obtained on aerodynamic models having a smooth surface ora rough 
one. A paradoxical result is reached on this point. The wind action on the 
ae) smooth shell gives a smaller resultant thrust but higher stresses than in the Be ay 
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BEHAVIOR OF STRUCTURES DURING EARTHQUAKES? 


- spectrum is computed at the time that gives the maximum value. - Because of 
+3 this, it is found that the average earthquake spectra indicate that stationarity 
has not been reached for the longer periods and smaller damping. 
2 Mr. Berg’s finding, that for elasto-plastic systems on the average the en- 
Ps ergy input does not exceed the maximum input for an equivalent, elastic, un- 
tai damped structure, is very interesting and helps to explain the observed a 
ior of structures during earthquakes. _ He is, of course, correct in stating that — 
a the maximum energy input should not be e computed from a low point on the un- 
a In his discussion, Mr. Derrick makes numerous edema of a rather gen- 4 
eral nature whose complete rebuttal would be beyond the scope of a closing dis 
cussion, following comments are offered in responseto the stated and 
plied points of difference set forth by him. 
The El Centro, 1940 ground motion is the strongest so far recorded. It is ae 
certainly not the strongest possible ground ‘motion. | That ground a has =f 
been estimated by the writer2 to be perhaps | twice as intense (63% g) without 
making any special allowance. for type of soil. Of course, in the special cases 
of structures located directly on top of a slipping fault, or some unusual geo- 
logical conditions, there is no telling what motions might be involved. 
‘The El Centro accelerometer was not located on sedimentary rock but was 
- locatedon very deep, sandy alluvium, and for this reason seismologists attrib- 
ute an amplifying effect to the site. For example, Gutenberg? compared seis- 


mograms obtained on the rock outcropping at the Pasadena Seismological Tabs ; 
oratory with those obtained at Long Beach and at El Centro. According to him, 
» °- the amplification factor at Long Beach was 2.5 and at El Centro was 8, or in > 
rs other words, , for or earthquakes | of the same > magnitude and | at the same ® distance, 


October, 1959, by George W. Housner, 
Asst, Prof., Applied Mech., Calif, Inst. of Tech., Pasadena, | 
: “Intensity of Ground Motion During Strong Earthquakes, * by G, W. Housner, Second — 


Tech. Report, Earthquake Research Lab., Calif. Inst. of Tech., August, 1952, 
“Effects of Ground Shaking in Earthquakes, ” by B. Trans, Amer. 
physical Union, ‘Vol. 37, No. 6, December, 1956. 


— — — 
—— ing and novel comparison between earthquake spectra-ratios andthe idealcases [im 
a single-impulse anda constant-energy input. The spread of these response 
computations can be attributed to differences in earthquake records, and also 
to the manner in which response spectra are computed. Inthe case of the con- > 
4 stant energy input with effect of damping being as V1/v, the ensemble average 
_ is computed at the end of a specified time that is sufficientlylong to permit the 
4 
2 
on 
— 


ox ground acceleration at El Centro would & in a the ratio of 8 to 2. 5, or 3. 2 


according to the maximum El Centro- 


oer have been some five or six times greater than at Long Beach. However, 


neither instrumental recordings nor damage observations indicated that the 


—— motion at El Centro was appreciably more intense than in the Long r 
Beach-Compton area. . The records obtained during the Long Beach earthquake © 
and those obtained during the El Centro shock are consistent with each aad 


if it is assumed that there was no significant effect of ground.2 


It should be noted that Gutenberg’s measurements were of extremely weak 


q 


a that were far below the threshold of human sensibility. His 
findings are not directly applicable to strong ground motions. For example, — 

.-. addition to the El Centro-Long Beach discrepancy noted above, the acceler-_ 
_ ograms obtained during the San Francisco earthquake of March 22, 1957 do not 


_ 4 agree with his findings. These accelerations were recorded at five different 


- locations that were on filled ground, on alluvium of various depths andon rock. 
; The deepest and softest ground did not amplify the motion over that recorded | a 
; 4, 5 While it is true that the local geology can affect the earthquake _ 
ground motion, it appears that it is not a simple relationship that depends only | 
_ upon the depth and softness of alluvium. More strong earthquake recordings | 
4 on different types of soils are required before the —_— of ground can be elu- — 


The average spec spectra v were obtained -accelerograms located on on different 
types of ground, however, despite this, the individual spectra were all similar 
7 in their general shapes and did not indicate any significant influence of ground. 
7 1940 as did those of Taft, 1952, which accelerograms were recorded on sedi- - 
“mentary rock several hundred miles distant from El Centro. 
.* _ The average spectra, when multiplied by the prescribed factors, do not re- 
- produce the actual spectra. The factors that were given in the paper, scale © 
the average spectra so that ‘the spectrum intensity corresponds to that of the 
<. average of the two. components of the actual spectra. Presumably, the point — 
that Mr. Derrick is making is that the ordinates of the actual spectrum curves 7 
have some statistical distribution about the average andif the range of this dis- 
tribution could be determined, it would be better than the average. There oa 
however, no reason to suppose that this would give significantly different re- 


— Specifically, Mr. Derrick points out that one of the spectrum curves for — Lo : 


the El Centro 1940 shock is higher than average in the low period region. This 


ean probably be attributed, at least partly, to the fact that the accelerometer Ea 


= within 10 miles of the surface trace of the fault and hence the high fre- 7 
- quency components of the ground motion were only lightly attenuated: It would iC 
, probably be correct to say that for points so close to a fault the average spec- | 
trum curves should be boosted in the low period end. — = other hand, =< 


4 


— 
case during strong ground motion. The 1933 Long Beachearthquakehada mag- 

= 
Bead 
> 
f 
4 
— 

— 

| q of March 22, 1957,” by D.E, Hudson and G, Ww. Housner, Bull. Seismo. Soc. Amer., Vol. ny a 
_ 5 “Response Spectrum Analysis of the March 22, 1957 San Francisco Ea: 
j. OD, E,Hudson and G, W,Housner, The San Francisco Earthquakes of March 
ial Report 57, Calif. Div. of Mines, 1959, 
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points at greater distances fromthe fault hes average spectrun 


_ perhaps, be depressed in the short periodend. If the average spectrum curves | aM 
=. are used for design it should, of course, be kept in mind that the spectra of fu- 
earthquakes will have statistical fluctuations about the 

oy Mr. Derrick states that, “the assumption that the distribution of mass does 

“not influence the amount of work performed by an earthquake is open to ques  « 

a tion.” The writer believes that the statements regarding this matter that esa 

_ made in the paper are correct so that Mr, Derrick’s question presumably a 


rises from a misinterpretation of the statistical implications of those state- 
ments. Mr. Derrick states that, “there is no positive evidence on which to base of 
_ either an analytical or an empirical solution involving plastic action or damp- _ 
ing it in modern type multistory buildings.” If by this statement he means that — 
we are not able to reduce the problem to | simple design rules and that further 
_ research is required, he would certainly find everyone in agreement with him. 
There is, of course, a great deal of positive evidence available on plastic be- 


_ havior of steel frames, on random excitations and on non-linear vibrations. i. 


‘Mr. Derrick makes the statement, ‘that the conservative designer has no al-- 
ternative but limitation of distortions to linear or almost linear levels; many 
designers are already doing this.” ” In comment | on | this, the writer would say | 
 thatno multi-story building has been designed to remain elastic when subjected 7 
to the strongest possible ground motion, as, for example, the most intense mo- 
_ tion produced by an earthquake of 8.6 magnitude. If the linear design is based a 
on less intense ground motion, as that at El Centro, 1940, then in the event of © 


_ the 8.6 magnitude shock the structure will go beyond the elastic limit. It il 
in the writer’s opinion, very important that designers recognize this fact. © atte 
e Mr. Blume has commented onthe reliability of the average spectrum curves. 
As this is a point that may be misunderstood, a short explanation will be given 
7 here. All of the many ground acceleration spectra so far computed are con 
sistent among each other, except two that will be discussed later. The four al 
ground motions thatwere recorded near the epicenters of relatively large earth- 


¢ ‘Taft records were obtained on sedimentary rock and the Olympia, Washington 7 
records" were obtained on terrain of quite different geology. Ground ‘motions 


. 4 ’ = are consistent with these, along the lines set forth inthe paper. The only spec- 
: a 7 tra that are not consistent with these are two cases where relatively strong — 


& ground motion was recorded on extremely soft, water-filled ground. These are 


7 the Seattle, Washington, 1949, and the Southern Pacific Building, San ‘Francisco, 


1957 records that are referenced in the paper. It may be noted, however, th 


_ in both of these cases , there was no evidence that the local geology had an am- | 

plifying effect on the ground motion. It may be stated that, other than such ex- i. 

-ceptional cases, the average spectra of future can be expected 
Although g geological conditions “may be expected ‘to have an influence on ‘the: 

ground motion, it is quite certain that large influences can be expected only if 
_ the geology is markedly different from the geology of those sites where most Z 

of our past earthquakes have been recorded. The precise shapes of individual _ 

spectra, will, of course, be subject to statistical fluctuations. The = 

_ intensities will depend upon the magnitude of the earthquake and the distance 

from the as well as affected such as and extent 
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 ~recorded at greater distances trom the epicenters Of large earthquakes and at 
dictancoc from the aniconterc of aemall have enoctra that 4 
4 
| 
at 
— 


1900 
faulting. It 1 must be in mind that ground motions are statis 


‘ earthquakes is that their average spectra can be expected to be like those of 
_ past shocks, and the statistical deviations about the average can also be ex- 


ne Mr. Blume a appears to have misunderstood the significance of the statements = 
made inthe paper regarding energy input and energy dissipation of over-stress- 
ea structures. The exact analysis of non-linear vibrations of multi-degree-of- __ 
freedom actual structures has not yet progressed to the point where its results 
; 4 can be ; applied in design, and, therefore, design based on inelastic behavior 
F must perforce make use of a more crude and approximate analysis. The crux © 
_ of any such method is the way of estimating the energy input that must be ab- _ 
; _ sorbed. Ih his 1956 paper6 the writer proposed the method set forth in the : 
c. present paper (Eqs. 20 and 21) which Mr. Blume also adopted in his paper. | It 
a = , of course, not obvious that this is a satisfactory method of computing the > 
energy input, for the input might actually depend upon the energy absorption of 


™ the system. Recently | the writer’s colleague, T. K. Caughey, completed a math- a 


ematical proof that for a random (white noise) excitation of a single- degree-_ 
— of- freedom structure with viscous damping, the proposed method (Eq. 21) is 
indeed correct as regards average values of energy input for different amounts 
4 of damping. This fact, together with Mr. Berg’s findings on inelastic action © 
= _ with earthquake excitations, makes it quite certain that the proposed method is 
= § reasonable to use solong as the increments of inelastic energy dissipation are 
4 small compared to the total vibrational energy. ‘The ability of the structure to | 
absorb energy, Mr. Blume’s so-called reserve energy, is of course a necessary 
_ piece of information in the application of the method. For a 
ie quantity can be computed, but for complex structures it is not yet kno 
_ how to evaluate it, for the various locations where yielding, and so on, will ; - 
velop depend strongly upon the precise properties of the structure as well as — 
on the prior history of yielding, 
_ In his Fig. 1 and Table 1 Mr. Blume compares the spectrum for the average 
of the two components with his curves that represent single components. ks: In 
Table 1 of his paper, the writer has prepentes the spectrum intensities - one 


how these differ from the averages. 


- - the paper are the statistical (ensemble) averages and they should not be con- | 
ou fused with t the statistical fluctuations of the ‘samples about the le average. Neither: 


= and ispectrum intensities that were presented inthe ‘paper. Any method of of spec- 


-. computation is, of course, subject to inaccuracies in recording, reading . 
and plotting, transcribing data, and so on. These errors will show upinthe | 
_ final results, so that in earthquake spectra computations there is no possibility F 
of obtaining 100% accuracy. In this case, the best check on the accuracy of a _ 


get of computations isto make a comparison witha completely independent set _ 


6 «Limit Design of Structures to Resist Earthquakes,” by G. W. Housner, Proceed- 
ngs, 1956 W World Conference on ce on Earthquake ‘Engrg., Earthquake Engrg. Research _Inst., 
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DISCUSSION 
culations. Sneha completely independent computation of spectra has been 
e by Mr. G. V. Berg | on the digital computer at the University of Michigan, 


pes 1 Arbor, Mich., and satisfactory agreement » was found between his results Daa 
and those of the writer. The differences were such as might be expected from — a ii 
ty independent handling of the raw data. The writer would also like to take this 4 
opportunity to clarify the regarding the accuracy of the spectra that 
= made by Mr. Blume.7 His paragraph that begins, “Another interesting 
' = was that for stiffness system C the 1940 and 1934 El Centro earthquakes © "3 q ys 
ne produced the same absolute maximum base shear for the low damping val- 
ues. ..,” states that there was poor agreement between the spectra and the a ” 
computed response of multi-story buildings, particularly for low damping. In- 
dependent computations now show that the discrepancy can be attributed to er- 
. 2 rors in the dynamic analyses presented in Mr. Blume’ S paper, and it is not the - ; 
result of inaccuracies in the spectrum curves. 
_ Mr. Blume and Mr. Derrick make statements and recommendations con- 
G _ cerning methods and values to be used in the design of structures. With some 
> of these the writer would agree, and with some he would ‘disagree. However, 
‘since the design | of structures is beyond the scope e of the paper, andsince a rig-— 
orous treatment of design would require more space than in aclos- 
ing discussion, these questions will not be treated 
a The writer agrees with Mr. Cohen and Mr. Weissman that the maximum = Ss 
ee earthquake is not the same as the maximum experienced earthquake ~ 
- and it was not intended that the or should imply that it is proper to assume Ae 
any engineering discussion CF. Richter’s “Seismic. Regionalization 
Maps” his method of constructing them should be pointed out. His maps are ay: : 
- combinations of seismic zone maps (showing where earthquakes occur) and = : 
geological maps (showing the type of soil). He makes the assumption that, _ 
other things being equal, the intensity of ground shaking will be proportional to *) . 
the depth and softness of alluvium. He makes no claim that his maps have an 
_ accuracy suitable for engineering purposes and, in fact, measurements of strong 


_ earthquake ground | motions show that lee: method of constructing the maps i is 


still some time off. _ These discussors also present an interesting description 


q Maps. As Mr, Cohen and Mr. Weissman state, engineering ‘risk zoning 


‘ the nuclear r explosion ground shock problem and | point out that, in a 


it is the same asthe earthquake problem. It appears to the writer that any dif- 

ferences between the two problems are probably such that the nuclear 

more difficult than the earthquake problem. For example, in the nuclear 

gan coupling between ground and structure is often important, whereas in| 
red the earthquake prion (in the United States) this appears tobe important only — 
exceptional cases. Similarly the influence of local “appears to be 
"much more pronounced in the case of the nuclear problem. OT x 
closing this discussion, the writer would like to thank Sennen. Berg, Der-— 
Set rick, Blume, Cohen and Weissman for their discussions. He would also like to 
: — that his paper was intended to give a seuniue (and necessarily in- 


account of the essential facts known about the behavior of structures 
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October, 1960 


on the actual design of of tens and building code 


E E as the ‘magnitude e of earthquake to design against, the e allowable stresses to of 
used, the appropriate factors of safety, the permissible degree of simplification © : 
eS and approximation, the empirical curves to be used, and so forth. In this re- or. : 
gard, the following two statements are, , perhaps, in order. First, it is not pos al 
3 ile to make adynamic analysis that will describe exactly the future behavior — 
TR of a proposed (but not yet existing) structure during a future earthquake. Sec- oe 
a 
jee to be satisfactory only for simple, , ordinary types of structures, and can- *evg 
net be expected togive good results when to structures that are unusual q 
q The following corrections should be made = the working of the author’s . 
iginal paper. The last sentence in Section 3 should be changed to read, ' ‘ror 
a such cases it has been proposed tousethe square root of the sum of the squares 
of the maxima.” The last sentence on page 126 should read, “If the design and 
= construction were such that the structure had only small damping during elastic = _ a 
; _ vibrations, the spectrum curves show that during strong ground motion the vi- - 
BF as brations of sucha structure could become worse than the equivalent 25% g load Bs 


128 should read, “ 


q 
, ; by building codes.” The fourth sentence from the end of the paper on page 128 + 


q _ should read, “One of these projects is a study of how to compute the ultimate | 
3 lateral strength of structures; another, being done at Stanford University, is a 5 
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BEHAVIOR OF BUCKLED RECTANGULAR PLATES 


Discussion by | Morris Ojalvo and Ahmet 


= MORRIS | OJALVO, 18 M. ASCE and AHMET CAKMAK, 19_ Mr. Stein has | pro- J 
an paper on the behavior of ‘flat plates in edge compression 


May 


discussion will compare the results obtained the for the in- 
_ finitely long plate with the results of a similar investigation, in which S. Le- a 7 
= method of solving the large deflection plate equations was used. 20 In 


Assoc. Prof., Dept. of Civ. Engrg., Ohio State Univ., Columbus, , Ohio, 
19 Graduate Student, Columbia Univ., New York, N. Y 
—— 0 “Bending of Rectangular Plates Under Compression, in One Direction Only, ‘With | 
_ Large Deflections,” by Ahmet Cakmak, thesis presented to Princeton Univ., Princeton Be 


_ * J., in 1958, in partial fulfilment of the requirements for the degree of Master of Science Be 
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Ed the c comparative investigation, the effective width was plotted against average 
_ unit shortening in the direction of the load for plates having length-width ratios 

between 0.4and 1.0. The effective width curves for rectangular plates of length- a 
width ratio below 0.8 do not appear inthe reference cited. 20 These were com- es 
puted after the first computations demonstrated that the envelope of the — 

5 tive width curves could give the solution for the infinitely long plate. al 


Tesults consider the first three terms of the s series for the deflections of the - 


to one of versus shortening so that it can be 
the ‘Tesults of ‘Mr. ‘Stein’ “second approximation” these figures 


4 
at and be refers to the effective width. As can be seen, there 
- is very little difference and it may be concluded that Mr. Stein’s perturbation 
method! is as accurate { as ‘Mr. Levy’s method, w which h requires the use of the dig 


- mm j Fig. 15 shows several of the effective width versus shortening curves and the = (7 
envelope curve for the infinitely long — Fig. 16 the envelope curve is 

4 
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COMMENTARY ON PLASTIC DESIGN IN STEE 


| that “It has been shown that if the plastic moment of resistance ‘is s adequate at 
the rolled section (section 1 Fig. 8.6) and at the re-entrant corner (section 2 
‘Fig. 8.6), there will be no necessity of checking at other sections along the | 
haunch.” statement, which was made when hen only bending of the haunch was 
considered, is questioned by this writer, 
_ For the haunch shown in Fig. 1, between sections 1 oe 2, “the same as- 
ne ~ sumption is made that the original paper makes, namely, that the cross section 
7 is symmetrical and that the centroid is at mid depth. _ The haunch then can be 7 
treated as a cantilever beam ( Fig. 2) fixed at A (relative to B) and acted upon 
by a concentrated load Q and a moment M at the end B. The moment M can be | - 
Clockwise or counterclockwise and and assumed positive in in in in 


— ‘The resisting plastic mor ment Mpx of the member at any cross section can 


is the yield point stress. The external moment Ma~ at 7 ae Ta 


M = Q(L +M 2) 


two curves of Fig. 2), which will be equal to p (the distance of the e plastic —_ > 
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— October, 1960, 
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pnt 


which - = substituting 
B tor = and { the of of Q given by Eq. an in 
= 


Eq. it is found that for values of 


L and that hinge develops 


q at an intermediate section between the ends A and B due to the applied mo- 
z (ah M and the corresponding Q (Eq. 4). Therefore, it is evident that even 
t bending at the ends, it may fail at any other 
section where the plastic develops. 


hough the member is safe for 


_ from the end A) it is in 
— we 
— 
— — 
— 
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jon by R R. J. P.Garden 


J. P. GARDEN.!—The report gives a a concise explanation of the 


binge method of computing deflections, and states that the other four methods 
_ provide a small increase in accuracy at the expense of of considerable inc increase : 


‘The plastic hinge ‘method achieves its practical validity because it obtains — 


iad deflections that would be correct for loading, slightly different from the stip-. 


oa: A method will be described of making truer computations of deflection with 
little increase in complexity. This is, there is little complexity in making the ~ 
actual co computation for deflection if the e designer is is provided beforehand with 

properties applicable to the shape of section and the nature of the material. z 

_ The work involved in providing the prerequisite properties is not great for [i 

_ simple shapes of section with simple s stress-strain relations; it is greater f for q 
materials with stress- strain ‘relations, such as exhibited by concrete. Apart 

7 from mention that the ‘method is applicable to other types of stress-strain re- 

= the illustrations given will be confined to plain rectangular sections of - 

material assumed to accord with the simple plastic theory. This holds for” 

Basically, the method relies onthe fact that a property, called by the writer, 
plastic curvature, can be so determined that the deflection ofa | struc- 


ture in the elastic- -plastic range is to be determined by considering the ; struc- os 
ture to remain ina sper elastic state under the influence of the actual elas- 


description of the method has been published elsewhere. 3 For the 
present discussion, it is sufficient to state two key values relating to commuted | 
curvature for rectangular sections behaving in accordance withthe simple plas-— 
tie theory. ‘Thus, for any member having imposed on it moments that rise lin- 
a early through My to Mp over a length of member called 1, the amount of com- 
mated curvature due to th spread will always be equal tos 
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curvature, 
at from the point where M = M (ges Fig. 1). As Mp = 1.5 My for a 
-— tangular section, 1 will be ae of the Line ength x on Fig. 1, and =" will be 1/27 x. 


* Fig. 3 indicates, by open circles, the locations of the centers of gravity of com- 
muted plastic curvatures, anda solid circle shows the location of the _ point 
The location of the point of concentrated hinging may be ail by the 
method of trial as used in the report, with additional terms used to take ac- 
count of t the commuted plastic curvature. To illustrate the use of the 
terms in the slope _ deflection equations, , the working for the second trial only, a : 
‘First, to illustrate how these commuted curvatures are included in the slope 
‘3 BO ron equations, . Fig. 4 (a modification of Fig. 9.3 of the report) is given. a 
The full circles indicate localized hinging 6", and 6"'p, and the open circles in- 
= dicate of Gravity of to of 
(8 
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“FIG. 2, EXTENT oF PLASTICITY | APPROPRIATE To 


— 
a 


9 43 and nd 45 and with gives: 
L2 
0. 095 in 


Mp 
H4 = 0807 (c, 0. 0303 in 
me commuted vebstions, 45 and "43, amount runt to [- 13/(12 x | x a7) - 13 
6x 27)] MpL/EI=-13MpL/108EL 


_ This spread plastic rotation is ignored in the plastic hinge en ‘it fol- ‘ 
lows that this emma, must tacitly have undervalued the positive plastic rotation 


on by and the difference in oy by tl the two methods be 


131 Mp I L108 EI (1/2) 
Mp Li p L2/432 El or 0.03 M, Mp L2/EI 
‘The actual difference of calculated values. is 0.0295 Mp L2/E1. 


ceding, the following points can be made: 

(1) The deflections derived in the report will very nearly be true_ values a 
loads that are only a little different from the stated loads. 


x _ The results from the two methods would show much less difference for Boi 

: 


— 
— 
7 : at andnear point 3 by an amount that will be very nearly equal when expressed oa mil 


this xa then, the method is good enough in practice. 


However, as referred to in Table 9.1 of the report, this simple method yields 
large errors in n curvature and deflection for members whose moment values of — 
Mp are constant along a finite length, and also for members that have, at hinge _ 
points, a curvature moment diagram where the tangent tothe curve, at the hinge 
- point, is parallel to the / member. For instance, in Fig. 9. 4 of the report, the 
curvature near point C in the rafter, is very great at full collapse load. Indeed, © we 
a rough calculation shows that, again for members of rectangular section, a 
rotation due to spread plasticity near C that would occur with Mg at only 0.9 a} ; 
a = is sufficient of itself to cause a vertical deflection at D of about the same 
3 magnitude as the quantity given in Eq. 9.17 of the report for the total deflec- = 
D for full collapse load. The following comments can be made: 


‘a (1) The simple plastic: ‘theory ee give a very great deflection of D for 
full nominal collapse load even for members of usual structural sections. a is 
(2) Occurrence of strain-hardening will reduce the discrepancy. 
(3) The deflection derived in the he report will be true for loadings appreci- 


ably, but not grossly less than full nominal collapse load. 


Constants, suchas used in this discussion, can be derived for any structural = 


1 shape and applied in the same way for any material assumed - to have simple 
plastic behavior. The effect of upper and lower. yield point can be incorporat- 
For strain-hardening and other continuously varying stress- -strain rela- 
aA eae, the parameters have values that vary with values of M, but again, similar © 
§ use can be made of constants that are appropriate to > any value of My chosen q 


to limit the desirable strain. ‘The: writer would be interested to hear the opin- 


Meat of committee members on the value of the commuted curvature method for 
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by O. C. Zienk nkiewi 


a problem for which visco- elastic cotunlone are of considerable interest. “The. 

bs general terms in which the analysis is carried out and the admirably clear in-_ 

_ troduction will serve as a starting point for other similar investigations. Itis | 
_ the object of this discussion to comment on some minor points omitted in the ea 
Peer and to suggest a possible alternative approach to this class of problem. — ‘ 
? — first of the points concerns the question of boundary conditions of the 
plate. Fora complete separation | of the variables used in Eqs. 14 through Po 


‘it is necessary for the boundary conditions to be similarly separate. This i 


which in 1 terms of of the separated va variables will vill correspond ond to 
It will be noted that a similar separation is 


_ edges of a plate in which the conditions _ me 


tt, however, a plate with a free edge is considered or a ere in which the ae 
aa is given some ge slope or deformation (which itself can be an in-- 
lo 


a question of notation. of the meaning suffices a 


lowing procedure is s at times advantageous. “ This, in 1 essence, isa a well- pteenl 


June, 1960, by George E. Mase, 
2 Prof, of Civ, Engr. Northwestern Univ., Evanston, Ill. we 
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input” of dé — at a time: T from which results 


Dre (t-7) d 
It is ‘that ) is in ‘fact t the “ creep function ” of the plate 
under constant loadand, in the particular example discussed by the author with 
 aplate of Kelvin material, corresponds with the bracketed term of Eq. 24. ee 
any particular form of 6(T) defining the load-time variation, the Duhamel in- 


This is essentially to the author’ s s Eq. 21. but has this this thats 
whereas it is often difficult to approximate the behavior of a visco-elastic ma-— 
terial bya simple mechanical model, it is a relatively simple matter to arrive 
at the creep function formula by an experiment. The integrations can then mS 
af tay _ ‘The integral approach is not limited to the “ “proportional loading ” case. It 
.- is possible to tackle with its use problems of the moving load type or, in fact, oo! 
_any in which the loading distribution is not of the simple type defined by Eq. a 
To approach this type of problem let us consider the deflections caused by an 
incremental step function of load. 


har; t= tx, yt). 


placed at time a ft, on on the plate. As this increment by itself does in fact 
_ obey the proportionality relationship, it is noted that the deflections can be ob- 
g tained by an | elastic solution of the plate with rigidity, D, loaded with a distri- 
buted loading of magnitude df (6(r) is a unit step function now), multiplied by 
The preceding statement is clearly true for all the increments of which the 
load, f, is composed, and the final solution can be synthesized by superposition, 
as the sum of such deflections. 
antes appears to require a series of solutions of the elastic problem for dif- 
= ferently distributed forms of loadat different times. If however the deflections 
_ ats some specified time, t, are ‘required, it is apparent that the same deflection 7 
hy 5 would be achieved by loading the elastic plate of rigidity D with a loading df 
: { & (t-7) at the instant, t, and that the total deflection could be obtained sim- 
— 


ply by solving the elastic to a distributed of 


or by the biharmonic| equation 


at any instant of time. ‘This latter form shows indeed that any of the 
£ known elastic techniques of selena flat plates can be easily adapted tothe treat-— 


y,t 


Of 
= nett -7) 5 


a ‘5 alent form of the loading given previously, is used. It is for instance possible 
. ¢. to apply numerical methods of solution in exactly the same way as done in elas - 
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plate deflection after the preliminary step of evaluating p has" 
out. The application of this to the moving load type of problem is evident. fe 
es” is interesting to note that in the examples cited by the author the distri- — 
bution of bending moments in the plate varies in 1 time ina ‘different manner 
fro om the assumed variation in load. In fact, ,for the case of a constant load dis- - 
tribution applied at time zero the moments aretime dependent forthe assumed | 
_ Kelvin and Maxwell types of bodies. With many visco-elastic materials the 


"compressibility is negligible while in others, such as concrete for instance, it _ 

 i-_—e that the equivalent Poisson’s ratio is constant. _ This in effect means x 
that the operators Q/P and K of Eq. 3 are of a form ‘differing by a constant ~ 

It is easy to show that in such cases with a constant load the q 


_ distribution of moments is also constant with time and, in fact, is indentical _ 
the appropriate elastic solution. 


multiplier only. 
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of Mechanics Division 


In 7 the s ‘second paragraph, paste 


173, the material in at the end of this item should 


be chan ed to the followin 


The Polytechnical School in Mexico City ane, as a result of | 7 
perusal of an article in CIVIL ‘ENGINEERING, October, 1957, p. 79, to 
have been designed by the dynamic theory. Further study reveals, how- 7 

ever, that this building had a short did that engineers 


98. Ea. 17 add a final parentheses. 


P. 93. NE Eq. 18b insert a just prior to th the first number one. 
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